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ABSTRACT 

Context. We present a combined photometric calibration of the Supernova Legacy Survey (SNLS) and the SDSS supernova survey, 
which results from a joint effort of the SDSS and the SNLS collaborations. 

Aims. Our primary motivation is to eventually sharpen cosmological constraints derived from type la supernova measurements by 
improving the accuracy of the photometric calibration. We deliver fluxes calibrated to the HST spectrophotometric star network for 
large sets of tertiary stars that cover the science fields of both surveys in all photometric bands. We also cross-calibrate directly the 
two surveys and demonstrate their consistency. 

Methods. For each survey the flat-fielding is revised based on the analysis of dithered star observations. The calibration transfer 
from the HST spectrophotometric standard stars to the multi-epoch tertiary standard star catalogs in the science fields follows three 
different paths: observations of primary standard stars with the SDSS PT telescope; observations of Landolt secondary standard stars 
with SNLS MegaCam instrument at CFHT; and direct observation of faint HST standard stars with MegaCam. In addition, the tertiary 
stars for the two surveys are cross-calibrated using dedicated MegaCam observations of stripe 82. This overlap enables the comparison 
of these three calibration paths and justifies using their combination to improve the calibration accuracy. 

Results. Flat-field corrections have improved the uniformity of each survey as demonstrated by the comparison of photometry in 
overlapping fields: the rms of the difference between the two surveys is 3 mmag in gri, 4 mmag in z and 8 mmag in u. We also find a 
remarkable agreement (better than 1%) between the SDSS and the SNLS calibration in griz. The cross-calibration and the introduction 
of direct calibration observations bring redundancy and strengthen the confidence in the resulting calibration. We conclude that the 
surveys are calibrated to the HST with a precision of about 0.4% in griz. This precision is comparable to the external uncertainty 
affecting the color of the HST primary standard stars. 

Key words. Cosmology: observations - Techniques: photometric - Methods: observational 



1. Introduction 

Substantial efforts have been spent in the last few years to im- 
prove the accuracy of the calibration of large photometric sur- 
veys (e.g. Ivezic et al. 2007; Padmanabhan et al. 2008, for SDSS 
Regnault et al. 2009, for SNLS Schlafiy et al. 2012, for Pan- 
STARRS), establishing large catalogs of stars with broadband 
magnitudes known with an accuracy of 1 or 2%. The scientific 



* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de l'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products 
available at the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative project of 
NRC and CNRS. 

** Tables 22 and El - E 10 are available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 



goals pursued cover a broad domain from stellar physics to cos- 
mology. One of the main drivers for precision photometry is the 
quest for precise cosmological constraints from the Hubble dia- 
gram of type la supernovae. 

In spite of the vast improvements in photometric calibration 
over the last decade, calibration uncertainties remain the domi- 
nant source of systematic error limiting the precision of cosmo- 
logical constraints obtained from more than 400 well-measured 
type-la supernova (SN la) light curves (see e.g. the discussion 
of systematics in Conley et al. 2011). As an illustration, the ac- 
curacy of the measurement of the dark energy equation of state 
parameter w provided in Sullivan et al. (201 1) is about 8% while 
it would reach 5.7% with a perfect photometric calibration. The 
calibration accuracy is likely to remain a serious issue in up- 
coming surveys expecting an order of magnitude increase in the 
number of supernovae. 

The photometric component of the SNLS was conducted at 
the Canada France Hawaii Telescope (CFHT) using the wide 
field MegaPrime/MegaCam camera. During its five years of sur- 
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Fig. 1. Schematic view of calibration data used in the present 
work. Each box represents a set of standard stars established in 
that photometric system. Each arrow shows the available cross- 
calibration data corresponding to the measurement of flux ratios 
between the two different sets of standards. The two new sets of 
SNLS measurements introduced in this paper are indicated with 
red dashed lines. 



vey operations (2003-2008), about 1000 multi-band SN la light 
curves were discovered in the redshift range 0.2 < z < 1, 500 
of which were spectroscopically identified. The first three years 
of the data sample are analyzed and published in Conley et al. 
(2011); Guy et al. (2010); Sullivan et al. (2011). The SDSS-II 
SN Survey was one of the three components of the SDSS-II 
project. During three three-month seasons of survey operations 
(Fall 2005-2008), the SDSS SN survey discovered multi-band 
light curves for 500 spectroscopically confirmed SNe la in the 
redshift range 0.01 < z < 0.45. The processing and analysis of 
the first year sample are given in Holtzman et al. (2008); Kessler 
et al. (2009). 

This work is part of the SDSS+SNLS joint analysis. It 
presents the collaborative effort made to improve the accuracy 
of photometric calibration. Previous calibrations for the two sur- 
veys are described in Holtzman et al. (2008) and Regnault et al. 
(2009) (hereafter R09) respectively. This paper supersedes both 
by providing a common and consistent calibration for the two 
surveys. We exploit dedicated CFHT observations that comple- 
ment the usual calibration data. The complementary data con- 
sists of direct observations of primary spectrophotometric stan- 
dard stars and of SDSS science fields. 

The advantage of these complementary observations is three- 
fold: 1) they provide a direct cross-calibration of the tertiary 
standard stars used to calibrate the supernova measurement in 
both surveys, 2) they enable an accurate consistency check of 
the photometry between both instruments, and 3) they provide 
redundant paths, subject to different systematics, to anchor the 
common cross-calibrated set of tertiary standards on the flux 
scale defined by the HST white dwarfs (Bohlin et al. 2001). The 
various calibration paths to the HST system are illustrated in 
Fig. 1, where the red dashed lines show the paths enabled by 
the complementary CFHT observations. 

Since the publication of the SNLS calibration in R09, we use 
additional calibration data including a factor of two more data to 
map the photometric response and data to calibrate the replace- 
ment i band filter that was installed in 2007. In addition, the con- 
sistency checks introduced in this work showed subtle problems 
in the flat-fielding of both surveys that required corrections. As 
a consequence, this paper also describes a reanalysis of the orig- 
inal calibration data for the two surveys. 



The paper is organized as follows. An overview of the sub- 
ject is given in section 2. The determination of the MegaCam 
flat-fielding solution, the determination of a consistent set of ter- 
tiary standards in the SNLS fields, and the calibration of the 
SNLS are detailed in sections 3, 4 and 5. Section 6 and 7 describe 
a redetermination of the photometric calibration of the SDSS su- 
pernova survey. Section 8 describes the analysis of the dedicated 
cross-calibration data that anchor the two surveys to the same 
flux scale. We gather all the relevant data to provide a combined 
calibration solution to both surveys in section 9. Finally, in sec- 
tion 10 we summarize the properties of the calibration products 
and discuss the consistency of the available calibration data. We 
conclude by identifying survey and detector design limitations 
and present perspectives for future improvements. 

2. Overview 

2.1. Definitions 

The purpose of photometric calibration is to relate the instru- 
mental fluxes measured in the image pixels to the physical fluxes 
of the observed objects. Instrumental magnitudes are defined as 
'«adu = -2.5 log 10 <Pad\j- The calibration is commonly divided 
in two steps. The first step consists in the compensation of the 
variation of the effective instrument throughput (in space and 
time) to deliver a homogeneous and consistent set of measure- 
ments. The second step consists in the delivery of a model relat- 
ing the homogenized quantities to the physical flux of the astro- 
physical objects. 

The results of these two steps are calibrated broadband mag- 
nitudes (m) of objects in a photometric system whose interpreta- 
tion in terms of physical fluxes is given by 



m = -2.5 log 10 



J A AT(A)S(A)dA 
J A AT(A)S le f(A)dA 



(1) 



S (A) is the spectral energy density (SED) of the object above the 

earth's atmosphere in units of ergs cm A , T(A) is the overall 
instrument throughput (including the atmosphere) that defines 
the natural photometric system, and S K { is a reference flux scale 
to which the object flux is compared. T{A) is expressed as the 
detector response to incident photons and its normalization is 
irrelevant because it cancels in the ratio. The factor A is required 
to convert the SED to a function proportional to the number of 
photons per unit wavelength. 

In what follows we refer to the AB magnitude system as 
defined in Fukugita et al. (1996), unless otherwise specified. 
Accordingly, S K f = S- d b where S a b is: 



S- db (A) = 10 



-48.6/2.5 



cA' 



erg- 



(2) 



with the wavelength A expressed in angstrom and c the speed of 
light expressed in A ■ s _1 . 

The model describing the physical interpretation of the in- 
strumental fluxes has two components. The first is the determi- 
nation of the effective passband corresponding to the instrument 
measurements up to a normalization factor. The second is the de- 
termination of the flux scale of the homogeneous measurements 
in each passband. The latter is obtained from a comparison to a 
flux reference, which we choose to be the HST-CALSPEC stellar 
library (see the CALSPEC website 1 and references therein). 



1 http : //www . stsci . edu/hst/observatory/cdbs/calspec . 
html 
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Our primary goal is the calibration of the SNLS and the 
SDSS supernova surveys. The main output of this work con- 
sists of calibrated AB natural magnitudes for large sets of stars 
selected in the science fields. These stars constitute in-situ pho- 
tometric references directly usable to calibrate supernova flux 
measurements. They will be referred to as the tertiary standards, 
while in the same terminology, the spectrophotometric standards 
established by the HST constitute the primary standards, and 
stars used in the calibration transfer from the HST to the science 
fields are referred to as secondary standards. 

We also deliver an estimate of uncertainties related to this 
calibration that can be propagated to the cosmological result. 
Accurately measuring the flux ratios between different photo- 
metric bands is of paramount importance for cosmological con- 
straints, as this enables the comparison of supernova luminosity 
at different redshifts. Since cosmological constraints are insensi- 
tive to a global calibration offset on the flux scale, we concentrate 
on the accuracy of the relative calibration between the different 
photometric bands in the two surveys. We do not include the un- 
certainty on the global flux scale of the calibration reference (see 
eg. Bohlin & Gilliland 2004). 



2.2. Description of Instruments and Surveys 

The SNLS covers the four low extinction fields of the CFHT 
Legacy Survey Deep component (hereafter Dl to D4). They 
were repeatedly imaged in the 4 photometric bands gM rM z'm Zm 
every 3—4 nights during dark astronomical time 2 . Photometry in 
the um band is also available for those fields but was not part of 
the supernova survey. The original im filter was broken in July 
2007 and replaced by a slightly different i M filter in October of 
the same year. We will refer to the photometric band defined 
by the replacement filter as Um- The MegaPrime instrument is 
mounted on the CFHT prime focus and hosts the MegaCam cam- 
era, a mosaic of 36 e2v 2048 x 4612 CCDs, that covers a field 
of 0.96 x 0.94 deg 2 (Boulade et al. 2003). Raw MegaCam im- 
ages are processed by the Elixir pipeline (Magnier & Cuillandre 
2004) that handles bias subtraction, flat-fielding and defringing 
in iM and zm bands. 

The SDSS-II Supernova Survey primary instrument is the 
SDSS CCD camera (Gunn et al. 1998) mounted on a dedicated 
2.5-m telescope (Gunn et al. 2006) at Apache Point Observatory 
(APO), New Mexico. The focal plane hosts a CCD array or- 
ganized as six columns of 5 CCDs. Each CCD of the col- 
umn records images in one of five broad optical bands: ugriz 
(Fukugita et al. 1996). The camera was used in time-delay-and- 
integrate (or drift scan) mode, which provides efficient sky cov- 
erage. The Supernova Survey scanned at the normal (sidereal) 
SDSS survey rate, which yielded 55 second integrated expo- 
sures in each passband. The survey covers a 300 square-degree 
region (2.5° wide over 8 hours in right ascension) that was re- 
peatedly observed (every fourth day in average) over the course 
of three three-month seasons. The region is centered on the ce- 
lestial equator and is referred to as "Stripe 82". 

Many aspects of the analysis are complicated by two in- 
strumental effects: spatial variation of the instrument passband 
response, and stray light which adversely affects common flat- 
fielding methods. Below we describe these two issues. 



2.2.1. Varying passbands 

It was realized in R09 that spatial variations of the passband 
response resulted in non-negligible color terms between pho- 
tometric measurements obtained at different positions on the 
MegaCam focal plane. The size of the color terms between the 
center and the edge of the focal plane is about 0.02 (depending 
on the band) for main sequence stars indexed by their g - i color 
(cf. Sect. 3 and Fig. 3). The major contribution to this variation 
is related to the manufacturing process of the large MegaCam 
interference filters (30 x 30 cm 2 ). The variation was found to fol- 
low an almost perfect radial pattern, with the filters being redder 
at the center of the focal plane than at their edge. The typical 
variation of the filter transmission's mean wavelength amounts 
to a few nanometers. Other expected contributors to variations 
of the effective passband are differences of quantum efficiency 
between CCDs (they were found small enough to be neglected), 
and variations in time of the spectral shape of the atmospheric 
extinction (that average out over multi-epoch measurements and 
will not be considered). 

The measurement of SDSS effective passbands are described 
in Doi et al. (2010). Small differences (at most 2 nm variations 
of the effective wavelength against the average transmission for 
ugri, 4 nm in z band) were found between the individual filters 
mounted on the different columns of the camera. The SDSS cali- 
bration strategy is such that this issue can generally be neglected 
in the calibration transfer. 

Varying passbands bring two complications to the analysis 
of photometric measurements. First, the flux interpretation of 
a broadband magnitude (Eq. 1) depends on the position of the 
measurement on the focal plane of the instrument. Defining x as 
the focal plane position and T\ x as the passpand response at x, 
magnitudes in the natural (and position dependent) photometric 
system are denoted by 



m\x = -2.51og 1() 



J A AT lx (A)S(A)dA 
j A AT\ X (A)S a b(/l)aU 



(3) 



Second, the variation of response depends on the object SED 
and therefore it is not possible to ensure the uniformity of the 
instrument response for objects with an arbitrary SED. There 
is thus a degree of freedom in the choice of the spectrum for 
which the response of the instrument is made uniform (refer to 
Sect. 3.3 for further discussion). 

2.2.2. Stray light in the wide field corrector 

In addition to the problem of the passbands varying with the po- 
sition on the focal plane, there are other position dependent op- 
tical effects that complicate the determination of the instrument 
photometric response. Parasitic reflections between filters and 
optical elements of the wide-field corrector are the most prob- 
lematic because they cannot be modeled and subtracted easily. 
These effects produce stray light that pollutes the images ob- 
tained by observation of uniform illumination sources such as 
twilight sky, dome screen, or night sky background. These tech- 
niques thus cannot be used directly to produce the instrument 
flat-fields. 3 



2 MegaPrime shares the rest of the lunation with other instruments. 
A "run" covers at least a 14 days period around the new Moon. 



3 Another effect that must be accounted for when computing the pho- 
tometric response from images of uniform illumination is the distortion 
of the plate-scale (see Fig. 4 in R09). This distortion can be easily com- 
puted from the astrometry. Other effects can also arise from the variation 
of the PSF across the focal plane. When not properly accounted for by 
the photometry method, the effective fraction of the flux captured by 



3 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 



2.3. Current state of the SNLS and SDSS calibration 

2.3.1. Absolute calibration transfer 

The original calibration strategies for both the SDSS (Stoughton 
et al. 2002) and the SNLS (Astier et al. 2006; Regnault et al. 
2009) photometric surveys follow roughly the same scheme. As 
neither survey telescope could directly observe the primary spec- 
trophotometric standard available at that time, they relied on 
third party observations, to which they anchored via a set of sec- 
ondary standard stars. 

In the SNLS, the secondary standard stars are those estab- 
lished in Landolt (1992). The SDSS secondary star network was 
built using a dedicated instrument called the photometric tele- 
scope (PT) (Tucker et al. 2006). The PT was used to determine a 
nightly photometric solution and to deliver calibrated secondary 
standards. These secondary standards are distributed in square 
patches of 0.36 sq-deg throughout the SDSS survey area. The 
supernova survey area (Stripe 82) is covered with an average of 
one PT-calibrated patch every 3 degrees in right ascension. 

These two-step calibration strategies have limitations. First, 
they accumulate the photometric and systematic uncertainties 
of two measurements. Second, and more importantly, the sec- 
ondary standard stars provide only broadband photometric ref- 
erences rather than spectrophotometric references. Since the sec- 
ondary star photometric system is different from that of the sci- 
ence system (i.e, SNe la), this difference ultimately limits the 
precision of the absolute calibration to the accuracy of photo- 
metric color transformations between those two systems (i.e. to 
some kind of interpolation between quantities integrated over 
different passbands). 

The passbands in the Landolt UBVRI system differ from the 
MegaCam um gM I'm im Zm, and the transmission functions of 
the Landolt filters actually used are not precisely known. The 
resulting color transformations are difficult to model and have 
significant uncertainties, ranging between 0.3 to 2% depending 
on the band. These uncertainties were the limiting factor in the 
previous release of the SNLS calibration. 

The strategy of the SDSS may appear more favorable as the 
PT was specifically designed to use a filter set identical to the 
2.5 m science telescope. However, the filter sets differ signifi- 
cantly primarily because of the differences in the environment 
(vacuum vs. air) as explained in Doi et al. (2010). In addition, 
repeated, in situ measurements enabled Doi et al. (2010) to de- 
tect a significant change in the w-band filter response after the 
2001 measurement. The PT filter responses, however, were not 
measured as frequently or as carefully, and we are again forced 
to rely on empirically determined color transformations to tie the 
two systems together. The problem with the determination of the 
color transformations is further complicated by errors in the PT 
flat-fielding as discussed in section 6. 

2.3.2. Survey uniformity 

The characterization of the MegaPrime response is based on twi- 
light flat-fields, corrected, at large scale, for the optical effects in- 
duced by stray light and plate-scale distortion. This photometric 
correction is measured on dithered observations of dense stel- 
lar fields following the scheme presented in R09. Although this 
measurement does not depend on assumptions about the tele- 
scope optics, we found that it is sensitive to variations in the at- 
mospheric conditions during the dithered observation sequences. 



the photometry can vary. This makes the effective throughput spatially 
variable in a way that depends on the photometry method. 



The determination of the instrument response presented in 
R09 covers the period 2003-2006, and a significant evolution of 
the response is attributable to changes in the telescope optical 
setup. An extension of this R09 study is thus required for the 
analysis of observations after 2006, in particular observations 
made with the Um replacement filter. In addition, since the de- 
termination of the photometric corrections are affected by ran- 
dom variations in the atmospheric conditions, we expect some 
improvement from including new independent observations. 

SDSS imaging data was obtained using drift scanning: each 
point on the sky was sampled by each CCD row so that the effec- 
tive response is averaged over all rows. The flat-field is therefore 
represented by a 1 -dimensional array of values holding the rel- 
ative average response between pixel columns. Such solutions 
were determined by different techniques in Ivezic et al. (2007) 
and Padmanabhan et al. (2008). The first technique relies on 
a combination of stellar locus analysis and comparison of the 
Sloan photometry to PT secondary patches. The second tech- 
nique solves for the flat-field vector by cross-correlating normal 
observations, scanned in the right ascension direction, with in- 
frequent observations that scanned the sky approximately along 
declination. 

The flat-fielding of the supernova survey currently relies on 
the PT-based solution, while the cross-scan based solution is ap- 
plied to the main survey since the eighth data release (DR8). 
Small but significant discrepancies between the two solutions 
can be found when comparing the photometry of stars in the 
Stripe 82 (see Sect. 6). 

2.4. Project for this work 

In this paper, we review, control and improve the original 
calibration paths of both surveys, and we provide an impor- 
tant cross-check by directly inter-calibrating the two surveys. 
We also describe direct observations of suitably faint primary 
standard stars (now available in the CALSPEC database) with 
MegaCam on the CFHT, providing an alternative path to anchor 
the SNLS to the HST scale. 

The various calibration paths exploited in this paper are sum- 
marized in the schematic view presented in Fig. 1. The two 
new MegaCam calibration and inter-calibration measurements 
(shown by red-dashed lines in Fig. 1) provide a redundancy in 
the calibration of both surveys to the HST standards. The orig- 
inal "Landolt" and "PT" calibration paths involve very large 
numbers of measurements accumulated over the course of the 
surveys. These calibrations are essentially free from statistical 
uncertainty, but, as already discussed, they are limited by the 
accuracy of color transformations. On the contrary, the direct 
path provides direct observation of spectrophotometric standards 
in the MegaCam photometric system and eliminates this domi- 
nant systematic for SNLS. The direct path, however, involves a 
smaller number of observations, potentially affected by system- 
atics related to the difficulty of observing bright stars. 

A general agreement between all the available measurements 
conducted with different instruments and subject to different 
sources of uncertainty would strongly support the conclusion 
that systematics are controlled and correctly estimated. In ad- 
dition, the data gathered for the direct cross-calibration of the 
SDSS and the SNLS provide high-quality photometry in both 
systems for a large number of stars. Those measurements sam- 
ple a wide range of stellar types and a large part of both focal 
planes. Such data can deliver two important by-products. First, a 
consistency assessment of the uniformity and quality of the pho- 
tometry in both surveys. Second, a precise determination of the 
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color transformations between the two systems. The latter can 
be compared with expectations to provide a sensitive check of 
the passband transmissions. 

This work is presented as follows. A new determination of 
the MegaCam photometric response extending over the entire 
survey is described in Sect. 3. It is used to build an average cat- 
alog of photometric tertiary standards in Sect. 4. On the SDSS 
side, the question of the discrepancy between the two available 
flat-field solutions is investigated in Sect. 6. Identifying evidence 
of an error in the nominal PT flat-field (appendix A), we opt for 
using the cross-scan solution which is independent of the PT. 
The stripe 82 tertiary catalog is corrected accordingly. 

Validation of the uniformity of the resulting tertiary catalogs 
is given by the analysis of the cross-calibration data (Sect. 8, see 
in particular Fig. 28 page 44). Constraints binding together the 
flux scales of both tertiary catalogs are also derived in Sect. 8. 

Measurements anchoring tertiary catalogs to the HST pri- 
mary standards are described in Sect. 5.1 for the path using 
Landolt secondary standards, in Sect. 5.2 for the direct measure- 
ments of HST standards with MegaCam, and in Sect. 7 for the 
path using the SDSS PT secondary patches. The combination of 
all available data to provide a common and improved calibration 
of the SDSS and SNLS surveys as well as the full estimate of 
the associated uncertainties is described in Sect. 9. The reader 
can skip the details of the analysis and jump directly to Sect. 10 
where we discuss the consistency between the different calibra- 
tion paths (see in particular Fig. 29) and describe the calibration 
products. 

3. MegaCam/MegaPrime instrument model 

This section is dedicated to the determination of the two key in- 
gredients characterizing the effective model of MegaPrime mea- 
surements: the effective transmission T\ X (A), and the photometric 
flat-field solution F(x) corresponding to our aperture photometry 
method. While these two subjects are usually treated separately, 
they are correlated in MegaPrime due to the continuous variation 
of the filter transmission with position. 

Figure 2 shows an overview of the data samples involved 
in the determination of the instrument model, and highlights 
the main steps of the analysis. In section 3.1, we gather inde- 
pendent characterizations of the optical components to build a 
model of the instrument passband transmissions. The resulting 
transmission functions T\ X {A) vary continuously with the dis- 
tance from the focal plane center x (hereafter referred to as the 
radius). In section 3.2 we use synthetic photometry to compute 
the small color transformations relating stellar observations ob- 
tained at different radii. We compute photometric corrections to 
the instrument response maps delivered by the usual twilight ob- 
servations in section 3.3. Those corrections are obtained from 
dithered observations of dense stellar fields under photometric 
conditions. We use the color transformations from Sect. 3.2 to 
make sure that the resulting photometric response maps F{x) 
make the instrument response constant for an (hypothetical) ob- 
ject whose SED corresponds to the AB spectrum S a j,. With this 
choice of reference SED, the relation between the instrumental 
and AB magnitudes is a single calibration constant (the AB zero 
point) that is valid across the entire focal plane. The uncertainties 
on the results are discussed in section 3.4. 

A comprehensive model of the MegaCam effective passband 
transmissions has already been released in R09. The main modi- 
fication since then is the replacement, in October 2007, of the in- 
terference filter defining band iu after its accidental destruction 
on July 20 th 2007. This reanalysis was triggered by this incident 



i Atmospheric transmission 
J above Mauna Kea 
1 (Butonetal. 2012) 



CCD quantum 
efficiency 



Filter 
transmission 




Optics & mirror 
transmission 



Twilight 
observations 




Color 
transformations 
Sect. 3.2 




Grid 
observations 



Fig. 2. Overview of the MegaPrime photometric response deter- 
mination. Boxes figure the various data sets involved in the con- 
struction of the model (dashed for external data). Ellipses figure 
the main steps of the analysis. 



and aims mainly at providing calibration data for the slightly dif- 
ferent replacement filter (hereafter i2 M ). We also propagate im- 
provements to other bands, mainly the benefits of an extended 
data set and better understanding of the instrument model. A 
summary of main differences is provided at the end of this sec- 
tion. 



3.1. Instrument passbands 
3.1.1. Passband model 

The effective transmission model for photometric band b results 
from the product of 5 components: 



T iM) = T f h (A,x)T c (A)R m (A)Ta(^)e ccd (A), 



(4) 



where T{ is the position dependent transmission of the interfer- 
ence filter, T c the transmission of the wide field corrector, R m the 
primary mirror reflectivity, T a the average atmospheric transmis- 
sion at Mauna Kea (at a typical airmass of 1.25) as measured by 
the SNFactory (Buton et al. 2012) and e cc[ j the measured quan- 
tum efficiency (QE) of the MegaCam e2v CCDs. Note that T a 
differs quite significantly from the preliminary version used in 
R09. 

The measurements of quantum efficiency for individual 
CCDs display slight variations between chips. They do not in- 
troduce noticeable modification of the effective passband trans- 
mission in any photometric band except for g M and u M . In these 
two bands, we measure a maximum shift of mean wavelength 
between the reddest and the bluest CCD of 1.5 nm and 2.3 nm, 
respectively. The use of per chip QE curves in place of the av- 
erage one published in R09 has a negligible impact on the cali- 
bration. In addition, they marginally improve the agreement be- 
tween synthetic and measured color terms between chips, which 
may indicate uncorrected weaknesses in the individual measure- 
ments. Since the per chip QE curves do not improve the cali- 
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bration, they are not used in this release of the calibration and 
instead we use the mean QE model. 4 



3.1.2. Laboratory measurements of filter transmission 

The transmissions of our interference filters under normal in- 
cidence were measured by their manufacturers at several po- 
sitions. Measurements for the original filter set, manufactured 
by SAGEM/REOSC, are published in R09. Measurements are 
available for each filter at several positions: at the center of the 
filter, and at 23, 47, 70, 93, 117, 140, 163, 186 and 210 mm 
away from the center. Transmission curves are linearly interpo- 
lated between measurements as a function of radius to provide 
the continuous model function, Ti (A, x) . 

For the Hm filter, measurements of the transmission provided 
by the manufacturer, BARR associates, are available at the cen- 
ter of the filter and at 20, 40, 60, 80, 100, 120, 140, 160, 180 mm 
away from the center. The potential existence of leaks outside the 
main passband is not excluded by these measurements that cover 
only the expected wavelength range of the main passband. The 
180 mm measurements were repeated for each corner of the fil- 
ter. The mean wavelength varies along the radius from 761.2 nm 
on one side of the filter to 764.9 nm close to the center, which 
makes it more uniform than the original filter. There are indica- 
tions that the filter transmission function may not follow a ra- 
dial pattern. Measurements at the 4 corners display small dis- 
crepancies, with variations between the two extreme sides of the 
filter reaching 1 .7 nm. However, without complementary mea- 
surements at intermediate radii, it was not possible to build an 
alternative to a radially variable filter model. Deviations from 
the radial symmetry are barely noticeable in the observational 
data (cf. Sect. 3.4). 

The normal incidence measurements are blue-shifted to ac- 
count for the average incident angle of the CFHT f/4 beam. As 
in R09, this shift is done using the approximate relation, 



T f (A) * t] 



^1- (sin 2 0)ln\ 



(5) 



where 8 is the incidence angle, and n is an effective index of 
refraction for the filter. 

The transmission curves of the original filters were measured 
again at CFHT in 2006 at 5 angles ranging from to 8° at a 
few locations. This allowed us to check that the above formula 
applies and to determine the refractive index n for each filter. 
We then integrated the transmission function over the telescope 
beam for each location in the focal plane. 

As the replacement filter i2 M transmission was never mea- 
sured under non-normal incidence, the exact value of its refrac- 
tive index n cannot be directly determined. We used the value 
measured on the original i M filter, n = 1.6. The validity of this 
choice is indirectly assessed in Sect. 8.3. 

Last, the comparison of the manufacturer transmission 
curves with the 2006 measurement at CFHT for the original filter 
set is provided in Appendix B. It revealed an important discrep- 
ancy (an 8 nm shift of the red cut-off) in Im and ru , while mea- 
surements in the other bands were found in nearly perfect agree- 
ment. As discussed in Appendix B, we believe that the CFHT 
measurements provide a better description of the survey filters 
and we therefore corrected the iu and ru transmission curves 



accordingly. The validity of this correction could not be directly 
tested on the whole spatial extension of the filters. An indirect 
assessment using stellar locus is provided in Sect. 3.4.2 and 8.3. 

The resulting transmission curves are given in table E.3-E.8 
at several radii away from the center. A sufficiently accurate 
mapping between object coordinates and position on the filter 
is obtained by considering the filter close to the focal plane and 
using the focal distance of the instrument, F = 14.89 m. 

3.2. Stellar measurements across the focal plane 

While slightly variable, the MegaCam natural photometric sys- 
tem remains sufficiently similar between different focal plane 
positions that precise transformations between stellar measure- 
ments can be derived for a large range of the stellar population. 
Passband transmission models have been used to compute syn- 
thetic values of the color transformations between the natural 
system at a given radius and the uniform system defined as the 
natural system at the center of the focal plane (denoted xq). The 
relations take the form 



m\ x 



m\ x - 5k(x)(c\ Xo - c (x)) 



(6) 



where c\ Xo is the color index of the star, c° is the reference color 
for the linear transformation, and 5k is the color term adjusted 
to a given stellar population. 5 The choice of c\ Xo is arbitrary, and 
we choose (g - i)\ Xo because it is well measured for most stars 
and it is a good proxy of the star temperature. 

The primary use of these transformations lies in the process 
of determining the photometric flat-field solution {cf. Sect. 3.3) 
for MegaCam, using measurements calibrated on Landolt stars. 
Since this calibration procedure delivers magnitudes in an ap- 
proximate Vega system, we determine transformations in that 
system as described in R09, and we make use of the primary 
spectrophotometric standard BD +17 4708 as a flux scale refer- 
ence. More specifically, in this system the synthetic magnitudes 
are defined as 



m 



-2.5 log 



id 



J A AT lx (A)S(A)dA 
J A AT lx (A)Sl f (A)dA 



(7) 



where 



SL(X) = S bd n(A)10 0Am >*" : 



4 http : / /www . cfht . hawaii . edu/Instruments/Imaging/ 
MegaPrime/ data . MegaPrime 



with S bd n the SED of the star BD +17 4708 6 and m bd n the 
MegaCam magnitude of this star inferred from Landolt observa- 
tions (see Sect. 5.1). As those scaling constants {mbdn) are close 
to the Vega magnitudes of BD +17 4708, this system is by con- 
struction an approximate Vega system. To avoid confusion with 
the AB magnitudes that we use everywhere else in this paper, 
magnitudes in this system will be, when needed, distinguished 
as mj x where the V superscript stands for Vega. 

The relations defined by Eq. (6) were adjusted in a limited 
color range on synthetic magnitudes computed from the Pickles 
(1998) and Gunn & Stryker (1983) stellar libraries. As an illus- 
tration, the adjustment of synthetic transformations between the 
center and the edge of the focal plane are shown on Fig. 3. The 
selected color range corresponds to a region where the linear ap- 
proximation is sufficiently accurate. 

5 Shall the selected standard S lc f be the spectrum of a common star, 
one would expect by construction that c? » 0. One may notice it is not 
strictly the case neither for a Vega or an AB system. 

6 The version used in this paper can be found at ftp : / / ftp . stsci . 
edu/cdbs/current_calspec/bd_17d47Q8_stisnic_QQ3 .ascii 
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Fig. 3. Synthetic color transformation between the center and the 
corner of the filters (17cm from the center) in each band. Dots 
are synthetic magnitudes computed for stellar spectra from the 
Gunn & Stryker (1983) spectral library. Linear color transforma- 
tion (solid line) are adjusted in the range delimited by vertical 
dashed lines. The color index is a Vega MegaCam color. The red 
square figures the synthesized transformation of the spectropho- 
tometric standard star BD +17 4708 computed from its STIS 
spectrum available from the CALSPEC database. 

The uncertainty on the slope of the transformation, evaluated 
from the difference obtained between the two libraries, is smaller 
than 2 mmag per color unit. The dispersion around the fitted 
transformation is typically smaller than 3 mmag. 

We also evaluated the sensitivity of the fitted transformation 
to extinction by interstellar dust. Spectra have been reddened ac- 
cording to a Cardelli et al. (1989) law to account for the mean ex- 
tinction of the stellar population. The extinction vector is nearly 
co-linear to the stellar locus in the selected color-range. For a 
value of E(B - V) = 0.4, the modification of the slope is typi- 
cally negligible in all bands except z, where it reaches 3 mmag 
per color unit. We conclude that, in fields with relatively low 
extinction, it is safe to translate measurements of stars lying in 
the selected color range across the focal plane using the fitted 
transformation. 7 



3.3. Determination of the photometric flat-field 

A common way to measure the instrumental response is the ob- 
servation of the twilight sky, assumed to be a uniform source 
of light. A first complication preventing the direct use of this 
technique is the variation of filter passband transmission across 
the Megacam focal plane. This effect makes the response vari- 
ation slightly dependent on the SED of the observed objects. 
Consequently, the strict uniformity of the response can only be 
ensured for a single class of objects sharing the same SED shape. 



Our first concern here is that the twilight images are functions of 
the twilight color which is not a useful reference for stellar pho- 
tometry. 

Putting aside this complication, the twilight images do not 
deliver satisfying results for the photometric uniformity of point- 
like objects, and is due to differences in the optical response to 
isotropic and parallel illuminations (see e.g. Manfroid 1995). In 
MegaCam, two major contributions are identified (Magnier & 
Cuillandre 2004): 1) distortions of the plate-scale and 2) par- 
asitic reflections on the various components of the wide field 
corrector. The combination of both effects leads to differences 
of around 10% between stellar fluxes measured at the edge or 
center of flat-fielded images. 

It is, however, expected that those purely instrumental effects 
are smooth and fairly stable in time, so that the determination of 
the photometric flat-field can take the form of a large-scale sta- 
ble correction to the twilight flat-field. The twilight flat-field thus 
corrects for small-scales and time variations of the instrument re- 
sponse, while the response at large scales is to be determined by 
other means. We will refer to this correction as the photometric 
flat-field correction in what follows. 

Smaller deviations (around 1%) also arise from variations of 
the PSF across the field that may not be correctly handled by the 
photometry method. While small, such effects make the flat-field 
definition ultimately dependent on the photometry method. We 
choose to measure fluxes in twilight flat-fielded images, 8 lead- 
ing to instrumental magnitudes written m' ADl] , where the t super- 
script denotes the use of twilight flat-field. Those magnitudes are 
related to natural magnitudes by 



OT ADU = m \* + 5Z ( X ) - Z 



(8) 



7 The parameters of the linear color transformation and their depen- 
dence on the focal plane position are available in electronic form at 
http : //supernovae . in2p3 . fr/snls_sdss/. 



where 6Z(x) is the photometric flat-field correction, and Z is a 
per-exposure global calibration offset (the zero point). Due to the 
variation of passbands the correction depends on the flux refer- 
ence chosen: roughly speaking, the color of the twilight spec- 
trum relative to the reference spectrum. 

3.3.1. The 'grid' observations 

A comprehensive model of this correction could in principle be 
built from the knowledge of optical distortions, reflections and 
of the spectral shape of the twilight. The alternative strategy that 
we follow in this work is to measure 5Z from a set of specifically 
designed observations (similarly to what was done in R09). To 
implement this strategy, the CFHT has gathered sequences of 13 
dithered observations of a dense stellar field in every passband. 
Since the beginning of MegaCam observations, 2 stellar fields 
have been observed in this manner on a yearly basis and also 
after each major intervention on the telescope optics. 

The characteristics of those 2 dense stellar fields are given in 
table 1 . Some attention has been paid to potential effects related 
to Galactic dust extinction in these low galactic latitude fields. 
The analysis of the average color locus of the stellar population 
showed that the typical value of the dust extinction corresponds 
to E(B - V) = 0.1 in the summer field (hereafter grid-2), and 
E(B - V) = 0.4 in the winter field (hereafter grid-1), with notice- 
able variations. 

A summary of those observations, hereafter the "grid" ob- 
servations, as well as associated telescope events is given in ta- 
ble A. 1. The dithering pattern followed by the observations is de- 
scribed in table A.2. It consists roughly of 13 exposures spaced 



8 The Elixir pipeline can apply photometric corrections to the images. 
When needed, we suppress these corrections at the catalog level. 
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Table 1. Dense stellar fields observed to model the photometric 
flat-field correction 



Name 


RA 


DEC 


E(B-V) a 








(mag) 


Grid- 1 (winter field) 


6h30m00 


14d20m00 


0.45 ± 0.08 


Grid-2 (summer field) 


21h00m00 


lOdOOmOO 


0.12 ±0.03 



Notes. 

(<I) Average value of the dust extinction derived from the analysis of the 
stellar locus displacement, assuming extinction follows Cardelli et al. 
(1989) law with R v = 3.1. 



by logarithmically increasing step sizes in right-ascension and in 
declination. 

3.3.2. Constraining a model of the instrumental response 
with grid observations 

An appropriate model for the twilight flat-fielded observations 
(m' ADl] ) is given by Eq. (8). The natural magnitudes m\ x of the 
grid stars are not known and have to be determined from the ob- 
servations themselves, up to a global offset. Also, due to the vari- 
ation of the effective passbands, the natural magnitude of a star 
is not expected to be constant across the focal plane. The only 
tractable way to obtain a closed system is to rely on approximate 
color transformation. Unlike what has been done in R09, we do 
not fit for color transformations and flat-field corrections at the 
same time. We choose to rely on the linear transformation from 
section 3.2. For stars lying in the validity range of those transfor- 
mations, Eq. (8) translates to the linear model already described 
in R09, 

m ADU - Sk{x)[g\ Xo ~ i\x - c°(x)] = m lxo + 6Z(x) - Z . (9) 

Free parameters in this model are the reference instrumental 
magnitudes of all stars at the center of the focal plane m\ Xo , the 
zero points of each exposure in the sequence Z, and the shape of 
6Z{x) that has to be parameterized. We also need the reference 
g\x a - i\x colors of the stars. It is however sufficient to have this 
color determined only approximately, because it appears multi- 
plied by 5k which is typically small (6k ^ 0.02). The construc- 
tion of an approximate catalog for the grid-field stars, neglecting 
in a first pass the impact of filter changes on the flat-fielding, is 
the subject of section 3.3.3. 

An irrelevant constant can be exchanged between the 3 terms 
of the right-end side of Eq. (9). We therefore assume 6Z — 
at the center of the focal plane and Z = for the first expo- 
sure. However, this system remains poorly constrained because 
of a partial degeneracy between the flat-field correction and the 
variation of Z between dithered exposures. The degeneracy is- 
sue making the construction of photometric flat-field (or "super- 
flats") for a CCD camera difficult is investigated in Manfroid 
(1996). A simple solution is to make the assumption that the 
whole sequence is photometric, i.e. to assume that Z is constant 
during the sequence. This assumption is hard to satisfy exper- 
imentally, even for the clearest nights, mainly because of vari- 
ations in the PSF between exposures inducing variations in the 
aperture corrections. 

To account for part of the aperture correction variations we 
used aperture photometry with reasonably large radii scaled with 
the image quality (IQ). The IQ is computed as the average sigma 



of a Gaussian fit to stars on each CCD. 9 We then perform the 
photometry in apertures of radius 7.5 x IQ. 

Second order variations of aperture corrections, as well as 
variations of the atmospheric extinction between successive ex- 
posures are to be expected at some level so that only a part of 
the zero point variation is corrected. Still, we rely on the as- 
sumption of photometricity during the observation sequence to 
alleviate the degeneracy problem and we expect that the result- 
ing error will average out with several determinations obtained 
from uncorrelated sequences. This issue is further discussed in 
section 3.4. 

The exact shape of the 6Z correction is unknown and thus has 
to be parameterized. The scheme adopted in R09 was to subdi- 
vide the focal plane into 1296 cells of 512 x512 pixels on which 
5Z is assumed constant in each cell, leading to a pixelized rep- 
resentation of the correction. Here we build a nearly equivalent 
model with a smaller number of parameters from the following 
considerations: 1) optical effects are expected to be a smooth 
continuous function of the focal plane position, 2) discontinu- 
ities between CCDs are introduced by the IQ-scaled aperture 
photometry method because different aperture radii are used to 
process different CCDs on the same image, 3) variations of quan- 
tum efficiency shape from one chip to another may introduce dis- 
continuities between CCDs if the twilight color is significantly 
different from that of the standard reference and 4) occasional 
jumps of the amplifier gains were noted in R09, and have to be 
accounted for. As a consequence we choose to expand 5Z as 
the sum of a smooth function decomposed on a uniform cubic 
B-Spline basis plus a discrete step function 5G accounting for 
differences between amplifiers: 

N 

8Z(x) = Y J b i B u3 (x) + SG a , (10) 

where a indexes the 72 amplifiers. We impose separately SG26 = 
for the amplifier 26 that lies close to the center of the focal 
plane and £ (=0 bjB^Xo) = 0. The choice of N, which sets the 
number of model parameters, is discussed in Sect. 3.3.4. 

Measurements are mostly affected by photon noise, which is 
approximately modeled as a Gaussian process whose variance 
is given by the flux integrated in the aperture. Only measure- 
ments with a signal to noise ratio exceeding 10 are kept and a 
weighted least-square minimization is used to estimate model 
parameters. Although the stellar instrumental reference magni- 
tudes contribute a very large number of nuisance parameters, the 
problem remains tractable if its sparsity is correctly taken into 
account. 10 

3.3.3. Construction of grid field catalogs 

The transformation of instrumental magnitudes to the reference 
system requires the approximate knowledge of the star color in 
the reference system. On a first pass, a calibrated catalog of the 
grid fields is built from all the gathered observations of those 
fields. These observations comprise a greater number of epochs 
than the grid dataset alone as those fields were routinely ob- 
served to measure the PSF of the CFHT instrument. The latter 
set of observations is usually slightly deeper than the nominal 
grid observations. 

9 More precisely: IQ = •y'det(M), where (M) is the average matrix 
of second moments for stars on the CCD. 

10 We use the sparse Cholesky factorization routine from the 
CHOLMOD library (Chen et al. 2008) to solve for the normal equa- 
tion. 
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Grid-1 



2.0 




Fig. 4. Stellar population within the two grid fields. The green 
dashed line is a synthetic stellar locus computed from the Pickles 
stellar library. Spectra have been reddened to match the average 
extinction measured in the field. The black arrow shows the cor- 
responding extinction vector. 



The flat-fielding of the observations is made using an approx- 
imate version of the flat-field correction 5Z 8 ignoring the effect 
of varying passbands. The simplified model is 



= m lxo + 5Z g (x) 



(11) 



Besides the approximate flat-fielding, the averaging procedure 
adopted to build grid field catalogs is identical to the one used to 
establish the tertiary star catalogs, and is extensively described in 
section 4. Once obtained, the average catalogs are calibrated to 
the Landolt system using same-night observations of the Landolt 
secondary patches as described in section 5.1. The catalogs are 
cleaned a posteriori to keep measurements with error below 5% 
and with reduced y 1 < 2 (which excludes some variable stars). 
The color-color diagrams of the resulting catalogs are shown in 
Fig. 4. 

Our primary interest lies in the measurement of the g M - i M 
color of the stars. This measurement is affected by the uncer- 
tainty on the gM and iy calibration. From the dispersion between 
the different calibration epochs, we estimate that the global un- 
certainty on gM — iM is below 0.008 and 0.004 for the grid-1 
and grid-2 catalogs, respectively. Another error is induced by 
the approximate flat-field applied to the observations and by the 
difference between the natural color g\ x - i\ x that we measure and 
the uniform color g\ Xa - i\ Xo that we use to index color transfor- 
mations. For stars in the color range considered the maximum 



difference between g\ x - i\ x and , 



'1*0 



is expected to be 0.01. 



Finally, the noise on the gM - Im color is expected to reach at 
most 0.07 rms for the faintest stars in the sample. As the 5k color 
term between the natural and uniform system is at most 0.027, 
all error terms have a negligible impact on the translation of star 
magnitudes to the reference system (at most 2 mmag rms for the 
noisiest stars). We will hence neglect the errors induced by the 
color transformation in what follows. 



3.3.4. Results 

The color transformations determined in section 3.2 are readily 
applicable to colors calibrated on the Landolt system. Grid mea- 
surements are matched to the color catalogs described above, 
and we use measurements of stars lying in the color range where 
the linear color transformation 5k is valid. The model (9) is then 
fit to the selected data. 

In all bands, we have computed solutions with increasing 
resolution of the spline basis. Increasing the resolution above 
N = 12 x 12 splines for the focal plane does not bring a signifi- 
cant^ 2 improvement nor significant changes to the solution. We 
thus selected this resolution to compute all flat-field corrections. 
The systematic error on the flat-field induced by the inability of 
the chosen parameterization to describe the actual shape of the 
photometric correction is below 0.001 mag. 

The main advantage of the parameterization adopted here is 
to decrease the number of parameters needed to describe the 
flat-field correction. The typical number of parameters for the 
selected resolution of the spline basis is N amv \i + A^ sp ii ne = 214, 
compared to 1295 in the scheme adopted in R09. The statisti- 
cal uncertainty on the reconstruction of the flat-fielding solution 
is below 1 mmag, compared to ~ 4 mmag in the previous R09 
scheme. 

We checked the sensitivity to outliers using simulations with 
2% photometric outliers 11 uniformly distributed in the range 
[-2; +2] mag. The simulation with and without outliers gives the 
same result, showing the efficiency of the procedure to discard 
outliers. 

As already stated, the model does not include varying photo- 
metric conditions during the observation of the grid sequence. 
Unfortunately, while high quality photometric time has been 
dedicated to the grid observations, very few sequences were 
found to exhibit a stable IQ. The flat-field solutions from the 
grid sequences with the most stable IQ are displayed on Fig. 5. 

The error induced by variation of the photometric zero point 
between exposures is expected to dominate the error on the de- 
termination of the flat-field correction. The selection of data to 
minimize the impact of zero point variations is discussed in the 
next section. 

Finally, the flat-fielding correction that satisfies Eq. (8) de- 
pends on the flux reference of the magnitude system. The cor- 
rections derived so far refer to the system defined by Eq. (7). By 
construction, those corrections make the natural magnitude of 
BD +17 4708 independent of focal plane position. Corrections 
needed to ensure the same property in an AB system are read- 
ily obtained (up to an irrelevant constant that can be chosen to 
ensure 5Z ab (xo) - 0) as: 

f AT [x (A)S ab (A)dA 
5Z a \x) = 5Z\x) - 2.5 log 10 J r J ' _ - (12) 



/ AT lx (A)Sl f (A)dA 



Once the corrections are obtained, photometrically flat- 
fielded instrumental magnitudes are given by: 



m A DU = m' ADV - 5Z 



ab 



(13) 



Since the correction is bound to a photometry method (IQ- 
scaled here) and a photometric system (AB in the rest of this 
paper), it is generally advisable to keep the image pixels scaled 
by the twilight flat-field and to apply needed corrections to cata- 
logs. 



11 The fraction of outliers introduced in the simulation is similar to 
what is observed in zm observations (1.7%) and significantly higher than 
the fraction of discarded measurements in the other bands. 
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Fig. 5. Maps of photometric corrections (relative to the twilight flat-field). Each panel displays a low resolution map (36 x 36 super- 
pixels of 256 x 256 real pixels, i.e. a resolution of about 2') of the 5Z corrections. Those corrections are determined on dithered 
observations of dense stellar fields in a given photometric band. Observations for the corrections displayed in um gM ?m im an d Zm 
panels have been obtained during the same night in 2005 under remarkably stable atmospheric conditions. The Um panel displays 
the result of an observation taken in 2008, about 1 year after the replacement of the filter. CCD boundaries are figured by white 
dotted lines. 



3.4. Discussion 

We now turn to a discussion of uncertainties associated with our 
determination of the instrument response. Noise does not con- 
stitute a limiting factor in any of the measurements required to 
build the response model. Instead, the most significant errors are 
expected to arise from imperfect verification of the assumptions 
that the model is built on. Explicitly, we have made the following 
hypotheses: 



The transmission of the instrument as a function of wave- 
length does not evolve with time. 

The transmission of the instrument as a function of wave- 
length depends solely on the focal plane radius, and this de- 
pendency is correctly mapped. 

The photometric response can be obtained as a stable photo- 
metric correction to the twilight flat-fields. 
The observations of dense stellar fields that constrain the 
photometric correction are obtained under stable photomet- 
ric conditions. 



In what follows, we quantitatively evaluate the validity of each 
assumption. We also discuss the resulting errors on the photo- 
metric response map. 



3.4.1. Filter aging 

The transmission model for the original Sagem/REOSC filters 
was based on early measurements by the manufacturer, prior to 
delivery in 2002. Later measurements of these filters have been 
conducted at CFHT in 2006 and display perfect agreement with 
the manufacturer findings in the u M , g M and z M bands. However, 
quite significant discrepancies between different measurements 
of the tm and iu filters were found, and thus a potential aging 
effect of these filters has been investigated (cf. Sect. 3.1.2). 

To investigate this effect, we formed independent catalogs of 
the D3 field on a yearly basis, and looked for time-dependent 
color terms between the partial and the reference catalogs. The 
computed color terms are displayed in Fig. 6. The possible varia- 
tion in filter response is less than 0.001 per unit of gM - im color 
and it can be safely neglected. This test excludes a significant 
continuous variation of the filters over the course of the survey. 
In particular, this negligible change in the filters cannot explain 
the discrepancy between the 2002 and 2006 measurements of 
the transmission curves; if these measurements corresponded to 
a real change, we would have seen a color term of about 0.005 
mag per unit of g - i. The filter measurements could both be cor- 
rect, however, if the filter response changed after the first mea- 
surements were taken but before the survey began. 

We conclude from this study that the passbands did not sig- 
nificantly evolve during the survey and that the 2006 measure- 
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Fig. 6. Color terms between tertiary catalogs averaged over dif- 
ferent periods. The survey observations in the D3 field have been 
divided into 5 periods covering about a year. Each data point is 
the fit of a linear color term between the catalog obtained from a 
single year and the catalog obtained from the whole survey. The 
stellar population is indexed by the g - i color. 



ments should be representative of the survey filters. This justifies 
our decision to correct the 2002 measurements so that they now 
match the 2006 measurements at the few positions where the 
latter are available (see appendix B). 

We recall, however, that this choice does not make the filter 
model entirely consistent. As an example, some spatial discrep- 
ancies between the predicted and measured color-terms for the 
rM filter are presented in the next section and are visible in Fig. 7. 
This discrepancy indicates that the model for this filter is still not 
fully accurate, which has consequences on the absolute calibra- 
tion (cf. Sect. 9) and on the response map as discussed below. 

3.4.2. Consistency of the model of passband variations 

A first argument in favor of the adequacy of our varying pass- 
band model is that the introduction of the color term in the model 
of grid observations given by Eq. (9) indeed improves the fit to 
the grid data set without the addition of free parameters. To go 
one step further and be more quantitative, one can leave 5k as 
a free parameter in the adjustment to the grid dataset. This was 
done in R09 using a coarse parameterization of the 6k variation 
in cells of 1536 x 1024 pixels. The 6k maps have been redeter- 
mined independently on all available grid epochs. The measured 
6k are compared to synthetic ones as a function of distance to 
the focal plane center in Fig. 7. 

We have displayed the results from the two grid fields sep- 
arately on this figure to distinguish potential effects related to 
reddening by dust. On can note that while the two grid fields suf- 
fer different dust extinction, the synthetic color terms, account- 
ing for extinction, are almost indistinguishable. The fitted color 
terms agree as well, and there is no reason to question the results 
from the dusty grid-1 field. 

The agreement between the model and the observations is 
generally satisfying. It deteriorates a bit in band r M at large ra- 
dius, where the synthetic color-term becomes the most depen- 
dent on the selected color range. The observations also suggest 
that the i2 M filter may be a bit more variable than measured. 

Quantitatively, the standard deviations of the residuals of 
the measured to the synthetic 6k maps are 0.014, 0.004, 0.003, 
0.003, 0.004, 0.002 in band u M , g M , r M , i M , i2 M and zm respec- 
tively. Residuals maps are shown on Fig. 8. The most prominent 
features are the radial deviations in band and Hm already no- 
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Fig. 7. Difference between measured and modeled color terms 
(relating the natural and reference photometric system) as a 
function of radius. Solid points show the average value of the 
6k map in bins of radius. We plot separately the value average 
on all epochs of the winter field (gray dots) and summer field 
(black squares). The error bars are estimated from the dispersion 
between epochs. Data points are not independent, and the error 
is dominated by the uncertainty on the overall normalization of 
the curve due to the constraint 6k(0) = 0. The synthetic color 
terms are computed on reddened spectra from the Pickles stellar 
library, to reflect the effect of the extinction measured in the two 
fields. 



ticed on the radially averaged plot. Otherwise the residuals are 
consistent with noise in the observations. The boundaries be- 
tween different chips are also visible in band u and g, which can 
be explained by slight differences in quantum efficiency curves 
from one CCD to another. 12 

Although statistically significant, the discrepancy between 
the calculations and measurements are negligible for our pur- 
poses. We stick to the modeled color terms as they provide some 
consistency with the flux interpretation of the magnitudes, and 
significantly reduce the noise of flat-field. 

One can also note that the examination of the various epochs 
independently exclude significant evolution of the pattern with 
time. The test is however insensitive to global variations of the 
mean passband with time which has been investigated above us- 
ing tertiary stars. 

3.4.3. Stability of the photometric flat-field corrections 

Twilight flat-fields are assembled independently for each run of 
MegaCam (every lunation). In contrast, the photometric correc- 
tion is determined about twice a year from a set of dedicated 



12 Although measurements of individual chips quantum efficiency are 
available, they are known to be affected by problems such as ice con- 
densation in the cryostat. Using them in the model does not improve 
the agreement with measurements and occasionally deteriorates it. As a 
consequence, we stick to the average QE model. 
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Fig. 8. Maps of residuals between the predicted and measured 
color term to the reference system. The synthetic color terms are 
computed on reddened spectra from the Pickles stellar library, 
to reflect the effect of the extinction measured in field grid-1. 
The measured 6k are independently adjusted on each grid se- 
quence on field grid-1 and averaged over all epochs to increase 
the signal-to-noise ratio. 



observations. We now want to assess the hypothesis that this cor- 
rection is appropriate to other runs than the one it was derived 
for. 

There are three main issues to consider. First, whether the 
differential response between point-like and uniform illumina- 
tion (what is mapped by the photometric correction) is stable. 
Second, whether the successive twilight flat-fields consistently 
map the instrument response to uniform illumination. Third, 
whether the variations of the instrument response occur on time 
scales smaller than a MegaCam run Below we qualitatively dis- 
cuss each of these issues. 

Changes in the instrumental setup : The origin of a difference 
between instrument responses to point-like and uniform illumi- 
nation is attributed to three main effects: 1) the spatial variation 
of the plate-scale, 2) the pollution of twilight images by stray 
light, and 3) the distortion of the PSF that varies the fraction of 
flux recovered by the photometry. Each of these effects is related 
to the optical configuration of the instrument and should be fairly 
stable. 

However, several important adjustments to the telescope 
setup have been made during the survey, and these may have 
modified its response. We have investigated the stability of 
the instrument response using the repeated observations of the 
SNLS deep fields. The compatibility between the observations 
conducted at different epochs provide a stringent monitoring of 
changes in the survey uniformity. The analysis is fully described 
in Appendix A, and shows that telescope modifications indeed 
affect the photometry. A single photometric correction is thus 
not suitable to flat-field the entire survey. 

We identified three telescope events with a noticeable impact 
on the flat-field corrections. The first is the flip of the L3 lens 
in the image corrector that produced a noticeable improvement 
of the uniformity of the image quality across the focal plane. 13 
The second is the unfortunate accumulation of metal shavings 
on the top optical surface below the filter mechanism until its 
identification and cleaning at the end of the year 2006. The last 



13 See http : //www . cfht . hawaii . edu/News/Pro j ects/MPIQ/. 



is the accidental breaking of the iu filter in July 2007. These 
events define roughly 5 epochs: the beginning of the survey with 
the original telescope setup, the improved setup, the period of 
noticeable dust accumulation, the short period between the dust 
cleaning and the filter breakout, and the end of the survey. The 
exact dates of the telescope events are given in table A. 1 . 

Grid observations were accumulated sufficiently often to ac- 
count for each of these major events. We can determine a photo- 
metric correction suitable for each epoch of the survey. However, 
this scheme can account only for abrupt changes in the telescope 
setup. The accumulation of metal dust is harder to handle be- 
cause this effect is likely to change slowly with time but only 
a single grid observation was taken during this period (shortly 
before the cleaning when the effect was at its maximum). That 
being said, we can use the monitoring of the instrument response 
provided by deep-field observations (see Appendix A) to argue 
that, once corrected with the available grids, the instrument re- 
sponse appears compatible across all epochs at the 5 mmag level 
and does not display a significantly stronger dispersion during 
the period of dust accumulation than during the other periods. 

Variations in the twilight flat-field : Aside from the twilight 
noise, which is negligible at large scales, there is one expected 
contribution to systematic differences between different runs: 
the twilight images contain a color-dependent term as noted in 
Sect. 3.3. If the effective color of the observed twilight varies 
significantly from one run to another, the photometric correction 
would have to be modified to compensate the variations induced 
on the twilight flat-field. 

The procedure of taking twilight flats is automated so we ex- 
pect the average twilight colors to be rather stable. We checked 
this by studying twilight flat-fields in the time period 04 Sep. 
2007 through 05 Sep. 2008 when no intervention on the tele- 
scope optics occurred. We computed the ratio of each flat-field 
image obtained during this time period to the average flat-field 
image over the period. Defining an effective twilight color as the 
correlation coefficient between the twilight flat and the 5k maps, 
we have measured variations of colors between successive flat- 
fields with standard deviations of 0.04, 0.15, 0.1 1, 0.15 and 0.45 
in units of the 6k maps for u M g M r M i M and z M bands respec- 
tively. 

The fraction of those variations unrelated to changes in the 
twilight color (but simply to actual changes in the instrument 
response following roughly a similar pattern) is difficult to esti- 
mate. Nevertheless, those variations can be used to set an upper 
bound on the effect: in all cases a variation of the twilight color 
by the quoted amount contributes a flat-field error that reaches 
at most 5 mmag. The effect is small and will be smoothed out 
by averaging different photometric corrections; hence it can be 
safely neglected. 

Variations of the instrument response on short time scale: 
Last, variations of the amplifier gain arise on time-scales shorter 
than the run and cannot be properly accounted for by the twilight 
flat-field. A prominent example of an amplifier gain variation is 
clearly visible on the 2005-10-09 grid displayed in Fig. 5. One 
amplifier (the left side) on the last CCD of the second row pops 
out. This feature is present in all bands for this epoch, but since 
it is not reproduced on other grid epochs we interpret this effect 
as a variation of the gain between the acquisition of twilight and 
the observation of the grid. 

Such features normally contribute a small extra noise that 
averages out on several lunations. However, when an ampli- 
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fier jump alter a photometric flat-field correction, all the im- 
ages flat-fielded using this correction are coherently affected. 
The jumps of the amplifier gains are estimated to be generally 
small from a preliminary analysis of a monitoring experiment 
of the MegaCam gains conducted with the SNDICE experiment 
(Guyonnet et al. in prep.). 14 As such deviations are expected to 
randomly affect different determinations of the photometric cor- 
rection, a practical way to handle this difficulty is again to rely 
on the averaging of several determinations. 

Compatibility between different epochs: From all the consid- 
erations above we conclude that the assumption of a stable pho- 
tometric correction is mostly valid, at the 5 mmag level, over 
large periods during the survey. We adopt the subdivision of 
the survey in five periods, as defined above, and we attribute a 
range of validity to photometric corrections determined during a 
given period that extends over the whole period. Remaining de- 
viations, due either to the instrument response evolving slowly 
within each period, to variations in the twilight flat-field, or to 
variations of amplifier gains are difficult to characterize with a 
high level of confidence. A reliable estimate is provided by the 
fact that the photometric corrections bring the different survey 
epochs in agreement at the 5 mmag level. An important point is 
that the remaining deviations should be smoothed out by averag- 
ing the observations over the 5 periods. 

3.4.4. Errors in the photometric flat-field correction 

The measurement noise contributes about 1 mmag to the flat- 
field correction uncertainty and is therefore negligible. The main 
weakness of the grid method is the inability to adjust variations 
of the photometric response during the ~ 17 minute observation 
sequences. 

The statistical properties of the variations of the overall zero 
point during short sequences can be obtained from the main sur- 
vey observation sequences that consist of half-hour long visits 
split into 5-6 exposures on the same field (with a negligible 
dithering). The zero point is found to be fairly stable during clear 
nights with an rms typically smaller than 5 mmag (see Sect. 4. 1 
for the determination of zero points and Fig. 10 for the distribu- 
tion of zero points variations). Most of this variation is uncorre- 
lated from one sequence to another. However, a small systematic 
shift of about +2 + 0.5 mmag per hour can be measured. 

We used the measured shift and rms as inputs to simulations 
in order to evaluate the impact of such variations on the determi- 
nation of the photometric correction. We checked for the influ- 
ence of a random variation of the photometric zero point with an 
rms of 5 mmag and a systematic shift of 2 mmag per hour. The 
systematic shift results in a negligible (below 1 mmag peak-to- 
peak) gradient on the focal plane. The random photometric zero 
point produces a more significant error on individual solutions 
taking the form of large-scale variations with a peak-to-peak am- 
plitude of about 1%. However, this random error is smoothed 
out as expected when several independent determinations of the 
photometric correction are averaged together. 

Some epochs display stronger (and time correlated) varia- 
tions of the observing conditions than what is typically observed. 
We took the conservative approach of discarding observations 
where conditions were not stable. Epochs were not used to flat- 
field the survey if the seeing variations were greater than 50% 



14 typically contained within 1%, most amplifiers being much more 
stable. 



Table 2. Internal consistency of the photometric flat-field. 
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Notes. 

(n) Estimated as the average of the rms of the relative zero point varia- 
tion on the focal plane (see appendix A). 

{b} Number of independent determinations of the photometric correc- 
tion. 

w This delivers a crude estimate of the rms of residual non-uniformity 
in the SNLS tertiary catalogs, which are obtained by averaging over the 
whole survey. 

peak to peak, or variations of measured aperture corrections 
were greater than 1.1% peak to peak. Those epochs are likely 
to display zero points variations greater than the 5 mmag rms 
assumed in the simulations. This cut removes 18, 23, 30, 28, 40 
and 20% of the available epochs in um^mSmJm, i%m an d Zm 
band, respectively. Valid observations are displayed in bold face 
in Table A.l. While this choice is somewhat arbitrary, exclud- 
ing or keeping the observations made under variable conditions 
has a small impact on the uniformity of the tertiary star catalogs. 
When included, the catalog are modified by about 1% peak-to- 
peak (or 3 mmag rms). 

3.4.5. Internal estimate of the flat-field uncertainty 

The typical r.m.s scatter of the selected grid solutions across 
epochs is 6, 6, 7, 7, 2 and 9 mmag in band u M ,g M ,r M ,i M and Zu< 
respectively. The structure of the variation is visible in Fig. 9. In 
all bands (but Hm that was installed after the cleaning), the accu- 
mulation of dust in the optics during the year 2006, visible at the 
lower right corner of the focal plane, dominates the variation. 
The pattern of noisy amplifiers also shows up on the variation 
map. Beyond these recognizable features it is difficult to distin- 
guish between what is attributable to actual modifications of the 
telescope response, and what should be considered as part of the 
flat-fielding error. 

As discussed above, the most reliable internal estimate of the 
remaining flat-fielding error is determined by the agreement be- 
tween different periods of the deep survey. The typical rms cry/ 
of the spatial discrepancies between individual exposures are ob- 
tained in Appendix A and summarized in Table 2. Assuming that 
the errors are smoothed out when different photometric correc- 
tions covering the 5 periods are averaged, one concludes that 
a rough estimate of the residual uncertainty on the tertiary star 
flat-fielding is given by JN^, where N g is the number of in- 
dependent determinations of the photometric correction. An in- 
dependent cross-check of the flat-fielding quality is obtained by 
comparison with the SDSS in Sect. 8.4; it gives similar numbers. 

3.5. Summary 

We have built a model of the effective MegaPrime passband cor- 
responding to the new Hm replacement filters. The im and vm 
models have also been significantly revised with respect to the 
previous R09 release of the SNLS calibration using the more re- 
liable 2006 filter measurement data. 

The flat-field processing has been slightly updated with re- 
spect to R09. The most significant change is the use of modeled 
color terms rather than measured color terms to translate mea- 
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Fig. 9. Variation of the photometric flat-field corrections. Each panel displays a low resolution map (36x36 super-pixels of 256x256 
real pixels, i.e. a resolution of about 2') of the rms between all the determinations of the photometric grid corrections in a given 
photometric band. CCD boundaries are shown by white dotted lines. Note that all photometric correction maps are normalized to be 
zero at the center of the focal plane. As a consequence the rms across epochs is null at the center by construction. Noisy amplifiers 
are easily seen in ru and zm bands. In all bands (but Hm that was not available), the variation of response induced by the 2006 
deposit of metal dust is clearly visible at the bottom right of the focal plane. 



surements when fitting the photometric corrections. This update 
is not expected to introduce significant changes in natural mag- 
nitudes except in vm at the corners of the focal plane, where 
the effect alters the flat-fielding by 0.5%. We also choose to 
make the instrument response constant for objects having an AB 
spectrum rather than for the primary spectrophotometric stan- 
dard BD +17 4708. This convention makes the derivation of AB 
magnitudes from instrumental magnitudes more straightforward. 
Another change is a smooth parameterization of the instrument 
response which is expected to provide a slight decrease in the 
flat-fielding noise. 

The main benefits for flat-fielding accuracy are related to the 
extension of the data sample (except for the broken im filter). 
Most sources of flat-fielding error are uncorrelated between dif- 
ferent epochs, and these errors are thus smoothed out with the 
addition of new grid observations. To avoid biasing the average 
we introduced a quality cut to remove epochs strongly affected 
by varying observation conditions. 

An upper bound on the flat-fielding error is provided by the 
variability of the photometric corrections determined on inde- 
pendent data samples. It is below 1% in all bands. Assuming 
the error is uncorrelated between epochs, this also sets an upper 
bound of typically 0.3 % r.m.s on the uniformity of the resulting 
catalogs. We study the consequences of flat-fielding uncertain- 
ties on calibration in Sect. 4.3.2. 



4. SNLS tertiary star catalogs 

4.1. Multi-epoch average photometry of tertiary stars 

A large number of exposures have been accumulated on the 
CFHTLS deep fields. A small dithering is generally applied to 
the observations for the purpose of filling gaps between CCDs. 
Its typical amplitude is between 80 and 120 pixels in right as- 
cension and 300 to 360 pixels in declination. This dither is small 
enough that position dependent changes in filter response can be 
neglected. We thus assume that a given star is always observed 
through the same effective filter. This assumption simplifies the 
construction of a multi-epoch photometric catalog using natural 
magnitudes for the selected tertiary stars in the deep fields. 

We use the stars selected as described in Sect. 11.1 of R09. 
We first determine aperture flux measurements for all stars in 
all the survey images, and then apply photometric flat-field cor- 
rections (Sect. 3.3) to the individual measurements according to 
Eq. (13). The survey is split into five periods, corresponding to 
noticeable changes in the instrument configuration. An indepen- 
dent determination of the photometric correction is applied to 
each epoch. The exact mapping between measurements and cor- 
rections is described in Appendix A and illustrated in Fig. A.2. 

The model for the measured flux Y Stt of non-variable stars 
after flat-fielding is: 

Y SJ = O s a, + n s4 (14) 
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where t indexes the exposure, s indexes the star, <J> S = KT 4 "' 1 * 
is the local broadband flux of star s, a, = 10° Az < is the photomet- 
ric zero point of exposure t, and n s t is the measurement noise. 
The overall scales of <$> s and a are not constrained from the data. 
Setting the scale of the O s to correspond to physical flux requires 
calibration data, and is the subject of section 5. At this stage, we 
are only interested in delivering an homogeneous set of measure- 
ments, and we will fit the model to the data with the constraint 
(a,) = 1 . We work in flux space rather than in magnitude space 
because at low S/N, averaging logarithms causes a bias. 

We opt for a two-step, iterative method. We start from a first 
estimate of $> s = <f>l \ obtained as the median of instrumental 
flux over all epochs for each star, and then we iterate the follow- 
ing steps until convergence: 

1. a^] = Y s<l /0, using the current estimate of <!>„ and solve for 
a, in the linear model: 



„(0 



a, + m s 



where the measurement noise on Y s t /<t>^ is modeled as an 
independent random Gaussian variable m s t centered around 
with variance cr 2 st . 

2. (D® = Y s t la, and improve the estimate of O s as the (non- 
weighted) mean of <I>® across a subset of photometric expo- 
sures 

In this procedure, +4cr J> , outliers are flagged and discarded. We 
model the relative flux uncertainty cr 2 s , as 



/, 2 + 



1 



7i 



(<D«)2 



(15) 



Accurate values of f}, and to a smaller extent, y,, are difficult 
to derive from first principles. We thus derived robust estimates 
of these parameters for each exposure by fitting the rms of af t 

computed in bins of O s , 

The major contributor to the y t term, which dominates the 
measurement error at low fluxes, comes from the Poisson fluc- 
tuation of the sky background. Its theoretically expected value is 
ft = -^pixCsicy^ where N p i x is the number of pixels in the aperture 
(< 700), cr 2 ^ is the variance of the background fluctuation, and 
G ~ 1.6e~/ADU is the gain of the readout electronics. This the- 
oretical estimate, however, is typically significantly smaller than 
the fitted y, value because it does not account for contamination 
by structured residuals in the background (PSF tails, ghost re- 
flections...), especially in the iu and zm bands where subtraction 
residuals from fringing can be important. 

The other important contribution to the measurement error 
is related to the flat-field error, The major components of 
this term are the error in the flat-field estimate, the variations 
of amplifier gains, and potential structures in the atmospheric 
extinction (clouds). Such errors are typically correlated among 
observations within the same night or the same run. This error 
term dominates the measurement error at the bright end: the er- 
ror on individual measurements of star magnitudes asymptoti- 
cally reaches 1.08/, la,, which typically amounts to 6 mmag for 
clear nights. 

The large number of epochs over a long time period makes it 
possible to control systematic errors by applying stringent selec- 
tion criteria (cuts) on the measurements entering in the averaging 
process. We applied two cuts based on the relative variations of 
the zero point: 1) the time stability of the zero point during each 
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Fig. 10. Top panel: Distribution of zero point variations dur- 
ing an observation sequence (typically 6 exposures of 220s in 
the same band in a row). Nights displaying a variability larger 
than 1% are discarded from the averaging. While the plot is 
cut at 0.025, numerous sequences (about 20%) taken under non- 
photometric conditions display a much larger variation. The zero 
point is defined as Z, = -2.5 log I0 a,. 



observation sequence, and 2) the spatial uniformity of the zero 
point in each exposure. 

The distribution of the standard deviation of a, from its mean 
value for the night is displayed in Fig. 10. The stability of the 
zero point is generally better than 1 % for photometric sequences, 
with typical values of 2-3 mmag. We discarded observations for 
nights with a relative variation of a, greater that 1% because they 
might have been affected by clouds; this requirement removes 15 
to 30 % of the nights. 

We also compute the relative variation of the zero 
point a* . on a per CCD basis, and compute AZj !) = 

max CC£ / 2.5 log 1() a { t ° ccd -mm ccd 2.5 log 10 or® cd . The distribution of 
AZ, is shown in Fig. 11: it is generally below 4%, with a typical 
value of 2%. This result shows that the various flat-field cor- 
rections are in agreement at this level, although there are oc- 
casional outliers related to atmospheric conditions or amplifier 
gain variability. These outliers are identified and discarded by a 
cut on AZ,. We discuss further the variation of AZ, with time in 
Appendix A. 

We did not attempt to correct the catalog for selection bias. 
Biases arise for faint objects due to the higher selection effi- 
ciency of positive fluctuations, and for bright objects that satu- 
rate some pixels on nights with good seeing, leaving only nights 
with poor seeing conditions. Completeness of the catalog, de- 
fined as the fraction of nights a given star is used, is shown in 
Fig. 12. Completeness does not reach 100% on average due to 
the dithering needed to fill gaps between CCDs, contamination 
of apertures by cosmic rays, and occasional defects in the instru- 
ment (part of the CCD mosaic being non operational at a given 
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Fig. 11. Distribution of the peak-to-peak zero point variation 
across the focal plane. Exposures displaying a peak-to-peak vari- 
ation greater than 3% (6% in um and zm) are discarded from the 
averaging. Vertical dashed lines show the cuts applied in the av- 
eraging process. 
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Fig. 12. Measurement completeness in the tertiary star catalogs. 
Each plot displays the ratio of selected measurements to the 
number of exposures as a function of the star magnitude in each 
band. The plateau at ~85% reflects mostly the ratio of the usable 
area in the focal plane to the total area covered by the catalog. 



time). We did not apply cuts to the catalog to discard stars po- 
tentially affected by selection bias. 

4.2. Anchoring to the Landolt system 

Nine fields containing a large subset of the Landolt equatorial 
secondary standards were routinely observed by the CFHT for 
calibration purposes. Those observations can be used to derive 
zero points for a large number of nights and tie the tertiary cata- 
logs to the Landolt system. 

The basic assumption underlying our strategy is that, be- 
sides the most obvious systematic differences between observa- 
tions (such as airmass differences or systematic differences in the 
PSF), 15 variations of observing conditions are uncorrected from 
one night to another, and therefore these variations will average 
out with many repeated measurement of the flux ratios. We are 
thus modeling the exposure zero point as: 

Z, = Z„ + 2.5 log 10 (* e *p) + (X, - l)k atm + 6Z,„ (16) 

where Z, is the zero point of a given exposure t, Z„ is the mean 
zero point for night n, (X, - l)k. dtm is the average extinction at 
airmass X, and SZ,„ accounts for a variation of the photometric 
response in exposure t with respect to the average of the night Z„. 
In the calibration process, we treat this latter term as a noise that 
is independent for each night, and whose average is zero over 
many nights. Deviations from this hypothesis, related to system- 
atic differences between calibration and science exposures, are 
discussed in Sect 4.3. 

Following the scheme described in R09, we assume that 
the transformation between the Landolt system and the refer- 



to be discussed in section 4.3 



ence MegaCam system follows an approximate piece-wise lin- 
ear color relation f\j. 

m lX0 =L + f L (C) + SL-6, b (17) 

where L is the Landolt magnitude of the star in the closest pho- 
tometric band, C is the Landolt color index of the star, and 5L 
accounts for small deviations of individual stars with respect to 
the single parameter relation /l. The calibration offset <5 a b sets 
the magnitude system to the AB flux scale, and is determined by 
referring to a primary standard which is the subject of Sect. 5. 
For now we ignore 6 d b and set the photometric calibration of the 
tertiary catalogs to a common but arbitrary scale, close to a Vega 
system. 

The large number of observations allows us to solve simulta- 
neously for the transformation and the exposure zero points. We 
impose the constraint that 5L — on average, and find no signif- 
icant degeneracy between parameters. The complete model for a 
photometrically flat-fielded Landolt star measurement is: 

m A DU - Sk(x)(c\ Xo - c°(x)) = 

L + f L (C) + 6L-Z„- 2.5 log 10 (f exp ) - (X, - l)fe atm 

The relation /l is parameterized as: 




aC if C < Cbreak 

acbreak + P(C - c break ) otherwise. 



The free parameters include the two slopes for the piece-wise 
linear color transform (a and /}), the airmass correction term k. dtm , 
the zero point for each night Z„, and the deviation 5L for each 
star. We impose the constraint < 5L >« by adding a prior term 
Yts dL 2 /o- 2 s to the^- 2 function. We choose cr, = 0.01 to match the 
observed dispersion of the MegaCam to Landolt relation. The 
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Fig. 13. Color transformation between the Landolt system and the reference MegaCam system. Open black circles show the 
MegaCam and Landolt measurements of the Landolt secondary standards. The overall scale is set by the magnitudes of the pri- 
mary spectrophotometric standard star BD +17 4708 (solid red square). Measurements of BD +17 4708 magnitudes in the Landolt 
system are provided by Landolt & Uomoto (2007). Since BD +17 4708 cannot be directly observed with MegaCam, its MegaCam 
magnitudes are inferred from the average color transformation (black line) determined from the secondary standard measurements. 
Potential deviations from this average color transformation are examined in R09 (see Sect. 5.1). Other primary standards from 
Landolt & Uomoto (2007) are displayed as red crosses; they are not close enough to the bulk of secondary standards to be reliably 
transformed with the average transformation. Observations of the Landolt standards in Hm were not numerous enough to deliver a 
robust measurement of the i2\ Xo - I transformation (refer to Sect. 5.3 for the Hm calibration). 

Table 3. Landolt to MegaCam average transformation parameters defined in Eq. (18). 



band 


color 


color range 


f-break 




a 




P 




t-atm 




k, b 


| vatm 


SL" 


Um 


U - 


V 


0-1 


0.40 


-0.2553 


+ 0.0082 


-0.3120 


± 0.0094 


0.3075 ± 


0.0036 


0.3052 ± 


0.0028 


0.3161 




gM 


B - 


V 


0-2 


0.45 


0.5320 


±0.0141 


0.4309 


± 0.0023 


0.1235 ± 


0.0006 


0.1296 ± 


0.0011 


0.1422 


0.007 


I'M 


R- 


1 


0-0.8 


0.65 


0.1583 


± 0.0058 


0.1493 


± 0.0295 


0.0629 ± 


0.0007 


0.0762 ± 


0.0010 


0.0785 


0.003 


hi 


R- 


I 


0-0.8 


0.40 


0.2200 


±0.0122 


0.1446 


± 0.0070 


0.0309 ± 


0.0007 


0.0402 ± 


0.0010 


0.0321 


0.002 


Zm 


R- 


I 


0-0.8 


0.35 


-0.1713 


± 0.0267 


-0.4483 


±0.0105 


0.0154 ± 


0.0014 


0.0367 ± 


0.0015 


0.0206 


-0.010 



Notes. 

(<I) Average airmass term measured on data. 

( *' Average airmass term measured on data for photometry in fixed aperture radii of 43 pixels. The difference between the two is attributed to 
correlation between airmass and aperture corrections. 

(c) Airmass term according to the average extinction model at Mauna Kea from Buton et al, (2012). 
(</) Offset affecting the BD +17 4708 color transform derived in R09 table 7. 



adjusted transformations are shown in Fig. 13. Corresponding 
values for the transform parameters are given in table 3. 

The fit of Landolt observations delivers measurements of 
zero points Z„ for a significant subset of the nights with obser- 
vations of the SNLS fields (see table 4). These zero points can 
be used to anchor the flux scale of the multi-epoch catalogs of 
the 4 fields to a common scale (relative to the Landolt system) 
according to the model given by Eq. (16). 16 Non-photometric 
nights that were discarded in the analysis of the multi-epoch cat- 
alog (see Sect. 4.1) were also rejected in this averaging. Since a 
significant time period can occasionally separate the calibration 

16 Anchoring to the Landolt system also enables the use of spec- 
trophotometric standards measured in the Landolt system to set the ab- 
solute scale of the catalogs (see Sect. 5.1). 



exposures and the science exposures, this cut does not exclude 
all calibration observations affected by non-photometric condi- 
tions. A more robust assessment of the photometricity is possi- 
ble when several calibration exposures are available for a single 
night, as is the case for about a third of the nights. We also dis- 
carded nights displaying a zero point that is significantly lower 
(by more than 0.05 mag) than the average zero point observed in 
a 10-night rolling average. Finally, we discarded nights display- 
ing 3 cr zero point outliers when applying the calibration model 
from Eq. (16) to calibrate the deep field exposures. Table 4 sum- 
marizes the number of photometric nights with concomitant ob- 
servations of the Landolt and science fields, the observed dis- 
persion between the independent measurements of the catalog 
scale provided by each night, and the resulting statistical uncer- 
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Table 4. Statistical uncertainty (in magnitudes) on the calibra- 
tion to the Landolt system. 







Dl 






D2 




band 


# nights 


u 


cr/ VTV 


# nights 


cr 


cr/ 


u M 


7 


0.007 


0.003 


1 


0.000 


0.000 


gM 


53 


0.008 


0.001 


44 


0.011 


0.002 


r M 


63 


0.009 


0.001 


64 


0.011 


0.001 


hi 


68 


0.010 


0.001 


57 


0.008 


0.001 


Zm 


30 


0.014 


0.002 


25 


0.012 


0.002 




D3 


D4 




8 


0.010 


0.004 


8 


0.015 


0.005 


8m 


56 


0.010 


0.001 


55 


0.009 


0.001 


I'M 


74 


0.011 


0.001 


64 


0.014 


0.002 




63 


0.012 


0.001 


51 


0.010 


0.001 


Zm 


25 


0.013 


0.003 


27 


0.013 


0.002 



tainty. The smaller number of zm nights is partly related to fewer 
Landolt observations in zm and cuts. Observations in u M were 
not part of the supernova survey and were not cadenced. 

4.3. Systematic uncertainties in the calibration to the Landolt 
system 

We now turn to a review of potential sources of systematic uncer- 
tainty in the calibration path through Landolt observations. Here 
we only discuss uncertainties that affect the determination of the 
absolute flux scale of the tertiary catalog. These uncertainties are 
summarized in Table 5. 



4.3.1. Systematic differences in aperture corrections 

The first is from potential systematic differences in the photom- 
etry of standard and tertiary stars. The brightness of standard 
stars requires a short integration time (2-3 s) to avoid satura- 
tion. During the course of the survey, a small defocus (smear- 
ing the PSF core by about 1 pixel) on Landolt fields has also 
been applied to observe the brightest secondary standards of the 
Smith et al. (2002) set. It is likely that this defocus introduces 
systematic differences between the effective PSF of tertiary star 
observations and primary or secondary standard observations, 
resulting in slightly different aperture corrections. 

To first order, PSF variations are approximately accounted 
for by the scaling of aperture radius in our photometry method. 
We statistically investigated deviations to this correction by com- 
puting aperture corrections between aperture of radius 7.5 and 20 
times the image quality. We found a remaining variation of aper- 
ture corrections with image quality, and offsets between Landolt 
and deep field exposures dependent on the defocusing of Landolt 
observations. 

The change in aperture corrections with image quality is 
illustrated in Fig. 14. The changes for tertiary standards and 
Landolt stars are similar, except for the offset related to defo- 
cusing. The IQ dependence is approximately linear with a small 
positive slope in the gM band and a negative slope in the other 
bands. 

The observation strategy of the Landolt fields slightly 
evolved over the course of the survey. Calibration stars were 
originally exposed for 3 s with the telescope in focus. The ex- 
posure time was reduced to 2 s in February 2005 and the tele- 
scope focus was offset by 0.5 mm with respect to the nominal 
focus to avoid saturation of the brightest stars. The amount of 
defocus was later reduced to 0.3 mm in August 2007. We there- 




L.5 2.0 

IQ(pixels) 



2.0 2.5 

IQ(pixels) 



Fig. 14. Variation of aperture correction from small to large aper- 
tures as a function of the image quality, IQ, defined as sigma of 
a Gaussian fit to the PSF (in pixels). Apertures are scaled with 
the IQ: small apertures correspond to a radius of 7.5xIQ and 
large apertures to a radius of 20xIQ. The aperture corrections 
are computed as if 75 /<f> 20 - 1 where (p 15 and <p 20 denote instru- 
mental fluxes in the small and large aperture, respectively. Data 
points correspond to the median of aperture corrections in bins 
of IQ over all the exposures in the survey. Results for short (and 
slightly out-of-focus) exposures on Landolt fields (black circles) 
and long exposures on deep fields (black squares) are displayed 
separately. 
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Fig. 15. Average aperture correction (between 7.5 and 20 IQ ra- 
dius) is shown for the 3 periods corresponding to different ob- 
servational setups. Landolt stars (short exposures, slightly out of 
focus) are shown in black circles, and deel fields (long expo- 
sures) are shown in black squares. The statistical uncertainties 
on the averages are negligible. 



fore have three distinct time periods corresponding to different 
observational setups. We computed the average aperture correc- 
tion offset between Landolt and deep field observations in bins 
corresponding to each time period. The measured offsets, shown 
in Fig. 15, were applied to the calibrated magnitudes to correct 
for the systematic differences on aperture corrections in the de- 
termination of the zero points. 

A reliable estimate of the error on this correction is difficult 
to obtain. The correction is based on the assumption that the to- 



18 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 




0.010 
0.005 
0.000 
-0.005 
-0.010 



' I I I I I - 


1 

: ♦ 


P + ! 


♦ 


t 


:, , 1 , , , , 1 , , , , 


• 



0.010 
0.005 
0.000 
-0.005 
-0.010 



■II 


I 




t 1 






U t t 




:, , 1 , , , , 1 , , , , 1 , , , , 1 , . , , 


,'} 



Fig. 16. Catalog zero point as a function of the grid solution. The 
error bars reflect the statistical uncertainty on the mean estimated 
from the dispersion of measurements in the period. 



tal flux loss is indeed proportional to the aperture correction, or 
in other words, that the fraction of flux captured in the larger 
aperture is stable. Whether this assumption is verified or not is a 
tricky question. A model independent hint on the accuracy of this 
correction is indicated by the agreement between the three time 
periods once the correction is applied, and by the residual dis- 
persion of the calibration as a function of seeing difference. The 
different periods are statistically compatible. To be conservative, 
we took the difference between the two periods with the greater 
statistical accuracy (2003-2005 and 2005-2007) as a reasonable 
upper bound on the systematic related to aperture correction. 

4.3.2. Flat-field error 

The first effect of the flat-field error is to degrade the unifor- 
mity of the tertiary catalogs as discussed in Sect. 4.1 and in 
Appendix A. It can also result in an error on the average cali- 
bration because secondary standards are observed preferentially 
at the center of the focal plane. The Landolt calibration is made 
with sufficient statistics to provide a strict limit on this effect. 
The dispersion of the calibration residuals as a function of the 
flat-field solution applied is shown in Fig. 16. Each point on 
the plot corresponds to a mostly independent realization of the 
SNLS calibration, as it involves disjoint data samples for both 
the determination of the flat-field and the zero points. The mea- 
sured dispersion is between 1 and 3 mmag in all bands. This dis- 
persion is dominated by the statistical uncertainty on each point, 
and can be used to set an upper bound on the contribution from 
the flat-field error. Since the flat-field error is likely to be un- 
correlated between different observations of the grid fields (cf. 
Sect. 3.4), we divided the measured dispersion by the square- 
root of the number of grid observations to estimate the flat-field 
error in Table 5. 



4.3.3. Shutter precision 

The precision of the MegaPrime shutter could bias the photomet- 
ric ratio between short and long exposures. The uncertainty on 
the exposure time was estimated in R09 to be less than 3 ms. For 
photometric standards exposed for 2 s, this uncertainty affects 
the flux by at most 0.15%. 



4.3.4. Residual background contamination 

The sky background subtraction is a difficult problem for aper- 
ture photometry. Our background estimate is based on the dis- 
tribution of image pixels in tiles of 256 x 256 pixels after ap- 
plication of a mask discarding objects up to a given contamina- 
tion level. The final estimate is a smooth interpolated version of 
the resulting low-resolution background map. Higher frequency 
structure remains in the background subtracted maps, in partic- 
ular, far in the PSF tails from bright objects. Since R09 showed 
significant hints for background residuals in the aperture pho- 
tometry of tertiary stars, we investigated the origin of this bias. 

The basic idea is the following: structure in the background 
is expected at scales smaller than 256 pixels. Important contribu- 
tors to such structure are PSF tails from bright objects and resid- 
uals from fringe subtraction in the i and z bands. These struc- 
tures remain in the background subtracted images and contam- 
inate aperture photometry. This is an intrinsic limitation of the 
aperture photometry method which assumes isolated objects on 
a background with noise expected to be zero on average, an hy- 
pothesis not verified in practice. 

PSF tails are particularly annoying in our case as they make 
the average level of aperture contamination dependent on the star 
location. For our tertiary stars, which are selected to be fairly iso- 
lated objects, one expects the average contamination of apertures 
to be lower than for randomly distributed objects. Therefore, our 
hypothesis to explain the bias spotted in R09 is that the back- 
ground subtraction algorithm is ill-suited to the specific loca- 
tions of our tertiary stars (typically an overestimate as they are 
isolated). 

We seek a quantitative estimate of this effect. Let us denote 
b(x) the residual background level at a given position x in the 
image after background subtraction, and C(x) the flux carried to 
this position by the PSF tails of neighboring objects. According 
to the discussion above, b and C are related. This relation is 
rather complex as the specifics of the background subtraction al- 
gorithm comes into play to suppress the large scales of the struc- 
ture. In addition, the value of C is difficult to obtain because it 
requires knowledge of the PSF far out in the tails. Estimating the 
effect from first principles is thus not feasible. 

We can characterize the residual background statistically 
by measuring an empirical relation between b and a quantity 
crudely related to C. We define a contamination index C(x) as: 



C(x) = ® t}f(x - x„) 



(19) 



where i// is an approximation of the PSF normalized to 1, o runs 
across objects detected in the image and is the flux of object 
o. The actual value of C does not have a clear physical inter- 
pretation because it depends on the number of objects included 
in the sum and on the approximation used to describe the PSF. 
We can expect, however, that positions corresponding to larger 
C correspond to positions with greater residuals. 

To determine the empirical relation between b and C, we 
sampled the deep field images in the survey at uniformly spaced 
positions. For each position x we compute two quantities. First 
we estimate b(x) from the average residual sky level (after back- 
ground subtraction) in a 32x32 pixel aperture centered at this po- 
sition. Second, the contamination index C(x) is computed under 
the assumption that the PSF follows a Moffat function (Moffat 
1969) with parameters adjusted on the images. The top panel of 
Fig. 17 shows clear evidence for the dependence between the 
residual background level and the isolation of the location. 
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The average contamination in tertiary standard measure- 
ments is obtained as follows: we first compute the contamination 
index (Eq. 19) for the positions of our tertiary standards, exclud- 
ing the contribution of the tertiary star itself. We then convert 
each contamination index into a residual background level ac- 
cording to the relation shown in Fig. 17, and compute the aver- 
age over all tertiary star positions. The distribution of the con- 
tamination index of tertiary stars is shown on the bottom panel 
of Fig. 17. As expected, the tertiary standards are well isolated 
objects with other image objects adding less than a fraction of an 
ADU per pixel. We found an average of 0.06, 0, -0.26 and 0.01 
ADU/pixel, respectively, in the gM im and zm bands. 

Those contamination numbers are in perfect agreement with 
the estimate of background residuals obtained in R09 Sect. (4.2) 
by studying aperture corrections as a function of flux. We thus 
conclude that the proposed mechanism explains the observed 
background subtraction bias. To account for the resulting bias 
on aperture magnitudes, we subtracted the average contamina- 
tion level to the aperture flux measurements. Note that this is 
only an average correction, accounting for the overall bias intro- 
duced by background residuals. Contamination of apertures by 
background structure still contributes noise on individual aper- 
ture measurements. The distribution of this noise has a positive 
skewness (whose significance depends on the magnitude range 
considered) and some caution should be exercised when using 
the aperture catalog for calibration. 

To illustrate the amplitude of the remaining effect, we give 
a comparison of catalog magnitudes with PSF photometry. The 
PSF photometry software independently adjusts and subtracts a 
background offset for each measurement, and is thus insensi- 
tive to residual background fluctuations. A known chromatic ef- 
fect (see Guy et al. 2010), attributed to wavelength-dependent 
changes in the PSF that were neglected in the PSF photometry 
method, 17 has to be accounted for before performing the compar- 
ison. The dispersion of the difference between aperture and PSF 
photometry of the tertiary stars as a function of the star color are 
displayed in Fig. 18. The chromatic difference is easily seen in 
the gM band. Smaller but significant trends are also visible in the 
rM and Im bands. 

After removing the chromatic trend and correcting the aper- 
ture magnitudes for the average background contamination com- 
puted above, the PSF and aperture magnitudes agree as a func- 
tion of the star magnitude. As illustrated in Fig. 19, there is ex- 
cellent agreement in the magnitude range 16-21 for the gM and 
r M bands, indicating that the average aperture contamination is 
mostly removed. The range of agreement appears a bit smaller 
in the im and zm bands: 16-19.5. To obtain the most precise cal- 
ibration, we suggest using stars only in the above confidence 
ranges. 

Applying the same analysis to the Landolt measurements, 
we found that the residual background level in short time expo- 
sures is consistently positive and varies around a mean value of 
0.1 ADU/pixels. This is attributed to a small biasing of our back- 
ground estimator, which is based on a combination of the mean 
and the median of pixel values. It delivers biased estimates on 
the asymmetric distribution of background pixels in short expo- 
sures. The impact of such a bias on the determination of zero 
points is smaller than 1 mmag and is included in the systematic 
error budget. 



17 This is handled in SNLS by interpreting the PSF magnitudes in a 
modified photometric system (see Guy et al. 2010 for further discus- 
sion). 
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Fig. 17. Background residuals as a function of the contamina- 
tion index. The contamination index at a given location for a 
given aperture is defined as the flux in the aperture expected from 
objects detected in the field, assuming they have a Moffat PSF 
profile (see Eq. 19). The top figure displays the average back- 
ground level (in ADU per pixel) in random apertures as a func- 
tion of this contamination index. We concentrate on random po- 
sitions with low contamination index (apertures expected to be 
mostly empty), sampling the background level far from detected 
objects. The shape of the curves reflect the remaining structure 
in the background due to PSF tails. The bottom plot displays 
the distribution of the contamination indexes at the position of 
tertiary stars (excluding the contribution from the tertiary star it- 
self). Horizontal lines in the top plot display the average value 
of the curve weighted for tertiaries, i.e. the expected per-pixel 
background offset for tertiaries. 
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Fig. 18. Difference between PSF and aperture photometry for 
tertiary stars as a function of star color. The black line displays 
a linear fit of the chromatic trend observed. 



20 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 



CD 
Q_ 

cp 

LL 

co 

0- 



0.02 
0.01 
0.00 
-0.01 



-O-0.02 



— 0.02 
j CD 
~3 9r 0.01 

CO 

£ ooofc- 

5^-0.01 Et 
-0.02 



1 I " " 1 1 



0.02 - 



16 17 18 19 20 21 22 

Km 

\<<><]<<<<\><><\<">\<> 



16 17 18 



19 20 21 22 



CD 






Q_ 


0.01 




cp 






LL 


0.00 


*■ 


CO 






CL 


-0.01 




.3 








-0.02 





0.02 r 



CD 

9r o.oi 

CO 



LL 
CO 
0- 



0.00 



-0.02 



- 1 1 1 


1; I . 
















:....l....l....r...:wft/;,;n,..: 



17 18 



19 20 21 22 



16 17 18 19 20 21 22 

ZM 



Fig. 19. Difference between PSF and aperture photometry for 
tertiary stars as a function of star magnitudes (with the chromatic 
trend from Fig. 18 and the average background contamination of 
apertures subtracted). Black circles display the average of data 
points in bins of magnitudes. Averages are computed with a clip- 
ping at 2.5cr as used in the calibration of SNLS supernovae light- 
curves (Guy et al. 2010). 



4.3.5. Survey uniformity 

The uniformity of the photometry between the 4 SN fields relies 
on the uniformity of the photometry of Landolt equatorial stan- 
dards. The dispersion of zero points between individual Landolt 
fields has been estimated in R09 to be 0.002 in gri and 0.004 
in z. As the true error is likely to have some spatial correlation, 
we use these numbers as estimates of the systematic uncertainty 
induced by non uniformity in the secondary network. Note that 
this error potentially affects the average calibration as the pho- 
tometric standard observations are calibrated against a subset of 
the secondary standards. 

The uniformity is confirmed by the generally good agree- 
ment of stellar locus between fields as displayed in Fig. 20. The 
dispersion of the stellar locus is compatible with the statistical 
uncertainty on the zero point in each field. 



indirectly provided by Landolt & Uomoto (2007) observations 
of HST standards in the Landolt system. 

However, recent work (Bohlin & Cohen 2008) has extended 
the CALSPEC database to include redder and fainter standards. 
Direct MegaCam observation of these new standards opens an 
alternate calibration path for the SNLS fields. 

We start this section by deriving AB offsets for the tertiary 
catalogs from the Landolt calibration observations. We then ex- 
plore the possibility of direct calibration through observations of 
three solar analogs. Cross-calibration with the SDSS, presented 
in Sect. 8, offers yet a third calibration alternative. We devote 
section 9 to the comparison of these various options. 

5.1. Absolute calibration through Landolt observations 

5.1.1. Transformation of Landolt standard measurements 

We have used MegaCam observations of Landolt secondary 
standards to adjust the color relation shown in Eq. (17). This re- 
lation can be used to transform Landolt system measurements of 
stars (particularly measurements of standard stars) to the tertiary 
standard-defined SNLS system. Calibration constraints can be 
obtained by applying this relation to a spectrophotometric stan- 
dard whose expected MegaCam AB magnitude fh\ Xo can be syn- 
thesized from the calibrated CALSPEC spectrum. This gives: 



6 ah = L + f L (C) + 6L- fh\ XB . 



(20) 



Any of the Landolt & Uomoto (2007) Primary 
Spectrophotometric Standard Stars (hereafter PSSS) can 
be used to set the unknown calibration coefficient <5 a b, provided 
the corresponding 5L are known. The red crosses on Fig. 13 
show primary standards from the CALSPEC database with 
accurate Landolt system magnitudes. 

The average color transformation / L is adjusted so that 
(6L) — on average for Landolt secondary stars. As shown 
in Fig. 13, the only spectrophotometric standard with Landolt 
magnitudes sufficiently close in color to the average of Landolt 
secondary standards (so that the corresponding 6L is expected to 
be small) is BD +17 4708. R09 has done a careful estimate of 
the BD +17 4708 5L values, reproduced here in table 3. We use 
these values to compute estimates of the calibration offsets ac- 
cording to Eq. (20). These calibration offsets are given in table 6. 



5. MegaCam absolute calibration measurements. 

The straightforward way to set the physical flux scale of the star 
catalog is to equate synthetic magnitudes of spectrophotometric 
standard stars with their corresponding magnitudes in the sys- 
tem defined by the catalog. Unfortunately, no such standards are 
currently present in the SNLS fields, although several pure hy- 
drogen white-dwarves have been identified (Lidman et al. 2012; 
Limboz et al. 2008) that might serve this purpose in the future. 
For now, our absolute flux calibration relies on external observa- 
tions of spectrophotometric standard stars. 

Measurements from the Hubble Space Telescope imaging 
spectrograph (STIS) provide the best available set of spectropho- 
tometric standard stars in the visible part of the spectrum. The 
calibrated spectra of HST fundamental standards are part of 
the CALSPEC database. Most of these stars are too bright for 
MegaCam to routinely and reliably observe. Therefore, until 
now, observations of the Landolt equatorial secondary standards 
(Landolt 1992) have been the preferred method to set the SNLS 
field absolute calibration, with anchoring to the HST flux scale 



5.1.2. Uncertainties on the MegaCam magnitudes of 
BD +17 4708 

The dominant contribution to the systematic error budget is the 
uncertainty on 6L, the term quantifying a star's departure from 
the average transformation between the Landolt system and the 
MegaCam system. Intrinsic dispersion around the average trans- 
formation is expected due to variations in the SED of individ- 
ual stars. Without any further assumptions, the dispersion of 
Landolt secondary standards around the average color transfor- 
mation gives an estimate of the related uncertainty. Dispersion 
ranges from 13 mmag rms in band gM to 20 mmag in band zm- 

As already stated, a more precise determination of 6L, i.e. 
of the exact position of BD +17 4708 with respect to the bulk 
of the Landolt secondaries in the Fig. 13 color-color diagram, 
was attempted in R09. This work used synthetic photometry on 
stellar libraries to model the small deviations induced by spe- 
cific properties of BD +17 4708 (extinction, metallicity, logg, 
binarity) with respect to the mean properties of the Landolt stan- 
dards, leading theoretically to a precise determination of 8L and 
hence, to accurate magnitudes for the standards in the MegaCam 
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g-l g-i g-i 

Fig. 20. Top: Color-color diagram of the 4 SNLS tertiary star catalogs. Bottom: Residuals to a linear fit of the average stellar locus 
in the range 0.2 < g - i < 1.2. The residuals for the 4 deep fields are displayed separately (Dl: red crosses, D2: green squares, D3: 
black circles, D4: blue asterisks). 



system. Unfortunately a few contributions, like the effect of a 
faint companion to BD +17 4708, were difficult to estimate accu- 
rately and quite large uncertainties remain on the BD +17 4708 
MegaCam magnitudes, especially in zm- The R09 6L values are 
given in the last column of table 3; R09 uncertainty estimates on 
6L are shown in the 'Landolt to Megacam transformation' row 
of table 5. 

In band um the value of 6L has very large uncertainty, large 
enough that its determination was not even attempted by R09. 
Consequently, we quote the dispersion of the secondary stan- 
dards around the color transformation as the uncertainty on the 
MegaCam magnitude of BD +17 4708. 18 It amounts to 0.02 
magnitudes in m m . 

The Landolt measurements of BD +17 4708 magnitudes are 
also affected by a small statistical uncertainty (about 0.002) that 
we include in the systematic error budget. The uncertainties af- 
fecting the Landolt observations calibration method (including 
relevant uncertainties discussed in Sect. 4.3) are summarized in 
Table 5. 

5.2. Direct calibration to the HST 
5.2.1 . Observations of 3 solar analogs. 

The addition of fainter standards to the CALSPEC database 
raises the prospect of a more straightforward calibration method. 
We investigated the possibility of using direct observations of 
the faintest PSSS to set the photometric calibration of the SNLS. 
Three stars, SNAP-2, P177D and P330E, lie close enough to the 
D3 field to be easily observed together in a single short observ- 
ing block with a potentially small airmass lever arm. 

Observations were conducted at CFHT in June 201 1 in short 
sequences, with exposures of standard stars bracketing a single 

18 The dispersion is measured on secondary stars in a limited color- 
range enclosing the color of BD +17 4708. 



exposure on D3. To avoid saturation and to study the behavior of 
aperture corrections on out-of-focus exposures, two exposures 
(one in focus, one out-of-focus) were taken of the brighter stars 
P177D and P330E. Exposures on SNAP-2 begin and close the 
sequence, enabling the monitoring of atmospheric extinction sta- 
bility. Standard fields are exposed for 3 s; D3 fields are exposed 
for 120 s. 

Images were processed by the standard Elixir pipeline and 
aperture photometry was performed in several predefined radii. 
Images were flat-fielded using the appropriate grid solution. 
Standard stars were observed on CCD #22. D3 flux measure- 
ments are matched with the tertiary star catalog and used to de- 
termined a zero point for the sequence. The zero point is then 
used to calibrate the standard star measurements, correcting for 
the effect of different airmasses. This allows us to determine the 
offset between the Landolt based calibration and an AB system 
by: 

£ ab = -2.5 lOg lo (0ADuAexp) + Z + k itm (X D -X s )- fh\ x (21) 

where 0ADu/ f exp is the instrumental aperture flux corrected for 
PSF variations between the current exposure and the deep field 
exposure on which the sequence zero point Z has been deter- 
mined, Xo and X s are the respective airmasses of the deep field 
and the standard star, lc atm is the average linear correction for 
atmospheric extinction and fh\ x is the natural AB magnitude of 
the standard star at the observed position synthesized from its 
CALSPEC spectrum. Again, we use the notation fh\ x in place of 
m\x to distinguish synthetic from measured magnitudes. 

5.2.2. Error model 

We build an error model for these measurements by consider- 
ing the following sources of uncertainty: the photon noise of 
the standard star measurements, the flat-field noise, the tertiary 
measurement noise, the time variation of the zero point during 
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Table 5. Uncertainty budget in the cross-calibration between SNLS and the HST via the observation of Landolt secondary standards. 





u M 


Sm 








Aperture corrections 


±0.0045 


±0.0023 


±0.0014 


< 0.001 


±0.0036 


Flat-field error 


±0.0012 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Landolt uniformity 


±0.0040 


±0.0020 


±0.0020 


±0.0020 


±0.0040 


Shutter precision 


±0.0015 


±0.0015 


±0.0015 


±0.0015 


±0.0015 


Background residuals 


±0.0010 


±0.0010 


±0.0010 


±0.0010 


±0.0010 


BD +17 4708 landolt measurements 


±0.0034 


±0.0027 


±0.0028 


±0.0028 


±0.0030 


Landolt to MegaCam transformation 


±0.0228 


±0.0023 


±0.0043 


±0.0026 


±0.0178 


Total 


±0.0239 


±0.0050 


±0.0060 


±0.0047 


±0.0189 



Notes. (a) The i2 M band can be calibrated relatively to i M with a precision of 2 mmag (cf. Sect. 5.1). 



the sequence and the uncertainty on the airmass correction. We 
are only interested in uncertainties affecting the flux measure- 
ments, for now. Any error on the calibration offset related to 
uncertainties affecting the synthetic magnitude m\ x , such as un- 
certainties on passbands and error on the reference spectra, will 
be discussed later (in Sect. 9) and are not considered in what 
follows. 

The photon noise can be readily determined from the sky 
level and the magnitude of the star. The typical contribution of 
background residuals to aperture contamination has been mea- 
sured on the numerous observations of the Landolt fields. We 
use the measurements from the D3 exposure to estimate the 
size of the flat-field noise for the sequence. We adjust the error 
model (15) to the dispersion of measurements with respect to the 
magnitude of the corresponding tertiary catalog star. Typically, 
the corresponding flat-fielding error term (/) is found to be 
smaller than 0.4%. 

The 72 tertiary stars on CCD # 22 allow a zero point determi- 
nation with typical statistical uncertainties of 2 mmag in bands 
giw^MiM, and 5 mmag in bands zm and um- This error affects all 
measurements of the sequence in the same way but is expected 
to be uncorrelated from one night to another. 

The photometric ratio between the two SNAP-2 exposures 
bracketing each sequence was used to estimate gray atmospheric 
extinction variation for the corresponding set of observations. 
The extra noise contribution of this component is found to be 
smaller than 2 mmag. 

Due to external constraints in the observing schedule, the 
deep field observation airmasses tend to be higher, on average, 
than the primary standard observation airmasses. This is unfortu- 
nate as it makes the final measurement sensitive to the variations 
of the atmospheric extinction around its average value k. dtm {cf. 
table 3). 

Computations of the night-to-night dispersion of the SNLS 
zero point were used to estimate the typical variation of extinc- 
tion around its average value. Magnitude variations were found 
to be 0.028, 0.009, 0.01, 0.01, 0.015 in ugriz bands, respectively. 
These estimates do incorporate contributions unrelated to atmo- 
spheric extinction such as instrumental effects, variation of the 
PSF uncorrected by the IQ scaling rule, and noise from the zero 
point determination itself. Therefore, effective variation of at- 
mospheric extinction on clear nights may actually be somewhat 
smaller. On the other hand, these numbers are in good agree- 
ment with similar SNFactory measurements (Buton et al. 2012), 
except for zm where our estimate is significantly larger than the 
SNFactory estimate (about 0.3%), and um where they found a 
variability closer to 4%. We used the largest numbers in each 
band as a conservative estimate of the uncertainty on the airmass 
correction term. 



5.2.3. Systematic errors on the solar analog measurements 

The solar analog-based measurement of the SNLS tertiary cat- 
alog AB offsets is also affected by systematic differences in the 
PSF between short and long exposures, uncertainty on short ex- 
posure durations, residual contamination of apertures, and, as 
the observations are coupled only with the SNLS deep field D3, 
non-uniformity in the survey. 

We carefully monitored the variation of aperture corrections 
between exposures by computing aperture photometry in a large 
panel of radii. We choose to use 16-pixel radius aperture mea- 
surements to limit the background noise contribution to individ- 
ual measurements. We also measure the median photometric ra- 
tio between these apertures and a 43-pixel radius aperture, which 
is assumed to capture a stable fraction of the total flux. All mea- 
surements were corrected by this measured photometric ratio. 
The statistical uncertainty on this aperture correction for individ- 
ual exposures is typically smaller than 1 mmag and is included 
as extra noise. 

The important matter is how well the 43-pixel aperture is 
representative of the total flux. A hint is given by the variation 
of the photometric ratio between the first and last exposures of 
the sequence. Figure 21 displays this photometric ratio for the 
43-pixel aperture as a function of the difference in image quality 
between the two exposures. While the seeing can vary a lot, the 
photometric ratio is typically stable, with a rms of 2.3 mmag. 
This sets an upper bound on the combined contribution of the 
atmospheric extinction variation, shutter noise and variation of 
the flux loss outside the 43-pixel aperture. We conservatively use 
this number as a (pessimistic) estimate of the systematic flux 
difference between photometry on short and long exposures. 

The background residual is measured for each exposure on a 
per CCD basis and all measurements are corrected for it. The ex- 
tra noise introduced by the contamination of individual apertures 
is included in the measurement covariance matrices. 

The uniformity of the survey has been discussed in 
Sect. 4.3.5. We use the statistical uncertainty on the average zero 
points of the individual SNLS D3 field, as given in table 4, as an 
estimate of the departure of this field from the average survey 
system. 

Finally, as for the previous calibration method, the shutter 
precision potentially biases the flux estimate of the short stan- 
dard exposures. However, no field has been exposed less than 
3 s, decreasing the impact of this systematic to less than 1 mmag. 

5.2.4. Calibration offsets from the solar analogs. 

Our error estimates were used to build a covariance matrix R 
of measurement errors accounting for correlation between the 
different error sources. The final AB offset measurement is ob- 
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Fig. 21. Stability of the flux measurements during the calibra- 
tion sequences in gMfM>-M bands. The difference of instrumental 
magnitudes in an aperture of 43 pixels between the first and the 
last exposure of the sequence is displayed as a function of the 
seeing evolution between the exposures. 



Table 6. Determination of the calibration offsets for the SNLS 
catalogs. 





Um 


8m 


rtt 


iti 


i2 M " 


Zm 


BD+17 


<5ab 


-0.439 


0.114 


-0.123 


-0.332 


-0.316 


-0.458 




0.024 


0.005 


0.006 


0.005 


0.005 


0.019 


SA 


<5ab 


-0.431 


0.111 


-0.137 




-0.335 


-0.447 




0.017 


0.005 


0.008 




0.007 


0.011 


^ 2 /d.o.f. 


8.4/8 


15/14 


8.1/9 


0/0 


10/15 


3/8 


A d 


0.008 


-0.003 


-0.013 




-0.019 


0.011 



Notes. 

For comparison purposes, cross-calibration of the i M and i2 M 
catalogs was used to provide a pseudo-measurement of BD +17 4708 
in il M (see Sect. 5.3) 

{b> Errors quoted from the systematic error budget (cf. Table 5). 

(c) Estimate of the uncertainty affecting the combination of solar 

analog measurements. The model includes the (potentially correlated) 

contributions from measurement noise, zero point determination 

uncertainty, flat-field noise, atmospheric variation, airmass correction, 

survey uniformity and aperture correction uncertainty. 

(rf) Difference between the two determinations of the offset. 



tained as the least-square average of the individual measure- 
ments 6^: 

i i 

where cr^ b = 1/ 2 (/ Ry 1 is the variance of the final estimate. 

The results are summarized in table 6. The offset can be 
compared to the offset inferred from Landolt observations of 
BD +17 4708. The results are further discussed in Sect. 9 where 
they are compared and combined with the other calibration mea- 
surements. 

5.3. Intercalibration ofi M and i2 M catalogs 

The number of epochs with coincident observations of Landolt 
and SNLS fields in band Um is small. Consequently, no accurate 
calibration of the i2 M band measurements can be obtained in 
this manner. On the other hand, no direct measurements of the 
HST standards were obtained before the breaking of the filter. 



Even so, a relatively precise cross-calibration of the Im and Um 
catalogs is possible due to the similarity of the two filters. 

We measured the average transformation between i\ Xo and 
/2| Xo for stars in the color range < g\ Xo - i\ Xo < 1 encompassing 
the BD +17 4708 color (g - i — 0.34). The measured color term 
between the two bands is i2\ Xg — i\ Xo = 0.025(+0.001)x(g| IO -i| J i :o ). 
The dispersion around this linear transformation is measured to 
be 4 mmag for stars brighter than im < 18. This is compara- 
ble to the expected measurement error of 3.5 mmag, dominated 
by the flat-field induced error. Therefore, intrinsic dispersion is 
expected to be no larger than 2 mmag. Applying this transfor- 
mation to BD+17 allows us to infer the Um magnitude of the 
primary standard from its gM and in magnitudes with a preci- 
sion better than 2 mmag, thus anchoring together the catalog flux 
scales in the Um and im bands. 

6. SDSS tertiary star catalog 

The SDSS follows a calibration scheme similar to what was pre- 
sented for the SNLS. To calibrate individual supernova measure- 
ments, a high-accuracy photometric catalog of non-variable stars 
was generated from a coaddition of repeated observations of 
stripe 82. The construction of this catalog is described in Ivezic 
et al. (2007). For brevity, this star catalog will be referred to as 
the "coadd" catalog. 

A comparison of MegaCam Stripe 82 observations (de- 
scribed in Sect. 8) and the coadd catalog revealed evidence for 
a percent-level residual non-uniformity in the final flat-fielding 
solution of the coadd catalog. While the effect on individual SN 
photometry is negligible and the errors tend to average out when 
cosmological parameters are determined, the non-uniformity 
triggered investigations of the origin of the flat-fielding error and 
its possible impact on calibration. We start this section by sum- 
marizing the key steps of the construction of the coadd catalog 
and then describe a revision of the coadd flat-fielding. 

6.1. Construction of the SDSS coadd catalog 

The coadd catalog was obtained by averaging SDSS measure- 
ments of 58 photometric runs from Stripe 82 delivering on av- 
erage 10 measurements per star. Quality cuts based on the re- 
peatability of the measurements were applied to discard variable 
stars. The version used in the SN processing consists of 681301 
stars. 

The flat-fielding of the SDSS camera is simplified by the drift 
scanning technique used to obtain imaging data. In drift scan- 
ning, each point on the sky is sampled by each CCD row; there- 
fore, the effective response is averaged over all rows. However, 
row averaging does not eliminate spatial and temporal variations 
of the atmospheric transparency. Long term drift is controlled by 
anchoring the catalog photometry to the PT secondary patches 
placed throughout the SDSS survey area. Stripe 82 has an un- 
usually high density of these secondary patches: on average, one 
every 4 degrees in RA but they are not uniformly spaced and 
there is a gap of 29 degrees between RA=60 and RA=89. 

Averaging over photometric runs improves measurement 
noise and several sources of systematics, most notably variations 
in atmospheric transparency. However, declination pointings re- 
main similar from run to run, such that Stripe 82 stars are only 
observed in a narrow range of pixel columns. Therefore, cross- 
scan variations (including flat-fielding errors and passband dif- 
ferences) are not smoothed out by run averaging. 

The original technique for flat-fielding the images based on 
sky levels proved to be problematic for reasons similar to those 
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preventing the direct use of twilights in MegaCam. A new pro- 
cedure that determined flat-field vectors for the SDSS camera 
based on the Ivezic et al. (2004) stellar locus was then imple- 
mented. 19 This technique stems from the observation that most 
of the stellar population lines up in a narrow locus in the 4D 
ugriz color-space. The stellar colors are projected onto 4 color- 
color diagrams where the bulk of the stellar population is tightly 
clustered around a line in color-color space. The requirement 
that the stellar locus be independent of the position on the focal 
plane can be used to determine differences in flat-fields between 
the bands. 

To apply the stellar locus technique to the Stripe 82 obser- 
vations, average principal colors have been computed in bins 
of 0.01 degrees in declination, over a region spanning Galactic 
coordinates. It is assumed that all bins probe mostly the same 
population of stars, thanks to the narrow range of declination of 
stripe 82, so that differences in principal color are not induced by 
changes in the stellar population. The flat-field differences were 
adjusted to ensure constant principal colors across declination 
bins. 

The technique only determines the difference in flat field 
between different filter bands. To determine the overall "gray" 
scale as a function of declination, Ivezic et al. (2007) relied on 
the calibration provided by the PT. The 2.5 m photometry in 
band gri is compared to the PT photometry of secondary stars, 
and the average difference is used as the reference flat-field vec- 
tor for r-band. In this procedure, any error in the PT flat-fielding 
is transferred to the 2.5 meter telescope as an error common to 
all the passbands. 

One last correction was applied to improve the uniformity 
of the catalog by taking into account the differences in the filter 
passbands for the different camera columns. In the published ver- 
sion, 20 all stellar measurements were transformed to the system 
defined by the official SDSS transmission curves 21 . The present 
work is based on an earlier (unpublished) version of the catalog 
that was used by (Holtzman et al. 2008) to calibrate the photom- 
etry of the supernova light curves. In this earlier version, no color 
transformation is attempted and all magnitudes are delivered in 
the natural system. 

The coadd catalog photometric scale is anchored to the PT 
secondary patches by the SDSS pipeline. The PT photometry, in 
turn, is tied to the USNO star network Smith et al. (2002). The 
calibration of the secondary patches is based on the PT obser- 
vations of 3 HST CALSPEC standards (Holtzman et al. 2008) 
so we do not rely on the calibration of the USNO star network. 
We do rely on its uniformity and suitability to provide an ac- 
curate measurement of atmospheric extinction. Observations are 
transformed from the PT natural system to the 2.5 m natural sys- 
tem using canonical color transformations. Standard star obser- 
vations have been conducted with the target object roughly at the 
center of the camera image, whereas the color transformations 
have been derived from measurements over the entire image. In 
this observational setup, non-uniformity in the PT photometry 
may result in some calibration bias. In addition, the color trans- 
formations were determined quite early in the survey and were 
not monitored closely. 



See http://www.sdss.org/dr7/algorithms/flatfield. 
html for a detailed description of the DR7 flat-fielding. 

20 http : / /www . astro . Washington . edu/user s/ivezic/ sdss/ 
catalogs/stripe82calibStars_v2 . 6 . dat . gz 

21 available from http://www.sdss.org/dr5/instruments/ 
imager/index.html#filters 



As variations of the PT response as large as 2% were found 
(cf. appendix C), the following sections describe our determi- 
nation of a flat-fielding correction to the coadd catalog which 
removes most of the sensitivity to PT photometry. In Sect. 7, we 
will use this uniform catalog to determine corrections to the PT 
flat-field, update the PT to 2.5 m transformation equation and 
revise the AB offsets for the corrected catalog from PT observa- 
tions of HST photometric standards. 

6.2. Corrected SDSS flat-field 

A separate effort (Padmanabhan et al. 2008) sought to enforce 
the uniformity of the SDSS survey with naturally overlapping 
fields and special crossing scans that observed the sky in di- 
rections that were approximately perpendicular to the SDSS 
Stripes. The observation of stars in different camera columns 
allowed a recalibration of the SDSS that did not rely directly 
on the PT flat-fielding. In addition, it calibrated each filter sepa- 
rately and did not rely on the stellar locus or other assumptions 
about the objects being observed. This technique was used for 
the SDSS Data Release 8 22 . This catalog will be referred to as 
DR8 and was retrieved from the SDSS database using the query 
given in appendix D. 

The DR8 catalog is based on a single observation of each 
star. Consequently it is statistically less precise than the coadd 
catalog. In addition, typical SDSS observations take place at low 
airmass making it difficult to determine the atmospheric extinc- 
tion from the SDSS 2.5 m data alone, requiring some additional 
assumptions. In contrast, the PT routinely measures standard 
stars at both low and high airmass and measures atmospheric 
extinction with much better accuracy. 

Since we know that the PT flat-field has significant errors 
we want to use the DR8 data to correct the coadd catalog as 
a function of declination. We correct the coadd catalog, rather 
than using DR8 directly, because the multiple star observations 
incorporated into the coadd catalog greatly decrease the photo- 
metric errors and make it possible to eliminate various types of 
outliers. In the following, we will refer to the resulting catalog 
as the "corrected" catalog. Corresponding magnitudes will be 
denoted with a c subscript. 

6.2.1. Matching the coadd and DR8 catalogs 

We analyzed the differences between the two catalogs. The re- 
sults of the analysis are restricted to the resulting catalog of 
matched stars. The matching criteria were: 

- The separation of two objects in the coadd and DR8 catalogs 
is less than 1". 

- There is no other object in the coadd catalog that is within 
3.6" in both right ascension and declination. 

- The object in the coadd catalog is not within 3.6"in both 
right ascension and declination of any other object in the 
DR8 catalog. 

A star must have 4 or more observations in the g, r, and i bands 
to be included in the catalog. In the following analysis, the mea- 
sured u and z magnitudes are used only if there are 4 or more 
observations. We have also chosen to require 14.0 < r < 20.4 



http : //www . sdss3 . org/dr8/. 
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Fig. 22. The difference in magnitude (coadd-DR8) for the 5 
SDSS filters is shown as a function of right ascension. 



6.2.2. Right ascension trend 

Figure 22 shows the difference in the two catalogs as a function 
of right ascension. The data show a clear trend that is approxi- 
mately linear in right ascension (a). We fit the data for each filter 
to the form 

Am — aa + b (22) 

The results of the fits for each filter are shown in Table 7. 

Given the significant trend in the difference, the question 
arises as to which catalog is more nearly correct. The SDSS DR8 
procedure (Padmanabhan et al. 2008) assumes the extinction to 
be of the form 

dk 

k(t) =ko + —{t- t ) (23) 
at 

where Icq and dk/dt are constants, t is the time of observation and 
fo a reference time (in this case 0700UT). The assumed mean 
values and dispersion of dk/dt are taken from Padmanabhan 
et al. (2008), converted into mmag/deg and shown in the "DR8" 
column of Table 7 assuming an airmass of 1.2, a typical airmass 
for a Stripe 82 run. The slope of the difference between the coadd 
and DR8 are similar to the values assumed by DR8, but not iden- 
tical. 

The mtPipe calibration (Tucker et al. 2006) used for SDSS 
data production does not include any time-dependent terms al- 
though occasionally different photometric solutions were used 
for different time intervals during the night. The PT calibration 
is imprinted on the SDSS 2.5 m data through the secondary 
patches. Thus, we would expect a time variation-induced cal- 
ibration error in the SDSS catalog data with the effect being 
proportional, on average, to the difference in time between the 
2.5 m observations and the PT observations. A typical scan of 
Stripe 82 starts in the early evening at approximately -60° and 
ends in the early morning at +60°, resulting in a tight correlation 
between RA and the time of night. In contrast, the 178 Stripe 82 
secondary PT patches used to calibrate the coadd catalog were 
taken throughout the night at more-or-less random times. 

Given the number of patches and their distribution in time 
of observation, it seems likely that the secondary patches should 
provide a good average over observing conditions throughout 
the night. Seasonal variations could produce a bias in calibra- 
tion as a function of right ascension since higher right-ascension 
patches tend to be observed later in the season, but any variation 
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Table 7. Trend line fits in RA 



Filter 


DKo - 


coad 




PT - coadd 




a 


b 


dk/dt 


a 




(mmag/deg) 


(mmag) 


(mmag/deg) 


(mmag/deg) 


u 


0.142 


3.1 


0.096 ± 0.200 


0.082 ±0.291 


g 


0.082 


3.7 


0.056 ±0.136 


-0.006 ± 0.027 


r 


0.075 


-3.4 


0.080 + 0.136 


0.013 ±0.017 


i 


0.087 


-6.6 


0.096 ±0.120 


0.004 ± 0.024 


z 


0.199 


-16.7 


0.176 ±0.136 


-0.022 ± 0.059 



Notes. The comparison with Padmanabhan et al. (2008) assumes an 
airmass=1.2, which is typical for a Stripe 82 run. 



that is slower than a night should be removed by the nightly cal- 
ibration. We therefore assume that the PT patches have no bias 
as a function of right ascension. 

In the processing of the 2.5 m data a single zeropoint is com- 
puted for each CCD so changes in atmospheric extinction could 
result in a bias in calibration as a function of right ascension. 
On the other hand, a typical photometric run in the stripe 82 cat- 
alog is often much shorter than the maximum of 8 hours: the 
average run length is 3 hours. Thus the effect of any change in 
atmospheric extinction should be substantially reduced for the 
shorter run. Regardless of what we might expect from the SDSS 
photometric processing, we can examine the agreement between 
the coadd catalog and the PT patches, which we assume to be 
unbiased with respect to right ascension. For each PT patch we 
determine the average difference between coadd and PT obser- 
vations of the same stars. We then fit a line to the data as a func- 
tion of RA. The results are shown in table 7. The individual PT 
patches show considerably more scatter than would be expected 
from the statistical errors; presumably this is due to small time 
variations in atmospheric extinction (which was assumed to be 
constant over the course of the night). We characterize the distri- 
butions by fitting straight lines whose slopes are reported in the 
last column of table 7. The quoted uncertainties are computed 
from the scatter of the data, not the statistical errors on the in- 
dividual points. Also given in table 7 are the slopes expected if 
the DR8 trend were correct. The fitted slopes are consistent with 
zero and, for all bands except u, are significantly different from 
the slopes suggested by the DR8 catalog. We conclude that this 
comparison disfavors the trends of the DR8 catalog and we re- 
move the slopes reported in Table 7 from the DR8 data, bringing 
it into agreement with the coadd. 

There are remaining trends in the data that are not well fit 
by the straight-line approximation. The origin of these variations 
is likely to be fluctuations in atmospheric transparency which 
affect the DR8 single epoch photometry more than the coadd 
catalog. However, we are primarily interested in the flat-fielding 
as a function of declination, and the major effect of remaining 
trends on the flat-fielding is a slight increase of the dispersion. 

6.2.3. Flat-field correction 

Figure 23 shows the comparison between the DR8 catalog and 
the coadd as a function of declination after the linear trend in 
RA has been removed from DR8. The data are binned into 226 
declination ranges (each bin is approximately 0.01° wide). The 
catalog differences display similar patterns from band to band 
(particularly in r, i, and z), reinforcing the notion that the PT 
flat-field is a significant source of error in the coadd catalog. 
However, the bands also exhibit significant trends and the pattern 
of residuals does not show the exact 0.6 degree periodicity that 
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would be expected if a gray scale error in the PT were the only 
effect. 

The solid lines on Fig. 23 are cubic B-spline curves, each 
of which has 7 1 control points determined by a fit to the binned 
data. The decision to use 71 control points was somewhat arbi- 
trary but puts the angular resolution at about 0.03. The analysis 
of the PT flat field leads us to expect a sharp discontinuity where 
the secondary patches abut one another. This would argue for a 
large number of control points in order to adjust those expected 
discontinuities. However, it was suspected that the binning of 
data to produce the coadd originally may have introduced some 
numerical noise, so the number of control points was chosen to 
be smaller than the number of data bins in order to produce some 
smoothing effect. The x 1 of the fit decreases slowly as the num- 
ber of bins is increased and doesn't seem to offer any useful clues 
as to an appropriate cutoff. 

We correct the coadd catalog by subtracting the fitted smooth 
correction. Because we redetermine a (small) absolute calibra- 
tion offset for the corrected catalog (described in the next sec- 
tion), we do not require the average correction to be zero. The 
residual difference between the DR8 and the corrected coadd 
catalogs shows a rapid 2 mmag rms fluctuation on small angular 
scales but is constant over larger scales. An assessment of the 
flat-fielding quality of the corrected catalog is given by the com- 
parison with MegaCam measurements in Sect. 8.4. In particular, 
Fig. 28 suggests that the flat-fielding is now correct at the 0.5% 
level in griz. 

7. SDSS PT calibration measurements 

We now turn to the determination of AB offsets for the corrected 
coadd catalog using SDSS PT observations of HST standards. 
The non-uniformity of the PT response has a potential impact on 
this calibration transfer and must be corrected for. The basic idea 
is to use the uniformity of the corrected coadd catalog to improve 
the flat-field of the PT, redetermine the color relation between 
the corrected coadd catalog and the corrected PT photometry, 
and use those corrected relations to transform PT observations 
of the primary standards to the 2.5 m photometric system. 

7.1. Correction of PT measurements 

The calibration of the 2.5 m telescope native system to the ab- 
solute (AB) system relies on repeated SDSS PT observations of 
the HST CALSPEC 23 standard stars (Bohlin & Cohen 2008) . 
The PT measurements are calibrated against the set of USNO 
standard stars and transformed to the SDSS 2.5 m natural sys- 
tem using fixed color transformations. The PT calibration is ro- 
bustly transferred to the 2.5 m telescope through the PT sec- 
ondary patches. 

With the discovery of photometric errors in the PT, it be- 
comes important to redetermine (or check) the transformation 
equations between the two systems. We did so by looking for 
corrections to the already transformed magnitudes taking the 
form: 



u c = u + h u (u — g — 1.362) + ci u (24) 

8c = g + b g (g-r -0.596) +a g (25) 

r c = r + b r (r - i - 0.241) + a r (26) 

i c = i + bi(r-i -0.250) + ai (27) 

Zc = z + b z (i -z- 0.128) + a z , (28) 



http://www.stsci.edu/hst/observatory/cdbs/calspec.html. 



Table 8. Color range used to compute the color transformations. 



Color 


Minimum 


Maximum 


u-g 


0.7 


2.7 


g~r 


0.15 


1.2 


r — i 


-0.1 


0.6 


i - z 


-0.2 


0.4 



Table 9. PT Color Transforms 



Parameter 


u 


8 


r 


i 


z 


a 


-0.0022 


-0.0021 


-0.0022 


-0.0014 


-0.0078 




0.0011 


0.0003 


0.0003 


0.0003 


0.0006 


b 


-0.0016 


-0.0040 


-0.0011 


-0.0145 


-0.0255 




0.0041 


0.0018 


0.0027 


0.0030 


0.0067 


x 1 


3739 


10003 


11427 


10820 


6215 


d.o.f. 


3012 


10100 


12215 


11391 


6636 



Here we denote the magnitude coming from the corrected coadd 
catalog with the subscript m c , to avoid confusion with the PT 
measurements transformed to the ugriz using the nominal color 
transformations. 

As we are fitting for corrections, we would recover b = a = 
in all filters if the nominal transformation described the data per- 
fectly. In order to compute these transformations the PT observa- 
tions of secondary stars are matched with the corrected catalog 
of 2.5 m observations with the following criteria: 

- Match within 0.001 degrees 

- Matches must be unique. 

Since we know that the PT is subject to flat-fielding errors 
we want to allow for the possibility that a varies with position 
on the focal plane. The color term b, however, is likely to be 
nearly constant. So, we fit our data to equations (24-28) in two 
steps. First, we determine the color term using the central portion 
of the PT CCD and the following criteria: 

- Valid PT magnitude and SDSS magnitudes for color 

- At least 4 observations in the corrected catalog 

- 2.5 m r-band magnitude in the range 14.0 < r < 19.0 

- The star must be in the center of the PT focal plane (400 < 
x < 1448 and 400 <y < 1448), where x and y are the pixel 
coordinates of the 2048 x2048 pixel CCD . 

- Color in the range given by Table 8. 

The fit is iterated and any magnitude that deviates by more than 
5cr is dropped from the fit. The data are shown in Figure 24 and 
the results of the fits are given in Table 9. In general, the color 
terms and offsets between the two systems are quite small. While 
some of the b terms are statistically significant their magnitudes 
are small. 

Next, we determine a position dependent a term by fitting 
the data in bins of a uniformly spaced 7x7 grid across the CCD 
focal plane. Although the b terms are not significant, the data 
is corrected for color prior to determining position dependent a 
terms. The resulting pattern of a terms is shown in Fig. 25. The 
main feature is a ~ 4% gradient between the 2 opposite corner 
of the CCD and has a similar shape in all the bands. 

7.2. Absolute calibration 

We can now use the PT measurements to update the SDSS su- 
pernova survey calibration. The original HST standard star PT 
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Fig. 23. Difference between the coadd catalog and DR8 as a function of declination. The symbols represent the binned data and the 
lines are a cubic spline function fit to the data with 71 control points. 




Fig. 24. SDSS 2.5 m to PT measured color transformations. The individual stellar measurements are averaged in color bins of 0.05 
width for clarity. The error bars display the lcr uncertainty on the average estimated from the dispersion of individual measurements. 
The straight black line is the linear fit to the data. 



measurements are displayed in Table 10. They are reported, as 
usual, in the 2.5 m photometric system by the application of 
the nominal (uncorrected) color transformations. Table 12 lists 
the same PT measurements in the corrected catalog system, ob- 
tained by applying equations (24-28) to the original (Table 10) 
magnitudes. As described earlier in the text, these corrections 
are PT focal plane position-dependent, and compensate for the 



non-uniformity of the PT response. AB offsets calculated from 
the corrected HST standard measurements are given in Table 1 1 . 

HST spectra and SDSS filter responses (Doi et al. 2010) are 
used to compute expected magnitudes for the HST standard stars 
via the synthetic magnitude formula given in Eq. (1). Results of 
this calculation are given in Table 13. As the CALSPEC ver- 
sion somewhat affects the resulting magnitudes, Table 13 also 
notes the specific files used. Given in Table 14, the differences 
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Table 10. PT measurements of HST standard star magnitudes 



Star 


m u 


6m u 


n u 


m g 


6m g 


n g 


m r 


6m r 


n r 


m, 




n i 


m. 


8m, 


n. 


G191B2B 


11.045 


0.031 


7 


11.461 


0.011 


7 


12.014 


0.011 


7 


12.397 


0.010 


1 


12.743 


0.005 


1 


GD153 


12.693 


0.020 


6 


13.055 


0.018 


6 


13.570 


0.006 


6 


13.943 


0.010 


6 


14.296 


0.021 


6 


GD71 


12.439 


0.016 


6 


12.753 


0.011 


7 


13.242 


0.012 


7 


13.610 


0.005 


7 


13.967 


0.023 


7 


P041C 


13.569 


0.015 


11 


12.260 


0.010 


11 


11.842 


0.009 


11 


11.718 


0.024 


11 


11.702 


0.011 


11 


P177D 


15.120 


0.017 


13 


13.746 


0.010 


13 


13.299 


0.009 


13 


13.158 


0.008 


13 


13.126 


0.009 


13 


P330E 


14.551 


0.018 


17 


13.281 


0.009 


17 


12.842 


0.007 


17 


12.701 


0.007 


17 


12.672 


0.013 


17 


BD+17° 4708 


10.559 


0.020 


169 


9.630 


0.029 


169 


9.350 


0.027 


169 


9.248 


0.030 


169 


9.235 


0.032 


169 



Notes. The HST CALSPEC stars are measured in the PT native system but reported in the SDSS 2.5 m (unprimed) system using canonical color 
transformations. For each filter x, we report the average magnitude m x , the standard deviation (not the error on the mean) 6m x and the number of 
observations n x . 



Table 11. PT Focal Plane Offsets." 



Filter 


X 


y 


a u 


a g 


a r 


<3i 


a z 


G191B2B 


816 


1013 


0.0018 


0.0026 


0.0015 


0.0016 


0.0039 


GD71 


261 


1251 


0.0153 


0.0149 


0.0080 


0.0058 


0.0100 


GDI 53 


914 


1426 


-0.0006 


0.0021 


-0.0010 


-0.0008 


0.0022 


P041C 


977 


898 


0.0021 


0.0004 


0.0003 


0.0015 


-0.0003 


P177D 


712 


1239 


0.0040 


0.0050 


0.0018 


0.0013 


0.0038 


P330E 


876 


633 


0.0001 


0.0018 


0.0046 


0.0031 


0.0024 


BD+17° 4708 


1024 


1024 


0.0024 


-0.0001 


0.0002 


0.0018 


-0.0015 


Error 


1024 


1024 


0.0028 


0.0009 


0.0009 


0.0019 


0.0018 



Notes. Values of the a term from Eq. (24-28) at the position of the PT focal plane where the star is measured. See Fig. 25 for a map of the 
correction. 

Table 13. Synthetic Magnitudes of HST Standard Stars." 



Star 


u 


g 


r 


i 


z 


File 


G191B2B 


11.008 


11.471 


12.011 


12.398 


12.746 


g 1 9 lb2b_stisnic_003.ascii 


GDI 53 


12.677 


13.067 


13.580 


13.959 


14.300 


gd 1 5 3 _stisnic_003 .ascii 


GD71 


12.433 


12.770 


13.258 


13.627 


13.970 


gd71_stisnic_003.ascii 


P041C 


13.507 


12.281 


11.852 


11.747 


11.734 


p04 lc_stisnic_003 .ascii 


P177D 


15.053 


13.771 


13.308 


13.178 


13.143 


p 177d_stisnic_003 .ascii 


P330E 


14.471 


13.295 


12.840 


12.708 


12.675 


p330e_stisnic_003 .ascii 


BD 17°4708 


10.499 


9.647 


9.351 


9.255 


9.239 


bd_17d4708_stisnic_003.ascii 



Notes. 

(fl) The SDSS filters are from Doi et al. (2010). 
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Table 12. Standard star magnitudes in the SDSS system after all 
corrections. 



Star 


m„ 


m g 


m r 


nii 


m z 


G191B2B 


11.048 


11.456 


12.014 


12.388 


12.735 


GDI 53 


12.699 


13.051 


13.573 


13.936 


14.289 


GD71 


12.429 


12.736 


13.236 


13.597 


13.946 


P041C 


13.569 


12.261 


11.844 


11.716 


11.707 


P177D 


15.118 


13.743 


13.299 


13.157 


13.128 


P330E 


14.553 


13.280 


12.839 


12.697 


12.675 


BD17° 4708 


10.560 


9.631 


9.352 


9.245 


9.241 



Notes. The color transformation and position dependent offsets de- 
scribed in the text have been applied to the measurement listed in 
Table 10. 



Fig. 25. PT flat-field corrections. 



between the synthetic and the measured magnitudes are the AB 
offsets of the SDSS photometric system. The standard calibra- 



tion uses the solar analogs because their colors are well within 
the range of validity of the PT— >SDSS system color transforma- 
tion. The dispersion in the results of the 3 solar analogs gives 
an estimate of the error in the AB offset - about 0.002 to 0.003 
magnitudes - although this does not include any systematic ef- 
fects that might result in common errors to the solar analogs. 
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Table 14. AB Offsets for SDSS." 



Filter 


u 


8 


r 


i 


z 


Solar Average 
Solar Error 


0.069 
0.006 


-0.021 
0.004 


-0.006 
0.003 


-0.021 
0.006 


-0.014 
0.008 


BD+17° 4708 


0.061 


-0.016 


0.001 


-0.010 


0.002 


WD Average 
WD Error 


0.020 
0.013 


-0.021 
0.006 


-0.009 
0.007 


-0.021 
0.006 


-0.015 
0.004 



Notes. 

(fl) Offsets should be subtracted from the SDSS native magnitudes to 
obtain calibrated AB magnitudes. 



We also show BD +17 4708, which serves as the primary ref- 
erence for the SDSS photometry as well as for several other 
experiments. The difference between BD +17 4708 and the so- 
lar analogs is another indication of AB offset errors. The WD 
result is also shown. However, WD colors are far beyond the re- 
gion of validity for the color transformation thus their computed 
AB offsets should not be regarded as reliable. 

7.3. Discussion 

Besides the hints given by the variation of the calibration as a 
function of the selected standard, we can try to review and quan- 
tify the sources of uncertainty affecting the calibration of the 
corrected catalog. As for the SNLS, our concern here are the un- 
certainties affecting the measured magnitudes of the HST stan- 
dards. The discussion of the spectra uncertainties is common to 
both surveys and delayed to section 9. 

7.3.1. Survey uniformity 

The PT calibration can be affected by a spatial non-uniformity 
in the USNO star network. In particular, non-uniformity in the 
network could result in a systematic difference between the cali- 
bration of the solar analog measurements (a w 15h) and the cal- 
ibration of the secondary patches (20/i < a < Ah). The question 
is thus how much error there might be in transferring the calibra- 
tion 12 hours ahead (or behind) in RA. Possible non-uniformities 
are not directly discussed in the literature but the USNO network 
is expected to be especially robust because of nearly continuous 
measurements of the primary standards over a two year period of 
observation. From the uncertainties quoted in Smith et al. (2002) 
for the primary standard star measurements in the global fit, we 
estimate that the relative uncertainty in uniformity is less than 
1.5 mmag in u band and 1 mmag in the other bands. 

The PT observation program was designed so that the USNO 
standard stars were observed at the center of the focal plane. 
In contrast to the measurements of the HST standards, the non- 
uniformity of the response thus has no impact on the PT calibra- 
tion to the USNO standard star network. 

7.3.2. Color transformation 

The color transformation from the PT to the 2.5 m system in- 
troduces some uncertainty on the precision to which the 2.5 m 
magnitudes of the standard stars are known. Two effects must be 
considered: the fact that individual stars do not fall exactly on 
the average color transformation, and the fact that in the SDSS 
camera each CCD has its own filter. 

We found the intrinsic dispersion around the average trans- 
formation by looking at the dispersion of magnitudes synthe- 



sized in both systems from a stellar spectral atlas (Gunn-Stryker- 
Bruzal). We use this dispersion to estimate how much spectral 
variations of the HST standards might vary from the average 
stellar population which determines the color transformations. 
This dispersion is probably an overestimate because we ignore 
the effects of noise in the spectral atlas. In addition, the spectral 
diversity is likely over represented by the spectral atlas. 

An alternative estimate of the magnitude of this error can 
be obtained from the variation of the AB-offset with the spec- 
tral type and color of the standard considered (available from 
table 14). We obtain similar numbers, except for u band where 
the variation is inflated by the fact that the blue WDs clearly lie 
outside the validity range of the linear color transformation, and 
g band where the dispersion in the atlas seems artificially high. 
However, we do not attempt to correct the dispersion measured 
from the atlas and accept this dispersion as a conservative error 
estimate. This error dominates the uncertainty budget for the PT- 
based calibration as shown in table 16. This is similar to the sit- 
uation in SNLS, where the transformation of the BD +17 4708 
Landolt magnitudes dominates the uncertainty on the Landolt- 
based calibration. 

The color transformation is adjusted using data covering the 
entire width of the SDSS camera. As a consequence it translates 
PT measurements to the "average" SDSS photometric system. 
Our corrected tertiary catalog stems from a version of the coadd 
catalog with star magnitudes reported in the native SDSS photo- 
metric system. The flat-fielding correction applied to match the 
DR8 photometry does not use any color-dependent term, and 
leaves magnitudes in the natural system. This calibration proce- 
dure should adjust the zeropoint of each camera column so that a 
star of average color (as is nearly the case for solar analogs) will 
have the same magnitude in each camera column in the native 
system. 24 

The filters have been measured independently. Table 15 
shows the synthetic magnitudes for the 6 camera columns for 
one of the HST white dwarfs and one of the solar analogs. 
The differences in magnitudes are small (1%) but not com- 
pletely negligible. The variation of the AB magnitude of the so- 
lar analog around its nominal value illustrates the variation of the 
calibration offset effectively applicable to individual columns. 
Neglecting the filter changes adds a calibration noise (centered 
on zero) of 6, 4, 1, 0.2 and 0.6 mmag rms in ugri and z bands, 
perfectly negligible for SN measurements. 

As far as we are only interested in accurate calibration on 
average (which is the case for SN-Ia science), the global AB off- 
sets are adequate, and we do not account for any related effect 
in our systematic budget. When considering only a subset of the 
camera columns however, some attention should be paid to the 
change of photometric system, and calibration offsets for indi- 
vidual columns may have to be determined from table 15. 

8. SDSS/SNLS direct cross-calibration 

8.1. Cross-calibration dataset 

In theory, several overlapping areas between the CFHTLS and 
the SDSS (in particular the CFHTLS deep fields D2 and D3) 
could be used to tighten both surveys' photometric calibra- 
tion. However, the utility of existing published data sets is lim- 
ited by two factors: the precision of the SDSS zero point de- 
termination for single epoch photometry and the difficult-to- 
evaluate precision of the inter-calibration between the D2-D3 

24 Of course, there will be differences in magnitude for stars of differ- 
ent colors. 
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Table 15. Synthetic and natural SDSS magnitudes of HST standard stars by camera column 



Star 


column 


u 


8 


r 


i 


z 


G191B2B 


1 


11.0087 


11.4721 


12.0166 


12.4006 


12.7460 


G191B2B 


2 


11.0061 


11.4686 


12.0152 


12.4036 


12.7595 


G191B2B 


3 


11.0133 


11.4752 


12.0087 


12.4039 


12.7520 


G191B2B 


4 


11.0110 


11.4742 


12.0148 


12.4020 


12.7542 


G191B2B 


5 


11.0058 


11.4749 


12.0141 


12.4018 


12.7470 


G191B2B 


6 


11.0112 


11.4750 


12.0147 


12.4006 


12.7430 


G191B2B 


nom 


11.0094 


11.4733 


12.0140 


12.4021 


12.7504 


P041C 


1 


13.5112 


12.2836 


11.8512 


11.7471 


11.7339 


P041C 


2 


13.5185 


12.2888 


11.8518 


11.7467 


11.7327 


P041C 


3 


13.5025 


12.2767 


11.8547 


11.7466 


11.7333 


P041C 


4 


13.5067 


12.2786 


11.8519 


11.7468 


11.7330 


P041C 


5 


13.5183 


12.2776 


11.8522 


11.7469 


11.7339 


P041C 


6 


13.5073 


12.2774 


11.8520 


11.7470 


11.7346 


P041C 


nom 


13.5107 


12.2805 


11.8523 


11.7469 


11.7336 



Table 16. Systematic uncertainty budget for SDSS 



Source 


u 


g 


r 


i 


z 


PT measurement error" 


0.0026 


0.0015 


0.0013 


0.0025 


0.0018 


PT transformation to 2.5* 


0.0029 


0.0009 


0.0008 


0.0009 


0.0017 


Dispersion around the transformation 


0.008 


0.007 


0.003 


0.004 


0.004 


USNO uniformity 


0.0015 


0.0010 


0.0010 


0.0010 


0.0010 



Notes. (fl) This is based on the variance of the PT observations. 

{b) This is based on the computed statistical error. The error is dominated by the zeropoint offset between the PT and 2.5 m. 



DR8 and the Stripe 82 observations. In this work we rely in- 
stead on a specially designed observation program of Stripe 82 
with MegaCam. 

The original program, called MAPC (for MegaCam absolute 
calibration program), imaged the CFHTLS deep fields Dl and 
D4 and 2 fields in Stripe 82 (hereafter SDSS36 and SDSS326) 
in short "observing blocks". Each observing block (OB) was ob- 
served in about half an hour in two bands. Since switching filters 
with MegaCam takes about 2 minutes, all the OB targets were 
first observed in one band, then another filter was selected to re- 
observe all the targets. To enable the monitoring of extinction 
stability, each observing block begins and ends with an expo- 
sure of the same field. Observations extended from December 
2006 to February 2008. The cross-calibration data set was com- 
plemented by an early-201 1 program of joint observations of Dl 
and SDSS36 alone (hereafter OB_Dl_ext) with a reduced scope. 

In the original program, 5 successive exposures of Stripe 82 
were taken with a coarse dithering applied. This allowed us to 
map a slightly larger field (-0.8° to 0.8° in declination, cover- 
ing two thirds of the width of Stripe 82) while providing a cross 
check of the MegaCam flat-fielding. The position offsets applied 
to each exposure are given in table 19 and roughly correspond 
to half of the focal plane size. Short exposures on standard fields 
were also included in the original program to provide an alterna- 
tive determination of the MegaCam absolute calibration, but this 
method proved more sensitive to observation systematics than 
the one presented in Sect. 5.2. Both of these complications were 
dropped from the complementary program to shorten the observ- 
ing block as much as possible. The definition and specifics of 
observing blocks for the two programs can be found in table 18. 

From this data sample we derive Stripe 82 star catalogs in 
MegaCam natural magnitudes, calibrated on the SNLS deep 
fields that can be readily compared to the SDSS catalog. 



Table 17. Field definition in the MAPC program. 



Field 


RA 


DEC 


Dl 


2:26:00.0 


-04:30:00 


D4 


22:15:31.0 


-17:44:05 


SDSS36 


2:25:59.9 


-0:00:00 


SDSS326 


21:46:00.0 


-0:00:00 


SA95 


3:54:13.9 


+0:02:41 


HZ4 


3:55:14.2 


+9:39:30 


GD71 


5:52:22.2 


15:45:24 


SA 113 


21:41:57.0 


0:19:32 


Feigel 10 


23:19:53.1 


-5:17:55 


BD+284211 


21:51:04.6 


28:44:01 



Table 19. The dithering pattern in MAPC. 



obs 


<5ra 


<5dec 


1 


0.00' 


0.00' 


2 


-24' 


+18' 


3 


-24' 


-18' 


4 


+24' 


-18' 


5 


+24' 


+ 18' 



8.2. Calibrated MegaCam measurements of the Stripe 82 
fields 

Images from the cross-calibration data are processed using the 
standard SNLS photometric pipeline. Each resulting catalog is 
then matched (regarding astrometry) to the tertiary star catalog 
of the corresponding field, namely the SNLS tertiary star cata- 
log described in Sect. 4 for SNLS deep fields and the corrected 
coadd SDSS catalog described in Sect. 6 for Stripe 82 fields. A 
matching radius of 1 arcsec was used and a few other selection 
criteria were applied to the SDSS corrected catalog: 
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Table 18. Observing group description in the MAPC program and its extension. 



Name 




gM 


r M 






Zm 


Field 


^exp 


Comment 
















Dl 


120 


















SDSS36 


120 


5 dithered exposures taken in a row 


OB Dl 


2 


2 


5 


2 


2 


2 


SA95 


3 


out of focus 1.5 mm 














HZ4 

GD71 

Dl 


2 
2 

120 


out of focus 1.5 mm 
out of focus 1.5 mm 
















D4 


120 


















SDSS326 


120 


5 dithered exposures taken in a row 


OB_D4 


2 


2 


4 


2 


2 


3 


SA113 
Feigel 10 
BD+284211 
D4 


3 
2 
2 

120 


out of focus 1.5 mm 
out of focus 1.5 mm 
out of focus 1.5 mm 
















SDSS36 


120 




OB_Dl_ext 


2 


1 


3 





7 


2 


Dl 

SDSS36 


120 
120 





Notes. <a) This column and the following give the number of times each observing block was completed in each band. Measures were subsequently 
selected. 



- The number of z band observations of the star must be greater 
than 4 

- The rms of the gri band star measurements must be less than 
30 mmag. 

From studying MegaCam photometric flat-field corrections, 
we know that the accuracy of single epoch flat-fields are limited 
by errors that are expected to average out in the construction 
of the deep field tertiary star catalogs. A simple way to propa- 
gate this improvement to single epoch photometry is to align the 
photometry of deep field exposures on the corresponding tertiary 
star catalog. For each observing group we compute a rough flat- 
tening correction to the photometry by computing a zero point 
deviation per chip in addition to the global zero point. The fol- 
lowing model is adjusted to the deep field observations available 
in the observing group. 

m\ x — otadu — Z, = 6Z(ccd) 

where nt\ x designates the magnitude reported in the tertiary cat- 
alogs, otadu the flat-fielded instrumental magnitudes and Z, has 
been adjusted for each exposure in a first pass neglecting the 5Z 
term. The fit is done in flux and measurements are weighted ac- 
cording to the error model introduced in Sect. 4.1. The SZ(ccd) 
corrections are then applied to all images of the observing group. 

The procedure described in section 4.1 is applied to build 
averaged MegaCam catalogs of the 2 Southern Stripe fields. 
Separate catalogs are built for different pointings. The catalogs 
are then calibrated according to the zero points computed on 
deep field exposures. The uncertainty on the catalog zero point 
is estimated as the night-to-night dispersion of the fit residuals, 
divided by the square root of the number of nights. 

Direct comparison with the SDSS is enabled by the transfor- 
mation of the measurements to the SNLS uniform photometric 
system. This is done using the transformation derived in Sect. 3.2 
for stars lying in the valid color range of the transformation. 

This procedure results in two sets of roughly 2000 and 6000 
stars with natural SDSS magnitudes (denoted mxs) and uni- 
form SNLS magnitude (denoted m\ Xo ). Depending on the point- 
ing and photometric band between 2 and 9 nights are averaged 
(see Table 18). For such short observing blocks, the rms of the 
nightly zero points is typically below 5 mmag. 



8.3. Passband consistency 

The first outputs obtained from the comparison of SNLS and 
SDSS measurements of the Stripe 82 stars are precise color 
transformations between the two systems. We determined linear 
color-transformations between SDSS and SNLS measurements 
for stars in the color range 0.5 < g — i < 1.5 according to the 
model: 

m ]X0 -m 2 . 5 =a(g-i-0.8)+p (29) 

Measurements are weighted according to the errors reported 
in both catalogs added in quadrature. To account for residual 
variation in the catalog zero point, independent /3 are determined 
for each MegaCam pointing. The transformations are illustrated 
in Fig. 26 and the measured slopes are given in table 20. We 
emphasize that these equations describe transformations of the 
reference MegaCam system (filters at the center of the focal 
plane) to the SDSS system. The reference MegaCam system is 
not equivalent to and should not be confused with the average 
natural MegaCam system. 

These transformations can be compared to synthetic 
photometry-based predictions as a relative assessment of the ac- 
curacy of both instruments' effective passbands. We synthesized 
stellar magnitudes for the 2 photometric systems from 2 stel- 
lar libraries, and determined a linear relationship for the result- 
ing color-color plots. The slopes are given in table 20 and over- 
plotted in red on Fig. 26. They are to be compared with the fitted 
slope in black. We report the statistical uncertainty on the slope 
determined on the Gunn-Stryker-Bruzal library, and the differ- 
ence of slope obtained with respect to the Pickles library as an 
estimate of the systematic uncertainty on this quantity. 

The discrepancy between the expected and measured slopes 
can be used to estimate the uncertainty on the passband knowl- 
edge. We provide in table 20 the wavelength shift to be applied to 
the SNLS passbands to explain the whole effect. As can be seen 
the agreement is good in all bands, except in r M and i2 M where 
a 3<t shift is detected. Note that the um measurement does not 
provide useful constraints on the filter model. 

The effective SDSS passbands actually vary slightly from 
one chip to another. In order to investigate whether passband 
differences may account for some of the observed discrepan- 
cies, we subdivided the sample into declination slices roughly 
corresponding to the 4 SDSS camera columns covered by our 
measurements. We compare these results to slopes determined 
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Table 20. Transformation equation between SNLS and SDSS systems. 







u 




g 


r 






a 




z 


Expected color term 


-0.343 ± 


0.081 


-0.128 ± 


0.004 


-0.032 ± 0.001 


-0.025 1 


: 0.002 


-0.007 ± 0.002 


0.008 ± 


0.001 


Measured color term 


-0.341 ± 


0.002 


-0.126 ± 


0.000 


-0.026 ± 0.000 


-0.027 i 


: 0.001 


-0.003 ± 0.000 


0.009 ± 


0.001 


Apparent shift (nm) 


-0.1 


±4.9 


-0.2 


±0.5 


-2.2 ± 0.5 


1. 


8± 1.7 


-3.0 ± 1.2 


-0.5 


±0.7 




0.04 
0.02 
0.00 
-0.02 
-0.04 



1 1 


1 1 1 _ 


^ 


— 






1 1 


1 1 1 " 



0.4 0.6 



0.8 
82.5 
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- '2.5 
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I 
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m 
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li 
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1 

H 
N 
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-0.02 






0.04 


1 1 


1 1 1 " 
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1 1 
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Fig. 26. Linear color transformations between the SDSS and 
SNLS uniform magnitude system. The black lines represent 
the best linear fit of the data, whereas the red lines are based 
on synthetic magnitudes of the Gunn-Stryker-Bruzal stellar li- 
brary and the SNLS and SDSS transmission models described 
in Sect. 2.2.1. 



separately for the camera columns using the independent pass- 
band determinations from Doi et al. (2010) in Fig. 27. There may 
be some correlation between the small observed and predicted 
slope variations from one camera column to another, but the flat- 
fielding uncertainty and other systematics make this comparison 
problematic. This is especially the case in um, where both the 
predicted and measured color terms are affected by large un- 
certainties, making the comparison inconclusive. However, the 
measurements consistently disagree with the models in tm and 
Um bands, indicating that SDSS camera column transmission 
variations cannot explain the discrepancies noticed in Table 20. 
It is more likely that these discrepancies originate from real de- 
ficiencies in the MegaPrime passband models. We further dis- 
cuss passband uncertainties and their impact on the calibration 
in Sect. 9. 



8.4. Uniformity of the photometry 

As they are subject to different systematics, a thorough assess- 
ment of the photometry uniformity can be delivered by the com- 
parison of the two measurement sets. Due to the drift-scanning 
of the SDSS observations, the uniformity of the SDSS photome- 
try is expected to be excellent at small scales in right ascension. 
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Fig. 27. Comparison between synthetic and measured stellar lo- 
cus slopes in the SNLS-SDSS color diagram for the 4 SDSS 
camcols covered in the observation program. 



In contrast, no specific direction is expected to be preferred in 
MegaCam observations. 

This study led to the discovery of the flat-field problem in 
the PT and triggered the recalibration of the SDSS coadd cata- 
log. The photometry residuals (after correction of the Stripe 82 
catalog) are shown as a function of right ascension, declination 
and MegaCam radius in Fig. 28. Errors in SDSS flat-fielding 
would be expected to be most prominently seen when plotted 
versus declination while SNLS flat-fielding errors could appear 
in any of the plots. The radial plots correspond to an approximate 
MegaCam symmetry and exhibit some residual structure in r and 
z bands. The rms of residuals in bins of right ascension are 4, 1, 
2, 3, 2 and 3 mmag in um 8m *m 'm '2m and zm- Those numbers 
agree with the expected large-scale flat-fielding accuracy of the 
SNLS tertiary catalogs. 

8.5. Average offsets 

We can now tie the SNLS and SDSS calibrations together. 
Ultimately, the quantity of interest is the flux interpretation of 
a calibrated measurement. As filters themselves have uncertain- 
ties, the error on this quantity depends on the SED of the object 
itself (to first order in its color). We discuss this question in the 
next section. We display in table 21 the offset between the two 
surveys for a star of color g — i = 0.8 (defined as the /3 term in 
Eq. (29)). 
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Table 21. ft values as defined by Eq. 29 between SDSS and 
SNLS surveys. 



SDSS36 SDSS326 Combined 



u 


-0.432 


± 0.002 


-0.447 


± 0.002 


-0.438 


± 0.001 


±0.014 


g 


-0.082 


± 0.002 


-0.074 


± 0.001 


-0.076 


± 0.001 


± 0.008 


r 


-0.017 


±0.001 


-0.016 


±0.001 


-0.017 


±0.001 


± 0.001 


i 


-0.001 


± 0.001 


-0.003 


± 0.001 


-0.002 


± 0.000 


± 0.002 


i2 


0.010 = 


= 0.003 


0.012 = 


= 0.002 


0.012 = 


= 0.001 = 


t 0.003 


z 


0.030 d 


= 0.001 


0.035 = 


= 0.002 


0.031 = 


= 0.001 = 


t 0.005 



Notes. All numbers are expressed in magnitudes. Those numbers should 
not be interpreted as calibration offsets, part of the difference is expected 
from the definition of the magnitude systems. A discussion about AB 
offsets can be found in Section 9. 



The offset is computed separately on each field. The statis- 
tical uncertainty is estimated from the night-to-night dispersion 
of measurements divided by the square root of the number of 
observation nights. We then combine the two numbers and use 
the difference between the two fields as an estimate of the uncer- 
tainty on the mean related to the residual non-uniformity of the 
surveys. 

9. Combined calibration 

The calibration products for both surveys consist of: 

- measured natural magnitudes for a set of tertiary standards 
m 3 , in an arbitrary magnitude system. 

- magnitudes of primary spectrophotometric standards m s in 
the same magnitude system, affected by some measurement 
and systematic errors n s . 

The primary spectrophotometric standards are to be used to set 
the actual flux scale (equivalently determining a magnitude off- 
set to the AB system <5 a b) of the first by equating the mea- 
sured magnitude of the standard m s , with its expected magnitude 
m s synthesized from the spectrophotometric reference measure- 
ment: 

m s = m s + 5- db . (30) 

with: 

[ 1 AT(A)S s (A)dA 

m s = -2.5 log 10 f 

J A AT(A)S ab (A)dA 

The result is affected by uncertainties on the effective instrument 
passbands T and on the SED S s of the reference, as well as the 
error on the measured magnitude n s . 
In our case we have at our disposal: 

- 59 direct measurements of 3 solar-analogs (P330E, P177D 
and SNAP2) performed with MegaCam and calibrated to the 
SNLS system (see Sect. 5.2). 

- Landolt & Uomoto (2007) measurements of the star 
BD +17 4708, color transformed to the SNLS system in 5 
bands (see Sect. 5.1). 

- PT measurements of 4 photometric standards (BD +17 4708, 
P330E, P177D, P041C) color-transformed to the 2.5 m sys- 
tem in the 5 photometric bands (see Sect. 7). We do not use 
the PT observations of bluer standards as they are outside the 
valid range of the PT— >SDSS system color transformation. 

In addition, the cross-calibration data described in Sect. 8 yields 
6 constraints on the relation between the two surveys' tertiary 
catalog flux scales, and a straightforward constraint on the iu 
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and Hm magnitudes is given by the comparison of tertiary star 
magnitudes in these two bands (see Sect. 5.3). Altogether, we 
have a total of 91 measurements (or pseudo-measurements) and 
associated uncertainty estimates with which to constrain the AB 
offsets of the tertiary catalogs in the 1 1 photometric bands con- 
sidered. 

Our purpose here is to optimally combine the entire set of 
calibration data, to check its consistency, and to derive a reli- 
able estimate of the uncertainty on the resulting calibration. As 
we want to combine several standard stars measured in different 
photometric systems related by linear transformations, we start 
by discussing and parametrizing the uncertainty on the instru- 
ment passbands and primary standard SEDs. 

9. 1 . Uncertainties on passbands 

The uncertainty on the effective instrument passband comes 
from various sources and may take some arbitrarily complex 
form. For instance, the portion of passband uncertainty from op- 
tical elements whose transmissions can be measured should ap- 
proximately reduce to the uncertainty of the wavelength calibra- 
tion of the monochromator (i.e. a global wavelength shift). On 
the other hand, the portion of passband uncertainty due to other 
elements, like the CCD QE, can be harder to characterize. The 
determination of the average atmospheric transmission is also 
affected by uncertainties, which primarily impact the u band cal- 
ibration. In what follows, we assume that we can parametrize 
the error affecting the effective passband defining the photomet- 
ric band b by a shift in the passband's mean wavelength SA b eff . 
In other words, when performing first order computations in the 
next section, we will assume that passband uncertainties result 
from a global shift of the filter. We now attempt to derive an esti- 
mate of the amplitude of the uncertainty on the mean wavelength 

9.1.1. SDSS passbands 

Measurement errors on the SDSS passbands are not completely 
characterized by Doi et al. (2010), but we will estimate the un- 
certainty from data in that paper. The accuracy of the monochro- 
mator calibration is given as ±3A, which we assume is a con- 
stant error for all wavelengths in a given filter band. The ability 
to reproduce an independent quantum efficiency measurement 
(Figure 2 in Doi et al. 2010) suggests that the amplitude error is 
less than 2%. Assuming a +2% error that is linear in wavelength 
across the band, we calculate an error in the mean filter wave- 
length that is 0.3% of the filter bandwidth. Using 1500A for the 
filter width and combining that uncertainty with the monochro- 
mator uncertainty, results in a 6A uncertainty in mean filter re- 
sponse. The 1500 A width is slightly too large for g, r, and i bands 
and a factor of 3 too large for u and z. However, we don't de- 
crease the error estimate for these bands because there are prob- 
ably other (unquantified) uncertainties for these bands (e.g. the 
difficulty of determining the UV transmission and the red cutoff 
in z). 

9.1.2. SNLS passbands 

The lack of monitoring for the MegaCam passbands proved to 
be a weak point in the SNLS survey. It was only recently re- 
alized that two of the interference filters may have undergone 
significant alteration in their early days. 
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The REOSC filters were measured twice (in 2002 and 2006) 
and both measurements appear to be in rough agreement in all 
bands but iu and vm- Two arguments support the hypothesis that 
the 2006 measurements yield a better representation of the sur- 
vey passbands: first, the stability of the color locus of tertiary 
stars during the 5 years of the survey, presented in Sect. 3.4.1, 
excludes the hypothesis of a slow aging of the filters between 
2003 and 2008. Second, recent lab measurements of pieces of 
the broken im filter were found to be in agreement with the 2006 
measurement while incompatible with the 2002 ones. Because 
the 2006 measurements only cover the outer edge of the filters 
(12-15 cm away from the center), it was not possible to use them 
directly. Instead, we opt for applying a correction to the 2002 
im and rM measurements (see appendix B). The unaltered 2002 
measurements are used as-is for the gM, and zm filters. We 
now try to quantify the size of the uncertainty on the mean pass- 
bands resulting from these choices. 

The mean wavelength shift between the 2002 and 2006 mea- 
surements provides an estimate of the error at a specific position 
(about 15cm away from the center). The measured differences 
amount to 5, 3, -37, -31 and 6 A in u M , gM, I'm, im and zm 
bands respectively. The spatial consistency of the filter model is 
then ensured by the agreement of modeled and measured color 
terms between different positions (as illustrated on Fig. 7). A 
noticeable exception is the tm filter for which the measured and 
modeled color terms start to disagree at radii greater than 10cm. 
Without proper measurements of the filter at radius smaller 
than 10 cm, it is thus not possible to trust the corrected filter 
model more than the uncorrected model. Therefore, we assign 
an uncertainty of 37A to the filter passband. The situation 
appears a bit better for the i M filter, since the accuracy of the 
2006 measurement has been checked (cf. appendix B), and the 
agreement between predicted and measured color-terms is quite 
good. Nevertheless, it seems premature to assign an uncertainty 
smaller than 3 1 A to the iu passband without better understand- 
ing the source of the filter change. For the u M , g M , and zm filters, 
we also use the quoted differences between the two available 
measurements, namely 5, 3 and 6 A, as an estimate of the mea- 
surement uncertainty. 

The Hm BARR filter was never completely measured. In par- 
ticular its effective refraction index is unknown and assumed 
to be similar to that of the i M filter. We base our estimate of 
the Hm passband uncertainty on the stellar locus in the plane 
Hm — i vs g - i which has been well-constrained by observa- 
tions of Stripe 82 stars and synthetic photometry predictions (see 
Sect. 8.3). The comparison between synthetic predictions and 
measurements suggests an Hm shift of -30 A(Table 20). 

The uncertainty on the average atmospheric extinction curve 
is estimated in Buton et al. (2012) to be smaller than 0.02 
mag/airmass. This has a negligible impact on the definition of 
any of the passbands. A noteworthy exception is the zm filter, 
which covers a region of strong H2O absorption (8916-9929 A). 
Errors in the modeling of the average extinction could cause 
the mean wavelength of the effective filter to shift by a few 
angstroms. However, the error is unlikely to exceed the uncer- 
tainty quoted for this filter. Given their high throughput and the 
smoothness of their transmission, we do not expect significant 
contributions to a wavelength shift from any of the instrument 
elements besides the interference filters. 



9.2. SED of the primary standard stars 

Five spectrophotometric standard stars from the CALSPEC 
database are being considered in our final calibration data: the F 
sub-dwarf BD +17 4708 and the 4 G-type stars P041C, P177D, 
P330E and SNAP2. All these stars have STIS spectra calibrated 
by the model spectra of the white dwarfs G191B2B, GD153, and 
GD7 1 . Any error in the spectral modeling of those three funda- 
mental standard stars would equally affect our 5 standard stars. 

The calibration reference is affected by both the white dwarf 
atmosphere model error and uncertainties in the determination 
of the primary standard surface gravity and effective temperature 
from observations of their Balmer line profiles. Following R09, 
we base our estimate on considerations from Bohlin & Hartig 
(2002) and adopt a 0.5% slope uncertainty (1 <x) over the range 
3000 A- 10000A. Relative to g band 25 , the resulting (correlated) 
uncertainties are 0.7, 1.1, 2.2 and 3.2 mmag in u, r, i and z bands 
respectively. 

In addition, some measurement error is expected to affect the 
individual measurements. Similarly to what was done in R09, 
our estimate of the covariance of broadband magnitudes syn- 
thesized from STIS spectra is based on the repeatability of the 
monitoring spectra of AGK +81 266. The individually measured 
spectra were integrated in the MegaCam and SDSS passbands 
and an empirical covariance matrix of the measurements was 
built. We found a repeatability of 5, 4, 5, 6 and 11 mmag in 
broadband magnitudes ugriz for both instruments, with signifi- 
cant correlation between neighboring bands, and (as expected) 
nearly perfect correlation between overlapping bands. The full 
covariance matrix is delivered in table E.10. Some of the stan- 
dard spectra were established with several independent measure- 
ments. When determining the expected broadband noise of the 
star measurement the covariance matrix is divided by the number 
of observations. 

9.3. Mixing calibration data 

For each available measurement (or pseudo-measurement in the 
case of a color-transformed magnitude) of a photometric stan- 
dard s in a photometric band b, denoted m s b , the corresponding 
expected standard AB magnitude m s , can be calculated from its 
CALSPEC spectrum. The model for the calibration data is: 

m s b ~m s b =6 b ab + e Sib , (31) 

where e is the error term affecting the estimate of mi - mf . It is 
further decomposed as: 

dm s , , 

e s ,b = n s , b + c s , b + —r-dA eS , (32) 

dA cS 

where n s j, is the measurement error affecting m s b (including sys- 
tematic uncertainties on the measurement), c s j, is the uncertainty 

affecting the broadband CALSPEC magnitudes, and ^jrSA^ s de- 
scribe the first order impact of the filter uncertainty on m s b . 

The cross-calibration data (between SNLS and SDSS as well 
as between bands im and /2m) adds further constraints. These 
take the generic form of approximate linear relations between 
very similar filters, applicable on average to the tertiary stars: 

m s b -m s b , (33) 

25 We are not interested in the overall flux scale. Otherwise the uncer- 
tainty would be dominated by the determination of the Vega absolute 
flux, which does not play any role here. 
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with some uncertainty on a and f3. 26 As before, we translate 
these color relations into cross-calibration constraints by assum- 
ing that those relations also apply to a spectrophotometric stan- 
dard lying in the color range of the transformation. The extra 
error coming from the dispersion around the color transforma- 
tion has to be estimated from stellar libraries. We thus add to our 
model the following terms: 



Pbb> ~ »bb'C s ~m s b + fh s b , = 6l h - 6 b ^ + e w 



(34) 



where (b,b') ranges over {(u M ,u), (g M ,g), (r M ,r), (i M ,i), 
(i2m, i), (zm,Z), 0'2m, im)}- Again the error term is compound: 



c s ,b - c St b> + 



yC + ( 



■+d, 



where d is an uncertainty term accounting for the dispersion of 
stars around the color relation. We obtain the expected rms of 
the d term by measuring the dispersion around the average trans- 
formation of magnitudes synthesized from a stellar spectrum li- 
brary. We find 48, 6, 2, 4, 3 and 2 mmag of dispersion around 
the relations between the SNLS and corresponding SDSS band 
for u M ,gM,rM,iM,i2 M ,z M respectively. 

To ease the handling of correlated uncertainties between all 
the measurements, we can gather the 91 calibration measure- 
ments into a single vector y and write the linear model for the 
calibration data in matrix form: 



y = Ax + e 



(35) 



We want to solve for the vector of calibration offsets x = 
(5"k , • • ■ ,d z f ,<5". , • • • , S z ,). The Jacobian matrix A relates each 
measurement to the corresponding offset (or a difference be- 
tween two offsets in the case of cross-calibration constraints), 
and e is the vector of error terms. 



9.4. Covariance matrix of calibration data 

Most of the measurement systematic uncertainties coherently 
affect different measurements (even across different bands). 
Neglecting the correlations would significantly underestimate 
the true error. We thus formed the full covariance matrix R = 
cov(c) of the calibration measurements accounting for all the er- 
ror terms in the model above. 

We already described the computation of the covariance of 
broadband CALSPEC magnitudes in section 9.2. The computa- 
tion of the resulting contribution to R accounts for the STIS mea- 
surement error and the uncertainty on white dwarf primary cal- 
ibration differently. The STIS measurement errors are assumed 
to be uncorrelated from one star to another and divided by the 
number of observations for each star. The calibration error is as- 
sumed to be a systematic error that produces the same error for 
all the measurements of all the stars. 

The error affecting the determination of the filter effective 
wavelength is assumed unrelated from one filter to another. Its 
contribution R F to R is obtained as: 



H+diagO^H 



26 A complete interpretation of these relations is not trivial as it re- 
quires some spectral model of the stellar population composing the ter- 
tiary standard stars. The estimation of resulting uncertainties would be 
somewhat complex. We prefer to rely on a star with available and re- 
liable spectrophotometry and model how it deviates from the average 
stellar population. 



Table 22. Correlation matrix of the final calibration. 
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where H is the matrix holding the terms corresponding to 

each measurement. 

Due to the large number of systematic uncertainties to be 
accounted for, the construction of the covariance matrix of the 
measurement error terms is less straightforward. In brief, we 
have assumed that all statistical uncertainties arising from mea- 
surement noise (photon noise or flat-field noise) are uncorrelated 
from one measurement to another, whereas part of the systematic 
uncertainties affect measurements coherently and are therefore 
taken to be fully correlated. Full details of the assumptions made 
in the covariance matrix computation are given in appendix E.4. 

9.5. Least-square estimate of the calibration 

We derived the combined estimate of the calibration offsets 
through weighted least-square minimization. The resulting off- 
sets are given in table 23, along with the diagonal component 
of their covariance matrix. The corresponding correlation ma- 
trix is given in table 22. The offsets are applied to the catalog of 
tertiary standard stars to provide our best estimate of their AB 
magnitudes. 

The global x 1 of the fit is 79.8 for 80 degrees of freedom. 
Note that we are mainly manipulating actual systematic uncer- 
tainties rather than Gaussian-distributed random variables. A 
better assessment of the agreement between the different sets of 
primary standard measurements is provided by Fig. 29. We defer 
the discussion of these results to the next section of the paper. 



10. Results 

10.1. Calibration products 

We deliver natural broadband AB magnitudes for a large set of 
selected stars (the tertiary standards) in the SNLS and SDSS sci- 
ence fields. These tertiary standards will be used to calibrate the 
science images of the surveys and may be useful for other, simi- 
lar purposes. 

The calibration of the catalogs to the AB system results 
from the combination of all the calibration data gathered in this 
work. The details of the combination are described in Sect. 9. 
Briefly, we combined calibration constraints coming from obser- 
vations of 5 different spectrophotometric standard stars from the 
CALSPEC database. Those observations were conducted fol- 
lowing 3 mostly independent paths presented in Sect. 2.4 and 
described respectively in Sect. 5.1, Sect. 5.2 and Sect. 7. In ad- 
dition, the SNLS and SDSS tertiary catalogs are tied together by 
a cross-calibration measurement described in Sect. 8. 
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We emphasize that the magnitudes are delivered in the natu- 
ral photometric systems of the corresponding instrument, as de- 
fined in Sect. 2.2.1. As a consequence, their correct interpreta- 
tion is given by Eq. (3) using the effective filters corresponding 
to the average position of the star observation in the field of view. 
We deliver a model of MegaCam passbands along with this pa- 
per and refer the reader to the work from Doi et al. (2010) for 
the SDSS effective passbands. 

Another consequence of the variation of the natural pho- 
tometric system with the position in the field of view is that, 
for calibration purposes, the reference catalogs are best matched 
with exposures sharing a compatible pointing (same pointing for 
MegaCam or same declination for SDSS). In other cases, small 
color terms accounting for the slight differences in the effective 
filters may have to be considered to obtain accurate zero points. 

One may also note that the catalogs have been established 
using IQ-scaled aperture photometry in the case of SNLS and 
the usual PSF photometry algorithm in the case of SDSS. While 
the result does not depend on the photometry method to first or- 
der, small spatial or chromatic effects may need to be accounted 
for when comparing to measurements obtained using different 
photometry methods (see e.g. Sect. 4.3.3). 

For most comparison purposes, it is also convenient to re- 
port MegaCam measurements in a single photometric system, 
close to the natural system. In this paper, we follow the conven- 
tion adopted in R09 and define the MegaCam uniform magni- 
tudes as the natural magnitudes that would have been observed 
for objects at the center of the focal plane 27 . We have derived 
color transformations between the natural and uniform system 
(see Sect. 3.2) and include in the catalog MegaCam uniform 
magnitudes for all stars that can be accurately transformed. 

Finally, we determine measured color terms between the uni- 
form MegaCam photometric system and the average SDSS sys- 
tem in Sect. 8.3. These facilitate comparison of other data specif- 
ically obtained in these photometric systems. 

Catalogs can be found in tables E. 1 and E.2. 28 AB offsets 
given in table 23 are already applied to the released catalogs. 
The covariance matrix associated with this calibration is given 
in tables 23 and 22. Finally the Megacam effective transmission 
models to be used in Eq (2) are given in table E.3-E.8. 29 What 
follows is a general discussion of these results. 

10.2. Agreement between the available calibration paths 

A key feature of this work is the redundancy between several 
calibration paths established independently and subject to dif- 
ferent systematics. An illustration of the agreement between the 
different paths is provided by Fig. 29. 

The agreement between the SNLS and SDSS calibration ap- 
pears better than one percent in the gri and z bands. There is 
an apparent offset of about 4% in u which is still compatible 
with the low precision of the cross-calibration data in this band. 
Internally, there is a small tension between the Landolt and so- 
lar analog calibrations of the SNLS in /mA'2m an d r M (~ 1-9 
and 1.3%). To a lesser extent, a similar trend is visible in the 
SDSS which suggests that at least a part of the tension may 
be explained by an actual calibration difference between the 

27 Note that this choice is not representative of the "average" 
MegaCam natural system. Passbands at the center are redder than the 
spatial average. 

28 Full catalogs can be retrieved from the CDS or the dedicated web- 
page http : //supernovae . in2p3 . fr/snls_sdss/. 

29 Higher resolution versions can be found on the webpage. 
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Fig. 29. Agreement between complementary determinations of 
the SNLS AB offsets. SNLS BD+17 and SNLS SA lines refer to 
the SNLS calibration on Landolt standards and direct solar ana- 
log measurements respectively (cf. Sect. 4). The SDSS SA and 
SDSS BD+17 lines refer to the calibration obtained by translat- 
ing the PT measurements of the solar analogs and BD +17 4708 
respectively, to the SNLS system via the relations obtained in 
Sect. 8. The data points correspond to the calibration offset rela- 
tive to the combined calibration including all the measurements 
and are not independent. The shaded region illustrates the uncer- 
tainty on the combined calibration. 



BD +17 4708 and the solar analog CALSPEC reference spec- 
tra. 

Overall, we reckon that the whole data sample is compatible 
with the modeled uncertainties. Consequently we opted for com- 
bining the whole data sample to set the flux scale of the tertiary 
catalogs and we believe that the uncertainties resulting from the 
combination are reliable. 



10.3. Uncertainties on the combined calibration 

A few points pertaining to the calibration uncertainties ought to 
be discussed. 

First, one can observe from the correlation matrix in table 22 
that the cross-calibration data brings significant correlations in 
the calibration of all the SNLS and the SDSS overlapping bands 
but u and g, where the accuracy of the cross-calibration is lim- 
ited by the less-similar filters and the higher dispersion of the 
stellar population about the color- tranformations. Consequently, 
flux measurements in the two surveys can now be compared with 
low cross-calibration uncertainty. 

It is also interesting to compare the relative importance of 
the various contributions to the global calibration uncertainty. 
The last row of table 23 shows that the uncertainty on the pri- 
mary standard spectra amounts for about half the uncertainty (in 
variance) for the redder bands (riz). The imperfect knowledge of 
filters impacts mostly the u and um band calibration due to the 
large variation of the standard star spectra in this wavelength re- 
gion and the tm band due to the larger uncertainty remaining on 
the transmission of this filter. 

Evidence for potential evolution of the tm and im filters has 
been found and the resulting uncertainty on the filter passbands 
is discussed in Sect. 9. 1 . We emphasize that the impact of filter 
passband uncertainty depends on the color of the object whose 
flux is being determined. Specifically, this error cancels for ob- 
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Table 23. SDSS and SNLS combined AB offsets. 
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With perfect passbands" 


0.0122 
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0.0030 


0.0044 


0.0086 


0.0047 


0.0022 


0.0024 
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Notes. Figures are given in magnitudes. 

(fl) Remaining uncertainty on the AB calibration when assuming that the passband measurements are perfect. 
(A) Remaining uncertainty on the AB calibration when assuming that the STIS spectra are perfect. 



Table 24. Calibration offset between R09 and the present release 
of the SNLS calibration. 
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offset' 1 
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0.004 
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Notes. 

(<I) The offset should be added to the previous release to match the 
present result, and comes in addition to a sign error correction in the 
previous release (see text). 

{h) Rms of the differences between the two catalogs for stars brighter 
than r M < 18. 



jects whose spectrum is close to the spectrum of our standard 
stars (main sequence stars of color g — i~ 0.55). This is why the 
poor knowledge of MegaCam i M and r M passbands only min- 
imally affects AB calibration accuracy: the primary standards 
we observe have SEDs close to the AB spectrum in those bands 
(stars with g—i~ 0.55 have r—i~ 0.1). This is not the case for 
type la supernovae; therefore, the impact of filter passband un- 
certainty on their flux measurements is bound to be much larger. 
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Fig. 30. Difference between SNLS-calibrated tertiary star mag- 
nitudes and previous R09 release as a function of the distance 
to the focal plane center. The 4 deep fields are compared sepa- 
rately (Dl: red plus, D2: green squares, D3: black circles, D4: 
blue crosses.) 



10.4. Comparison with SNLS3 

For the four photometric bands already calibrated in R09, we 
compared uniformity and absolute calibrations between the new 
and prior tertiary catalogs. A sign mistake affecting the R09 cat- 
alog was discovered in this comparison. It contributes to a slight 
deterioration of the catalog uniformity and shifts the average flux 
scale of the catalog by 8, 5, 6, 3 mmag on average for band 
gMi r M^M an d Zm respectively. Once the error sign is corrected, 
the remaining offsets are given in table 24. The comparison of 
calibrated magnitude differences as a function of the distance to 
the center of the focal plane is given in Fig. 30. 

10.5. Comparison with the previous calibration of the SDSS 
SN survey. 

Following the discovery of a non-uniformity in the SDSS PT 
flat-field, this work introduces two changes to the SDSS ter- 
tiary catalogs. The first is a correction of the tertiary cata- 
log flat-fielding based on the Data Release 8 SDSS catalog 
(Padmanabhan et al. 2008) flat-fielding solution. The second is 
the new determination of the corresponding AB offsets. The pre- 
vious calibration of the SDSS SN survey was based on the un- 
corrected version of the coadd catalog scaled to the AB system 
by the application of AB offsets from Holtzman et al. (2008, ta- 
ble 1). As supernovae are uniformly distributed in declination, 
the calibration offset to be expected from this release is readily 
given by the weighted average of magnitude differences between 
the two calibrated catalogs. This offset is given in table 25. 



Table 25. Calibration offsets between first and third year release 
of the SDSS-II SN survey. 
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0.007 
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0.035 


0.011 
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-0.031 


-0.004 
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-0.000 


-0.006 



Notes. 

(a) Calibration offset in mag (average difference between the present 
work and Holtzman et al. (2008)). 

<w Standard deviation of the photometry difference between the two 
catalogs (in mag). 

w Peak to peak difference between the two catalogs (in mag). 



The difference in the average photometry is fairly small com- 
pared to the size of the PT flat-field correction. This is a fortunate 
consequence of the fact that the effect is close to symmetric with 
respect to the center of the PT focal plane, where standard stars 
were observed. The row labeled rms in table 25 gives an estimate 
of the gain in photometric precision that is to be expected from 
the improved flat-fielding. 

10.6. Comparison with other work 

While we have achieved excellent agreement with the SNLS and 
SDSS supernova surveys we did not attempt to compare our 
photometry with any other survey. The SDSS SN photometry 
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has been compared with CSP (Mosher et al. 2012), where good 
agreement was found. The Pan-STARRS survey (Hodapp et al. 
2004) has recently reported a small discrepancy with the SDSS 
stellar photometry (Tonry et al. 2012), but larger than might be 
expected on the basis of this work. Their analysis used the pub- 
licly available SDSS catalogs which lack the corrections applied 
in this paper and they did not use the filter transmission mea- 
surements of Doi, et al. However, the discrepancy in photometry 
might very well persist even with our revised catalog. 

11. Conclusion 

We have delivered updates to the tertiary standard star cata- 
logs that define the photometric calibration of the SNLS and the 
SDSS supernovae surveys. The catalogs provide standard star 
calibrated AB magnitudes, in the natural photometric system of 
both instruments. The products are fully described in Sect. 10.1. 

The flux calibration is based on the observation of 5 dif- 
ferent HST primary spectrophotometric standard stars from the 
CALSPEC database. In addition, the SNLS and SDSS tertiary 
standard stars are now tied together by specifically designed 
cross-calibration observations. Three independent routes, each 
subject to different observational systematics, were used to trans- 
fer the HST primary standard star calibration to the tertiary stan- 
dard stars. Good agreement was found between these 3 routes, 
which were then combined to provide the final products. In par- 
ticular, the offsets derived in Sect. 9 from simultaneously fitting 
all available information have been applied to the delivered cat- 
alogs. These offsets have been determined with CALSPEC ver- 
sion 3 spectra, and may be updated as new CALSPEC versions 
are released. See Appendix E for more details. 

With respect to previous releases, the photometric calibration 
evolved by less than 1% in all photometric bands but u. Major 
improvements arise from better flat-fielding processes (that also 
improve the uniformity of the resulting catalogs), and the intro- 
duction of more data to constrain the calibration transfer. 

The calibration transfer from HST standard stars to stars in 
supernovae fields is thought to be accurate to about 3 mmag in 
the riz photometric bands. This is comparable, or better than the 
error related to the calibration uncertainty of the primary stan- 
dard stars themselves. The error in band g is about 4 mmag and 
remains dominated by the error in the calibration transfer. The 
uncertainty on the u band calibration amounts to about 1% in 
the SDSS and 2% in the SNLS, dominated by the uncertainty on 
the instrument effective passbands. 

The photometric calibration was complicated by the need to 
overcome several instrumental effects as well as weaknesses in 
instrument or survey designs. The most important complications 
are summarized here: 

- Uniform illumination of wide-field instruments is not 
enough to achieve accurate flat-fielding; internal reflections 
in the optical path must be accounted for. At least two ele- 
ments, CCDs and filters, are reflective enough to feed sig- 
nificant amounts of light back into the wide-field corrector 
optics. Reflected by convex lens surfaces, this light causes 
ghosting as well as diffuse and structured contamination in 
the flat-fields. For this reason, photometric uniformity of 
both the SDSS and SNLS instruments ultimately relies on 
the analysis of dithered stellar observations. The variation 
of atmospheric conditions during the observation sequences 
ultimately limits the flat-field solution accuracy. 

- Similarly, small-angle scattering of light caused by dust 
or defects in the optics can affect the photometry of point 



sources while leaving the flat-field twilight images un- 
changed. A remarkable illustration of such a phenomenon 
occurred during the SNLS survey: the progressive deposit of 
metal dust from the filter jukebox on one corner of the top 
corrector lens caused a local 6% deterioration of the effec- 
tive throughput of the photometry, and was only detected and 
measured through dithered stellar observations. Such effects 
can also be detected by monitoring apparent star brightness 
in often-visited fields. 

- For SNLS, spatial variations of the instrument effective pass- 
bands had to be accounted for in order to accurately conduct 
the calibration transfer and deliver relevant flat-fielding solu- 
tions. An extensive survey of the filters' transmission would 
make the calibration more accurate. 

- The characterization of the SNLS photometric system would 
have been simplified by the implementation of a well- 
planned filter passband monitoring program. Nonetheless, 
no significant chromatic differences have been observed in 
repeat stellar photometry over the course of the survey, pro- 
viding direct evidence of the photometric system's long-term 
stability. 

- Slight readout electronics instabilities cause percent-level 
noise in the MegaCam flat-fielding. For the SNLS, this effect 
is mitigated by repeated measurements over the long dura- 
tion of the survey; for shorter surveys, this effect would have 
been difficult to overcome. Reproducible illumination pro- 
vided by artificial light sources might provide more reliable 
monitoring of the instrument throughput. 

- The SDSS calibration relied on a separate telescope (the PT), 
but neither the flat-field nor the filter throughput was well- 
characterized enough for the precision photometry goals 
of this paper. Instead we had to rely on empirical color- 
transformations and correct the flat-field from external data. 

Most of these effects were overcome by introducing redun- 
dancy in the measurements, or by the analysis of specifically 
designed additional observations. While such strategies allowed 
us to reach an accuracy of 0.3%, it would be prohibitively time- 
consuming to reduce this accuracy further. Future photometric 
surveys such as SkyMapper (Keller et al. 2007), DES (Bernstein 
et al. 2012) and LSST (LSST Science Collaborations et al. 
2009) might take advantage of current developments regard- 
ing instrumental calibration with artificial calibration sources. 
Information provided by the SNDICE illumination prototype 
(Barrelet & Juramy 2008; Guyonnet et al. in prep.) installed 
at the CFHT has already proved immensely valuable in under- 
standing MegaPrime instrumental effects such as flat-field pollu- 
tion with ghosts and readout electronics gain variations. Ongoing 
analysis of SNDICE data may also give us a final understanding 
of the MegaCam filter passbands. A comparable device (Stubbs 
et al. 2010) has aided in the calibration of the Pan-STARRS in- 
strument; the Dark Energy Survey is taking a similar approach 
with DECal (Rheault et al. 2012). 
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Appendix A: Monitoring of the MegaCam flat-field 
variability from deep fields 

Variations of the instrumental response are determined by the 
time distribution of the photometric coefficient 5Z tccd (where t 
indexes the exposures, and ccd indexes the 36 MegaCam chips. 



We recall that 6Z tccd is computed for each deep field by build- 
ing an average photometric catalog of non-variable stars as de- 
scribed in Sect. 4. We then compute for each available exposure 
a zero point relative to the average catalog. Once the exposure 
t is scaled by this zero point, we recompute a zero point offset 
for each CCD (5Z uccd ) relative to the average catalog. Non-zero 
offsets indicate that the throughput of the instrument has varied 
non uniformly. This can be due either to clouds or to a variation 
of the instrument response that is not corrected by the twilight 
flat-fielding. 

For each exposure t we compute the rms of SZ, ccd of the 36 
CCDs. The statistical uncertainty on the determination of 5Z tccd 
is typically 2 mmag, so that sets the minimum expected value. 
The rms values observed throughout the survey are displayed in 
Fig. A. 1 . The measured rms is typically about 5 mmag. The ob- 
served excess relative to the minimum value of 2 mmag is to be 
expected for at least 3 reasons: 1) structured variations of the at- 
mospheric transparency, such as clouds, affect the observations, 
2) the gains of some CCDs amplifiers vary on time scales shorter 
than a single run, and therefore can't be calibrated by a single 
twilight flat-field, and 3) telescope modifications can have an 
impact on the photometric response, while leaving the twilight 
observations virtually unchanged. The effects of 1) and 3) can be 
observed directly in Fig. A. 1, while 2) is a noise that contributes 
to a global increase of the measured rms. 

Clouds are occasionally responsible for a large increase in 
the measured rms. These non-photometric conditions are easily 
detected and excluded from the averaging of the tertiary cata- 
logs. Exposures detected as non photometric are displayed as 
black dots on Fig. A.l. 

We have displayed telescope modifications that had a notice- 
able impact on the uniformity of the photometry as red vertical 
lines. The beginning of the survey was marked by continuous 
adjustment of the optical setup. We retained only the main mod- 
ification (that had the greatest impact on the image quality): the 
flip of the third lens in the wide field corrector. It is denoted as 
event A on the graph. The rms of exposures obtained prior to 
this event can be seen to be higher. The second and most sig- 
nificant phenomenon was the accumulation of metal dust on a 
corner of the wide field corrector optics. Its effect was the scat- 
tering at small angles of a growing part of the light, decreasing 
progressively the flux fraction contained in the core of the PSF 
(and hence the throughput of the photometry), while leaving the 
twilight images essentially unchanged. The effect is clearly vis- 
ible in Fig. A.l as a progressive increase of the rms during the 
year 2006. Its detection and cleaning occurred at the end of 2006 
and is denoted as event C. The beginning of the accumulation is 
harder to determine. We approximately determined that the ef- 
fect becomes sufficiently important at the middle of 2006 (event 
B). The last events are the destruction (D) and replacement (E) 
of the im filter in 2007. Exposures in gM and ru were apparently 
also affected by the same event (D), but we do not understand 
why they changed. 

The events are summarized in table A.l. They define 5 pe- 
riods for which independent photometric corrections can be de- 
termined. The application of our set of photometric corrections 
to the survey flat-fielding is illustrated in Fig. A. 2. This figure is 
similar to Fig. A.l, but made after the application of photomet- 
ric corrections. After correction, the measurements are consis- 
tent within an error cr/f, which is estimated as the average rms 
of 6Z ccd . Values of <Tff were given in table 2 in the main text. 
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Table A.l. Grid observations and associated telescope events. Table A.2. Dithering pattern of the grid exposures. 



Run" 


Date 


Bands'" 


Field 


Events 


03Bm02 


2003-10-02 


i 


2 






2003-10-03 


r 


2 




03Bm03 


2003-10-22 


g,u,z 


2 






2003-10-23 


i,r 


2 




04Am02 


2004-03-18 


u,g,r,i,z 


1 




04Bm05 


2004-12-02 
2004-12-02 


g,r,i 


1 


(A) L3 lens nipped 
upside down. 




2004-12-10 


g,u,z 


1 




05Bm02 


2005-10-09 


u,g,r,i,z 


2 




06Am04 


2006-05-19 






(B) Growing metal 
dust deposit. 


06Bm02 


2006-09-13 


i 


2 






2006-09-14 


r 


2 






2006-09-15 


z 


2 






2006-09-16 


u 


2 






2006-09-19 


g 


2 




06Bm05 


2006-12-27 






(C) Top lens clean- 
ing. 


07Am02 


2007-03-21 


r 


1 






2007-03-22 


gu 


1 






2007-03-23 


' 


1 






O f\f\'~l no 1 A 

2007-03-24 


z 


1 




07Am06 


2007-07-20 






(D) i M filter broken. 


07Bm01 


2007-08-24 






(E) Recoating. 


07Bm03 


2007-10-18 






(F) New i2 M filter 
installed. 


07Bm04 


2007-11-08 


i2 


1 






2007-11-11 


r 


1 






2007-11-12 


g 


1 




07Bm05 


2007-12-03 


u,z 


1 




08Bm02 


2008-09-03 


u,g,r,i2,z 


2 




09Am06 


2009-07-20 


g,z 


2 






2009-07-21 


r,u,i2 


2 




lOAmOl 


2010-02-06 


u 


1 






2010-02-18 


r.z 


1 






2010-02-19 


i2 


1 






2010-02-20 


g 


1 




10Bm02 


2010-10-05 


r 


1 






2010-10-09 


g,i2 


2 






2010-10-13 


r 


2 




10Bm03 


2010-11-08 


g.u 


2 




HAm06 


2011-07-30 


z 


2 






2011-08-02 


g,r,u,y 


2 




HBmOl 


2011-08-28 


g,r,z 


2 






2011-08-29 


g.u.y 


2 





Notes. (o) Run id in the elixir pipeline nomenclature. 

{b> Observations selected as valid are bold-faced in this column. 



Exposure 


a (deg.) fl 


8 (deg.) fl 


1 


0.000 


0.000 


2 


0.028 


0.000 


3 


0.056 


0.000 


4 


0.112 


0.000 


5 


0.225 


0.000 


6 


0.451 


0.000 


7 


0.902 


0.000 


8 


0.000 


0.027 


9 


0.000 


0.055 


10 


0.000 


0.111 


11 


0.000 


0.222 


12 


0.000 


0.444 


13 


0.000 


0.888 


14* 


0.000 


0.000 



Notes. (fl) Pointing deviation relative to the reference pointing of the 
field. ib) A concluding 14th exposure, that comes back to the first point- 
ing has been added to grid sequences since September 2008. This re- 
dundancy provides a way to monitor the the stability of the observation 
conditions during the sequence. 
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Fig. B.l. Comparison of available measurements of the REOSC 
filter transmission curves at a distance of 15cm from the center. 



Appendix B: MegaCam filters measurements 

The transmission curves of the original Sagem/REOSC 
MegaCam filters were measured twice: the first time by the man- 
ufacturer before their installation in MegaCam in 2002, and the 
second time in 2006 by the CFHT operational team. The 2006 
measurements are available only at positions close to the bor- 
der of the filter due to limitations in the bench clearance, and 
their spectral resolution is coarser. Comparison between the two 
measurements at 15cm from the center (after convolution of the 
REOSC measurements to match resolutions) is given in Fig. B.l. 
The position of the and filter cut-offs can be seen to be 
shifted ~ 8 nm redder in the 2006 measurements. The origin 



of these discrepancies and the reason why only two filters are 
affected are unclear. 

A confirmation of the 2006 measurements was recently 
given by the measurement of the transmission of one of the 
broken pieces of the i M filter. The measurement was conducted 
on the calibration bench of the SNDICE illumination device 
(Barrelet & Juramy 2008), and display a perfect agreement with 
the 2006 transmission curves 12 cm away from the center (the 
appropriate curve for the broken piece). 

The possibility of a continuous (slow) aging of the filters is 
excluded by the examination of color-terms between the begin- 
ning and the end of the survey (cf. Sect. 3.4.1). It is thus likely 
that the 2006 measurements are the most representative of the 
state of MegaCam filters during the survey. The lack for mea- 
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Table C.l. The PT data with secondary patch overlaps in SDSS 
stripe 82. 



Table C.l. The results of PT fits. 



Field 


Date (MJD) RA (deg) 


Dec (deg) 


Frame 


82A12 


52551 


315.0 


-1.00 


161156 


82A16 


52551 


25.0 


-1.00 


161196 


82B12 


52577 


315.0 


-0.33 


165233 


82B16 


52551 


25.0 


-0.33 


161246 


82C16 


52577 


25.0 


0.33 


165263 


82D12 


52577 


315.0 


1.00 


162033 


82D16 


52584 


25.0 


1.00 


166668 


82E12 


52557 


315.0 


-0.66 


162058 


2E16 


52553 


25.0 


-0.66 


161813 


82F12 


52557 


315.0 


0.00 


162068 


82F16 


52553 


25.0 


0.00 


161823 


82G12 


52578 


315.0 


0.66 


165524 


82G16 


52557 


25.0 


0.66 


162098 



surements at the center of the filters in the 2006 data requires 
additional data or assumptions. We opt for correcting the iu and 
ru REOSC measurements to make them match the 2006 trans- 
mission curves. The application of a linear stretch in wavelength 
of the form: 

T*(A) = T^A 1 + ^^-(A-A 1 ) 

was found to reconcile both measurement sets. We 
used (AuA 2 ,A 3 ) = = (5300A,6837A,6891A) for r M and 
(6600A, 8464A, 8539A) for i M . The resulting transmission 
model at 15cm is shown as a dashed black line in Fig. B.l; 
similar agreement is obtained with the 12cm data. 

Appendix C: SDSS Photometric Telescope (PT) 
flat-field 

There are some overlapping secondary patches in stripe 82 that 
can be used to check for deviations from a uniform response of 
the PT. This section provides an analysis of that data. 

The data sample consists of the 13 secondary patches listed 
in Table C.l. The fields labeled A, B, C, and D are the standard 
locations in declination for secondary patches in stripe 82. The 
field of view is 40 arc-minutes (2/3 degree), so there is no over- 
lap between the A, B, C, and D secondary patches. The fields 
labeled E, F, and G are special patches that were placed approx- 
imately half-way in declination between the standard locations, 
providing a 50% overlap with the standard patches. The camera 
pointings are all in the same relative orientation with East- West 
being in the CCD column direction. 

The analysis involved matching stars in the overlapping 
pointings and computing the difference in magnitude. The data 
come from the SDSS data processing (mtPipe) (Tucker et al. 
2006). The differences in magnitude measure the difference in 
response of the focal plane between a given column and another 
one that is displaced by 1/2 the size of the CCD. If the response 
is uniform all the differences would be zero. It is convenient to 
measure the distance along the focal plane in terms of degrees 
of declination. We determine this distance as the declination of 
a matched star minus the reported declination of the telescope 
pointing. The matching star is nominally at a declination lower 
by 20 arc-minutes. 

The results are summarized in Table C.2. The table gives the 
total number of comparison stars that were found and a fit to the 
equation 

Am = a6 + b (C.l) 



Band 


Stars 


x 1 


a 




b 


0> 


;( 


699 


5.6 


-0.129 


0.022 


0.041 


0.004 


8 


2514 


2.2 


-0.039 


0.009 


0.027 


0.002 


r 


2971 


6.5 


-0.037 


0.008 


0.022 


0.001 


i 


2581 


8.0 


-0.019 


0.009 


0.023 


0.002 


z 


1127 


3.0 


-0.053 


0.016 


0.038 


0.003 


gri 


8066 


39.5 


-0.032 


0.005 


0.024 


0.001 



where Am is the magnitude difference, 5 is the position on the 
focal plane in units of degrees of declination, and a and b are fit 
parameters. Ax 2 is given for each fit. The fit is made to 6 binned 
data points so there are 4 degrees of freedom for each fit, except 
for the combined fit to g, r, and i (last line in Table C.2). This fit 
is performed on 18 binned data points from g, r, and i so there 
are 16 degrees of freedom for that fit. The chi-squared's are gen- 
erally good although they could probably be improved by the ad- 
dition of a quadratic term. However, it was felt that a quadratic 
term was unnecessary and possibly sensitive to unknown sys- 
tematic errors. Every filter shows a consistent trend of positive 
residuals and a negative slope with declination. However, the fits 
are quite precise and there are statistically significant differences 
between the filters. 

The results above only give the difference in response be- 
tween the two halves of the PT. Determining the actual response 
requires additional data or assumptions. There seem to be two 
obvious models that could explain some features of the data: 

- The CCD readout uses different amplifiers for the two halves 
of the chip we are considering. It is possible that the gain 
difference between the two amplifiers has not been properly 
calibrated. Such a mis-calibration could occur, for example, 
if the light levels in the flat-field exposures resulted in a non- 
linear amplifier response that was different for the two am- 
plifiers while normal exposure levels remained in the linear 
regime. This scenario could explain the offset between the 
two halves of the CCD but doesn't seem very plausible and 
would not explain the observed slopes. 

- The illumination of the flat-field exposures may be non- 
uniform. This would lead to an error in the calibration of 
the CCD response. The non-uniformity is most likely to be a 
slowly varying function of position. 

In any case the data are not compatible with a uniform photomet- 
ric response in the PT flat-fielded exposures. The PT flat-field 
error has two consequences. First it is transferred to the 2.5 m as 
a flat-fielding error. It should imprint a periodic pattern 0.6 de- 
grees wide in declination on the photometry. Second, as primary 
standard stars are observed preferentially at the center of the fo- 
cal plane, it can bias the calibration relative to the tertiary stars 
that are observed throughout the focal plane. It is thus necessary 
to provide a correction of the 2.5 m flat-field independent of the 
PT and to characterize the PT flat-field error to derive accurate 
calibration from its measurement. 



Appendix D: DR8 catalog retrieval 

The Stripe 82 DR8 catalog has been retrieved from the DR8 
database using the following query: 

select psfmag_u ,psfmag_g ,psfmag_r ,psfmag_i , 
psfmag_z , ra , dec , run , camcol , psfmagerr_u , 
psfmagerr_g ,psfmagerr_r , psfmagerr_i , 
psfmagerr_z into mydb . MyTable from 
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phototag where dec between -1.3 
and 1.3 and ((ra between 300 and 360) 
or (ra between and 60)) and type = 6 
and mode = 1 and psfmag r < 22.5 



Appendix E: SNLS and SDSS calibration products 

Our main delivery are catalogs of natural magnitudes for a large 
set of tertiary standard stars calibrated to the AB flux scale. We 
also provide the material required for the interpretation of those 
magnitudes in term of physical fluxes (instrument transmission 
functions) as well as the material required to update the AB cal- 
ibration as new CALSPEC spectra for our primary standards are 
released. 

All the material described below is available from the SNLS 
and SDSS joint calibration webpage: http://supernovae. 
in2p3 . fr/snls_sdss/. 

E.1. Tertiary stars catalogs 

We provide catalogs of natural AB magnitudes in each of the 
MegaCam filters for selected tertiary standard stars in the four 
CFHTLS deep fields in the natural MegaCam photometric sys- 
tem. This system is described in Sect. 3.1 and the corresponding 
transmission curves are given below (Sect. E.2). Note that the 
transmission depends on the location on the focal plane. A sub- 
sample of the Dl catalog is given in table E.2. 

Similarly, we deliver natural SDSS AB magnitudes for se- 
lected tertiaries in the Stripe 82. The magnitudes are given for 
the average filter response as given by Doi et al. (2010), i.e., 
they should not be adjusted for the small differences in the filter 
responses of the different CCD columns. An excerpt from the 
catalog is given in Table E. 1 . 

In both cases the released magnitudes are already placed on 
the AB system by the addition of the zero points given in table 
14. Full catalogs in electronic format can be retrieved from the 
CDS or from the SNLS and SDSS joint calibration webpage. 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 























1 


















+ 
















". 
















+ 


f++ + 


f + + 




....+ 






....+. 


t- 

+ + 




- + ++ + ++ 


++ 




h+ + 

1 


+ 


++ 

H 1 

i 






- 




i '1 
1 1 

+ 


++ 










- + 
































1 










1 




























































































f- — 




— 


















+ ^ 


■+--•- . 








_,_+ 






-t— — " 








































































































































































-+- 
— - ~ t ~—*- 






















-F 



















































0.02 
0.01 
0.00 
-0.01 
-0.02 

0.02 
0.01 
0.00 
-0.01 
-0.02 

0.02 
0.01 
0.00 
-0.01 
-0.02 

0.02 
0.01 
0.00 
-0.01 
-0.02 

0.02 
0.01 
0.00 
-0.01 
-0.02 

0.02 
0.01 
0.00 
-0.01 
-0.02 









































































4 


■+■ , 


. ■*■ 








- 


-+- 










. . . — 


+ . 














H 


-+ 











































































































































































-+- 

_ _~ . 










-t- 








, -*- 

— 




+ 











































































































































.-+■ 














































- 












-t- 


+ 








. . . . ! -t 














-1- 






■+- 

























































































-0.6 



-0.3 0.0 0.3 

6 (Deg.) 



0.6 



-0.6 



-0.3 0.0 0.3 

a (Deg.) 



0.6 



9 

(cm) 



12 



Fig. 28. Relative uniformity of SNLS and SDSS photometry. 
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Fig. A.l. The standard deviation of the relative CCD zero points for each exposure in the survey. Black dots denote exposures that 
were excluded from the average because of cuts on the stability of the observing conditions or on the uniformity of the effective 
instrument response. Exposures are ordered by increasing observation date. The vertical gray dotted lines indicate the separation 
between different MegaCam runs (roughly one run per lunation). The twilight flat-field is common to all observations in a run. The 
vertical red lines indicate significant changes to the telescope. The letters refer to the events described in table A. 1 . In this figure, 
the same photometric flat-field correction, determined from observations on 2005-10-09 (black vertical line), is applied to all the 
observations. 
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04A 
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Fig. A.2. This figure is similar to Fig. A.l, but different photometric flat-field corrections were used for the different periods of the 
survey. The vertical black lines indicate the date of grid observations used to compute the photometric corrections. The changes of 
data point color indicate the different flat-field corrections applied. These corrections were used to produce the multi-epoch tertiary 
catalogs described in Sect. 4. 
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M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 



S C3 



5 =3 



o 
o 


s 

c 


O 

o 


s 

C 


O 
O 


CO 
O 


o 
o 


CN 

c 
c 


CN 
O 


cn 
c 
c 


CN 
O 

o 


cn 
c 
c 


r- 
o 


s 

C 


O 

o 


o 

C 


Tt 
o 


SO 
CN 
C 


o 
o 


cn 
o 


o 
o 


cn 
c 
c 


CO 
O 


P- 

CN 

c 


cn 
o 
o 


c 
c 


so 

o 


CO 

c 

C 


CN 

o 
o 


c 
c 


o 


OS 
C 
C 


CO 

o 

O 


cn 
c 
c 


cn 

CN 

o 


o 
c 


CN 

o 
o 


c 
c 


o 


CN 

c 
c 


o 


CN 

c 


O 


o 


CN 
O 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


+i 


+1 


+i 


+1 


+l 


+1 


+i 


+1 


+l 


+1 


+l 


+1 


+l 


+1 


+l 


+1 


+1 


+1 


+i 


+1 


+i 


■H 


+i 


+l 


+l 


+l 


+1 


+1 


+1 


+l 


+l 


+1 


+l 


+l 


+l 


+l 


■H 


■H 


+l 


+1 


+l 


+l 


+1 


+l 


+1 


Os 

as 
m 


cn 

IT; 
OS 


SO 
cn 
CN 


Os 
SO 


o 
Tt 


cn 

CO 


Tt 

CN 

m 


CO 

Tt 


o 

CN 
Tt 


IT, 

c 


SO 

m 


so 
CN 
C 


SO 
CO 


CN 
Os 


o 
o 


CN 

cn 


tt 

SO 

r- 


o 
■* 


o 

SO 
t~ 


CN 
C 

r- 


t — 

Os 
Os 


s 

CN 


cn 
q 


V") 

in 

Tl 


cn 

CN 


o 

Os 

m 


Os 
Os 


SO 
CO 

c 


SO 


cn 
SO 


in 

r~ 


SO 
Os 
CN 


Os 
so 

cn 


o 

s 


in 

CO 

tt; 


CNl 


5 


CN 

r» 
m 


Os 

o 

cn 


Os 


CN 
CO 
CO 


Os 
SO 
CN 


in 

CN 

m 


CN 

m, 
it 


c 

CO 


so 




SO 


so 




d 


so 


CO 


d 

CN 


Os 


CO 


Os 


d 

cn 


CN 


CO 




d 

cn 


CN 




d 

CN 




Os 


CN 


CN 


o< 


SO 


d 

CN 


d 

CN 


CO 


r~ 


so 


d 

CN 


d 

CN 


Os 


CN 


SO 


CO 


d 

CN 


Os 


CO 


SO 


CN 


SO 




d 

CN 


o 
o 


o 
c 


CN 

o 
o 


s 

c 


o 
o 


O 

c 
c 


o 
o 


o 
c 


s 

O 


CN 

c 
c 


o 

O 


CN 

c 
c 


SO 

o 
o 


C 
O 


o 
o 


o 
c 


m 
o 
o 


CN 
O 


o 
o 


SO 
C 
C 


o 
o 


CN 

o 
c 


CO 

o 

O 


O 


CN 

o 
o 


o 
c 


P~ 

o 
o 


s 

o 


CN 

o 
o 


o 
c 


CN 

o 
o 


SO 

c 
c 


cn 
o 
o 


CN 

c 
c 


CO 

o 

O 


o 
c 


o 

O 


■* 

o 
c 


CN 

o 
o 


o 

C 


o 
o 


c 
c 


o 

O 


o 
o 


m 
c 
c 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


+l 


+1 


+1 


+1 


+l 


+1 


+l 


+1 


+1 


+1 


+l 


+1 


+i 


+l 


+l 


+1 


+l 


+l 


+l 


+1 


+l 


+l 


+l 


+l 


+1 


+1 


■H 


-H 


+l 


+l 


+1 


+l 


+1 


+1 


+1 


+i 


+l 


+l 


+1 


-H 


+1 


+1 


+1 


+l 


+1 


Tt 

s 


f- 
o 
o 


SO 

in 
Tt 


cn 

SO 

r- 


cn 


CN 
SO 


o 

SO 
Os 


cn 
t~ 
m 


CO 
m 


IT; 

m 

cn 


o 

OS 
CO 


CO 
CO 

cn 


r- 
cn 
r- 


c 

r~ 


o 
cn 
O 


o 

SO 


o 

CN 


CN 

r- 
■n 


cn 
Tt 
co 


cn 
Tt 

CO 


Tt 

Os 


Tt 
r- 


Tt 
Tt 


CN 

cn 
so 


o 


o 
cn 
so 


Os 
CN 

q 


Os 
Os 


CO 
so 
CN 


CN 

cn 
r» 


SO 

cn 
cn 


CO 

■* 


so 

cn 


r» 

Os 
cn 


m 

Os 
CN 


Os 
m 

CO 


SO 
Os 

o 


o 

SO 
SO 


■* 

cn 


so 


m 

cn 

Os 


c 

cn 

CN 


o 

so 


CO 


CO 






SO 


sd 


CO 


CN 


SO 


CO 


d 

CN 


OS 


CO 


OS 


d 

CN 


CN 


CO 




d 

CN 


CN 




d 

CN 


CO 


Os 


CN 


CN 


Os 


SO 


CN 


O 

CN 


OS 


r~ 


so 


d 

CN 


d 

CN 


Os 


CN 


SO 


OS 


d 

CN 


Os 


Os 


sd 


CN 


sd 


CO 


d 

CN 


o 
o 


o 
c 


o 
o 


o 
c 


o 

O 


Tt 
C 

c 


o 
o 


o 
c 


cn 
o 
o 


o 
c 


O 
O 


o 
c 


cn 
o 
o 


VI 

c 
c 


O 
o 
o 


c 
c 
c 


cn 
o 
o 


r- 
c 
c 


o 
o 
o 


cn 
c 
c 


o 
o 


o 
c 


m 
O 
O 


SC 

c 
c 


o 
o 


o 
c 


o 
o 


CN 

c 
c 


o 
o 


c 
c 
c 


O 
o 


cn 
c 
c 


CN 

o 
o 


o 
c 


SO 

o 
o 


o 
c 


O 

o 


cn 
c 
c 


o 
o 


o 
c 


o 
o 
o 


■* 

c 
c 


O 
O 


o 
o 

O 


cn 
o 
c 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


+l 


+1 


+l 


+1 


+l 


+1 


+i 


+1 


-H 


+l 


+l 


+1 


+l 


+1 


+l 


+1 


+i 


+l 


+l 


+1 


+i 


-H 


+l 


+1 


+l 


■H 


+l 


+1 


+l 


+] 


+l 




+l 


-H 


+l 


■H 


+l 


-H 


+i 


+! 


■H 


■H 


+l 


+l 


+1 


r- 

Os 


co 
co 
O 


CN 


OS 

Tt 


IN 
IN 
O 


cn 
Tt 

CO 


m 

OS 
CO 


CO 
IT; 
IT; 


so 

CO 

m 


cn 
c 
cn 


r- 

Os 
f- 


m 

cn 


cn 
Tt 
t~ 


m 
m 


Tt 

CN 

o 


cn 

SD 
IT, 


o 

r- 


CN 
OS 

in 


cn 
cn 

CO 


cn 
cn 

CO 


m 
o 

CO 


m 

s 


m 
Tt 


CN 
SO 


in 

CN 
SO 


CN 

so 


Os 


Os 
CO 


CO 

o 

CN 


CO 


r~ 
o 

cn 




CN 

cn 


CN 

m 
cn 


cn 

CN 

in 


t— 

CO 
SO 


CO 
CN 

o 


m 

r- 

so 


m 
cn 


C 
m, 
c 


CN 
CN 
Os 


cn 

CN 
CN 


Os 

m 


SO 
t- 

o 


CO 
CO 


so 


r-^ 


SO 


SO 


CO 


d 

CN 


so 


CO 


d 

CN 


OS 


co 


OS 


d 

CN 


CN 


CO 


r~ 


d 

CN 


CN 


t~ 


d 

CN 


CO 


Os 


CN 


CN 


o< 


so 


d 

CN 


d 

CN 


o< 




so 


d 

CN 


d 

CN 


o< 


CN 


SO 


o< 


d 

CN 


os 


Os 


SO 


CN 


sd 


CO 


d 

CN 


o 
o 


o 
c 


o 
o 
o 


c 
o 
c 


o 
o 


CN 

c 
c 


O 

o 
o 


o 

C 


CN 

o 
o 


o 
c 


o 

O 


o 
c 


CN 
O 
O 


cn 
c 
c 


o 
o 


c 
c 
c 


CN 
O 
O 


Tt 

c 
c 


O 

o 
o 


cn 
c 
c 


CN 
O 
O 


CN 

c 
c 


cn 
o 
o 


in 

c 
c 


CN 

o 
o 


C 
C 
C 


cn 
o 
o 


O 

c 


o 
o 


c 

C 
C 


o 
o 


■* 

c 
c 


o 
o 


o 
c 


m 
o 
o 


c 
c 
c 


o 
o 


O 

c 


o 
o 


o 
c 


o 
o 
o 


cn 
o 
c 


o 
o 
o 


O 
O 


CN 
C 

c 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


+l 


+1 


+i 


+1 


+l 


+1 


+i 


+1 


-H 


+1 


-H 


+1 


+l 


+1 


+l 


+l 


+l 


+1 


+i 


+1 


-H 


+1 


+l 


+1 


+1 


+1 


■H 


+] 


+1 


■H 


■H 


+l 


■H 


+l 


+l 


+i 


+l 


-H 


+1 


+l 


+1 


+l 


+l 


+l 


+1 


r- 

00 
CO 


cn 
Tt 


Tt 
m 

CO 


c 
Tt 

Os 


o 

SO 

m 


OS 

cn 

CO 


Tt 

Os 

r- 


Tl- 
OS 

r- 


Tt 
so 
Tt 


CN 
C 
Os 


so 

Tt 

CO 


SO 

r- 
o 


o 

cn 

CO 


SO 

cn 

CN 


SO 
m 


CO 

cn 


m 
Os 
r- 


c 
f- 
in 


SO 

s 


CN 
CN 


m 
m 
r- 


r~ 

In 

SO 


CO 

cn 
tt 


CN 

c 

CO 


Os 
CN 
SO 


t- 




CO 

■* 


CO 

cn 
■* 


m 

Os 


o 

s 


cn 
■* 
so 


o 
o 

■* 


cn 
c 


m 

OS 


so 
c 


o 
o 
■* 


CN 
-J 
SO 


cn 


■* 
m 


o 


Os 


o 
o 

CO 


t— 

SO 


SO 

cn 

Os 




co 


so 


SO 


Os 


d 

CN 




CO 


d 

CN 


Os 


OS 


d 

CN 


d 

CN 


CN 


CO 


CO 


d 

CN 


CN 


CO 


CN 


d 

CN 


d 

CN 


CN 


CN 


d 

CN 


SO 


CN 


d 

CN 


d 

CN 


r~ 


SO 


CN 


d 

CN 


d 

CN 


CN 


CO 


d 

CN 


d 

CN 


Os 


d 

CN 




CN 


SO 


Os 


d 

CN 



CO 

- CD 

Q Q 



^ Cl> 

Pi Q 







o 


c 




CN 






CN 


CN 


CN 


CN 


CN 


cn 






CN 


Tt 


o 


Tt 


t~ 


Tt 


m 


m 


cn 


c 


CN 


CN 


cn 


c 


o 


c 




CN 


cn 


c 


cn 






cn 


o 


cn 


o 




CN 


o 


o 


o 


c 


O 


c 


o 


o 


O 


C 


o 


c 


O 


c 


o 


O 


O 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 




o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


O 


C 


o 


c 


o 


c 


O 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


C 


O 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


c 


o 


O 


c 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


+i 


+1 


+l 


+l 


+1 


+1 


+l 


+1 


+i 


+1 


+l 


+1 


+i 


+1 


+l 


+1 


+l 


+1 


+i 


+1 


+i 


+! 


+i 


+1 


+l 


+1 


+l 


+1 


■H 


■H 


+i 


+i 


■H 


+! 


+l 


-H 


+l 


+l 


+i 


+1 


■H 


+l 


+i 


+l 


+1 




c 


m, 


cn 


cn 


Os 


CO 


CN 


c^J 


Tl- 


Os 


CN 


SO 


so 


CN 


CO 


m, 


CO 


Tl- 


Os 


cn 


cn 


CO 


c 


CO 


cn 


o 


SO 


o 


cn 


TT 


cn 


Os 


OS 


m 




Os 




Os 


c 




m 


CN 


CN 


Os 


cn 




Os 




CN 


f- 




c 




Os 


CN 


m 


m, 


Os 


cn 


s 


m, 


OS 


OS 






m 




CN 


cn 


m 


Os 


Os 




m 






m 


cn 


CO 


CO 


o 




m 


Tt 


so 


Os 


SO 








m 


t- 




CO 


Os 


s 


■* 


r- 


Os 




CN 


q 


■* 


m 




q 


r-; 


m, 


o 




OS 


o 




CO 






q 




Tt 


cn 


Os 


^o 


cn 


CO 


CN 




Os 




CO 


m 


V~l 


cn 


CO 




Os 


CO 


t~ 




o 


d 


Os 


OS 


d 


d 










CO 


OS 






CO 


CN 


CN 




CN 


CN 




r~ 








CO 


r~ 


CN 


d 






OS 




d 


Os 




r~ 




r~ 


d 












CN 


CN 






CN 


CN 


CN 


CN 


CN 


CN 






CN 


CN 




CN 


CN 


CN 


CN 


CN 


CN 




CN 


CN 


CN 






CN 


CN 


CN 


CN 




CN 


CN 




CN 




CN 




CN 


CN 


Os 


CN 


cn 


CN 


SO 


m 


CO 


SO 






r- 


SO 


r- 


OS 


CN 


OS 


cn 


c 


cn 


Tt 


CN 


cn 


cn 


SO 


cn 


CN 


OS 


m 


o 


CN 


CN 


CO 


CO 


cn 


CO 


CO 


Tt 




m 


CO 








Os 




o 


C 


O 


c 


■* 




o 


c 




SO 


m 


m 






O 


C 




CN 


o 


SO 


SO 


c 


m 


cn 


Os 


c 


o 


CN 


CO 


c 


o 




o 


r- 




C 




c 


o 


CO 


o 




o 


CN 




o 


c 


O 


c 


o 


o 


O 


c 


o 


C 


o 


c 


o 


O 


o 


c 


o 


C 


o 


c 


ri 




o 


c 


o 


c 


o 


O 


O 


c 


o 


Tt 


o 


c 


o 


C 




c 


o 


c 


o 


c 


o 


o 


o 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


+i 


+1 


+1 


+1 


+i 


+l 


+l 


+1 


+l 


+1 


+l 


+1 


+l 


+1 


+1 


+1 


+l 


+1 


+l 


+1 


+l 


+l 


+i 


+! 


+1 


+1 


+l 


-H 


+l 


+1 


■H 


+1 


+1 


■H 


+1 


+l 


■H 


+l 


+l 


+1 


+1 


-H 


+l 


+l 


+1 


cn 


so 


m, 


CO 


o 


Tl- 


CO 




CO 


r- 


CN 


Tf 




SO 


'sO 


r~ 




Os 




c 




Os 




C 


CN 


o 


Os 




TT 


SO 




SO 


CO 


m 


cn 


CN 




CO 


Os 


c 


m 




so 


CO 


CN 


r- 


c 


o 


CN 


o 


OS 


CO 


CO 


CN 


c 


CN 


CN 




OS 


CN 


r-~ 


m, 


r-~ 


CN 


Os 


r- 


CO 


cn 




cn 




m 


m, 


o 




cn 


cn 


o 


cn 


cn 


OS 


Tt 


CO 


m 


Tt 


s 




Os 


CO 


C 


CN 


cn 


r- 




o 


■* 






cn 


Os 


cn 


CN 




C 


Tt- 


CN 


t-; 


Os 


CO 


cn 


so 


C 


CN 


CO 


Os 


s 


cn 


cn 


CO 


Tt 


Os 


OS 


cn 


m 


CN 


cn 




q 


TT 


O^ 




CN 


so 


t~ 


Os 




Os 


Os 


co 


cn 






d 




CN 


cn 


cn 




CN 


OS 








OS 


cn 


rn 


Tt 


cn 


CN 


cn 


CO 




CN 


cn 


OS 


CO 


Tt 




cn 


CN 




Tt 




d 


cn 


CO 


CN 


CO 


CN 




CN 








CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 




CN 


CN 


CN 




CN 


CN 


CN 


CN 


CN 


CN 




CN 


CN 


CN 






CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 




CN 




CN 


CN 


CO 


CN 


o 




O 


cn 


m, 


m 


O 


so 


CO 




Os 


Tl- 


Os 


c 


OS 


m 


O 


c 


'sO 


SO 


cn 


Tt 


cn 




SO 


SO 


CN 


Tt 


SO 


Tt 




m 


in 


CN 


SO 


SO 


m 


cn 




m 


in 


Tt 


m 


CO 


SO 


SO 


CO 


SO 


CN 


CN 


CN 


SO 


c 


■* 


m 


r» 


co 


CN 


Tl- 


t- 


Tt 


Tt 


SO 


CN 


CO 


m 




m 


m, 


m 


c 


TT 


CN 


CN 


cn 




c 


m 


m 


OS 


cn 


CO 


c 


CN 




m 


r- 


Os 


SO 


SO 


c^J 


CN 


■t 


so 




cn 


■t 


cn 




CN 


<N 


CN 


CO 


OS 


o 


c 


t- 


OS 


CO 


CO 


Os 


Os 


CO 


m 


CO 




SO 


Tt 


o 


in 


m 


Tt 


Tt 


Os 








c 


SO 


c 


m 


so 


CO 


OS 


OS 


Os 


OS 


Os 


Os 


Os 


OS 


Os 


Os 


Os 


OS 


Os 


Os 


CO 


CO 


Os 


Os 


co 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 




CO 


CO 


CO 


CO 




CO 




t- 


r~ 


OS 


OS 


Os 


OS 


Os 


Os 


Os 


OS 


Os 


Os 


OS 


OS 


Os 


Os 


OS 


OS 


Os 


Os 


OS 


OS 


OS 


OS 


Os 


Os 


Os 


OS 


Os 


Os 


Os 


OS 


OS 


OS 


Os 


OS 


OS 


Os 


Os 


OS 


OS 


Os 


Os 


OS 


Os 


Os 


OS 


■* 

1 


1 


it 
1 


1 


it 
1 


TT 

1 


Tf 
1 


Tt-' 
1 


it 
1 


TT 

1 


■* 

1 


TT 
1 


■* 

1 


TT 

1 


Tt-' 
1 


TT 

1 


Tt 
1 


Tt 

1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 

1 


Tt 
1 


Tt 
1 


Tt 

1 


Tt 
1 


Tt 

1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 

1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


Tt 
1 


o 


m, 


Vi 


SO 


O 


C 


m 


CN 




cn 


Vi 


CN 


CN 


so 


cn 


Tt- 


cn 


cn 




CN 


t- 


CN 


CO 




so 






CN 


CO 


r» 


m 


CN 


CN 


CO 


o 




o 


Tt 


SO 


cn 






in 




in 




f- 


CN 


CO 


m, 


Os 


cn 


c 


m, 


SO 


'SO 


m 






t- 


r~ 






m, 


c 


cn 


OS 


Os 


CN 


o 


m, 


m 


CO 


CN 


cn 




CN 


SO 




OS 


o 


so 






Tt 






in 


CN 


cn 


m, 


CN 


cn 


cn 


SO 


so 






O 




CO 


m 


O 


CO 


Os 




O 


so 


CO 


cn 


cn 


C 


OS 


m, 


OS 


CN 


SO 


c 


Os 


c 






SO 


Tt 




CN 


CO 


so 


TT 


c 


o 




TT 


in 


so 


CO 


cn 


os 


r~ 


it 


CO 


so 


s 


** 


OS 


m, 


TI- 


Os 


m, 


Tt 


m 


Tt 


m, 


m, 


so 


r- 


SO 


cn 




CN 


c 


o 


CO 


TT 




CO 


CO 


CN 




CN 




r~- 


CO 


co 


m 




m, 


o 


m 


c 


cn 


m 


o 




m, 


m 


CN 


cn 


CO 


m 


CN 


CS 


CO 


so 




CN 


m, 


Os 






CN 


CN 


CO 


cn 


m 


SO 


m 


OS 


o 


cn 


m 


OS 


m 


CO 


TT 


cn 


Os 


Tt 


Os 




CO 


CO 




Tt 


cn 


so 


sd 


SO 


sd 


SO 


SO 


so 


so 


SO 


sd 


so 


sd 


SO 


sd 


SO 


SO 


so 


sd 


so 


so 


SO 


SO 


sd 


so 


sd 


SO 


sd 


SO 


sd 


so 


sd 


SO 


sd 


so 


sd 


so 


SO 


so 


SO 


so 


^d 


so 


sd 


so 


sd 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 



48 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 



E.2. MegaCam transmission functions 

Low resolution version of the Megacam effective transmission 
curves are tabulated in tables E.3-E.8. Higher resolution versions 
are available on the joint calibration webpage. 

As described in Sect. 3.1, the transmission functions depends 
on the location in the MegaCam focal plane. Assuming that the 
transmission vary continuously with the radius from the focal 
plane center provides an effective description. We deliver tab- 
ulated measurements at regularly spaced radii. Those measure- 
ments can be linearly interpolated to other positions. 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 
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E.3. HST standard measured magnitudes 

The tertiary catalogs define a magnitude system. This system 
has been anchored to the AB flux scale according to a calibra- 
tion reference. We choose to use, as our reference, 5 stars from 
the HST CALSPEC database for which direct or indirect mag- 
nitudes in the system defined by the tertiary catalogs were avail- 
able. Therefore, our calibration reference is the current release 
(stinic_003) of the CALSPEC spectra for those 5 stars. 

The natural magnitudes of the measured HST standards are 
given in Table E.9. Those measurements can be used to up- 
date the AB calibration of the catalogs, given revised synthetic 
magnitudes for these stars. Measurements in table E.9 are deliv- 
ered in the same photometric system than the published tertiary 
standards catalog, i.e. with the AB offsets derived in this work 
(Table 23) already applied. Therefore, any offset computed, ac- 
cording to Eq. (31), from the difference between a measurement 
given in Table E.9 and the synthetic AB magnitude for the cor- 
responding star readily applies to the published catalog. 

When computing synthetic AB magnitudes for a star in the 
MegaCam focal plane, the peculiar position of the star in the 
focal plane must be taken into account. Also, we describe in 
Sect. E.4 the constitution of the covariance matrix for measure- 
ments in Table E.9. This matrix can be used to combine all mea- 
surements similarly to what has been done in Sect. 9. 
Table E.9. Measured magnitudes for CALSPEC standards. 
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" Xfp and jfp are the coordinates of measurements in the MegaCam 
focalplane (in cm, with the origin taken at the center of the focalplane). 



Instrument Band Star Xf p yt v " Mag 



" and y fp are the coordinates of measurements in the MegaCam 
focalplane (in cm, with the origin taken at the center of the focalplane). 

E.4. Covariance matrices 

The assumptions made to build the covariance matrix of system- 
atic measurement errors in the calibration data follows. 

In SNLS, we assume that the shutter bias affects identically 
all measurements of flux ratios between short and long expo- 



1\/T f (m in m 

lVlClid.v^d.111 


I'M 


SNAP? 


-1 04. 


-1 Ql 


16 O^S 


lVTf* trn in m 


I'M 


SNAP? 


-1 04. 


-? Ql 


16 01 8 

1U.VJ1 o 


1 VIC Id d.v^d.111 


I'm 




-0 78 
w. / o 


1 1 


1 ? 80S 


lvAfatTii i n m 

1 VIC Id d.v^d.111 


I'm 


P1 77D 


-0 59 




1 ^ 767 


lVlCldd.v^d.111 


I'm 


SNAP? 


-1 0^ 




1 S Q8Q 


lVlClid.v^d.111 


I'm 


SNAP? 


-1 0^ 


-? Ql 


1 S QQQ 

u . yyy 


lVTf* trn in m 

lVlCiid.v^d.111 


I'm 




-0 78 


1 1 


1? 81 S 


lVT^trn in m 

lVlClid.V^d.111 


I'm 


P1 77D 


-0 S1 

W.J 1 




1 ^ ?7Q 

1 J .Z, / 7 


lVlCldd.v^d.111 






-0 78 
w. / o 


1 1 


1 ? 688 


iVlCldd.v^d.111 




P1 77n 


-0 59 


J .'-tJi 


1^ 160 


lVlCldd.v^d.111 




SNAP? 


-1 04 


-? Ql 


1 S 874 


A/Tf^trn in m 

lVlClid.v^d.111 




SNAP? 


-1 04 


-? Ql 


1 S 8Q0 


A/If* tin in m 

lVlClid.V^d.111 






-0 78 


1 1 


1 ? 686 


A/If* tin in m 


Or, 


P1 77D 


-0 S9 






A/T poji in m 

lVlCldd.v^d.111 


Or, 


SNAP? 


-1 0^ 


-? Ql 


1 S 88Q 

IJ.OO" 


A/I pna i n m 

lVlCldd.v^d.111 


Or, 


SNAP? 


-1 0^ 


-? Ql 


1 S 846 


lVlCldd.v^d.111 


Or, 




-0 78 
w. / o 


1 1 


1 ? 70^ 

1Z. / UJ 


\/| fnyn in m 

lVlClid.v^d.111 


Or, 


P1 77D 


-0 S9 


J .'-tJi 


1^ 161 


A/T trn in m 

1 vie id d.v^d.111 


Or, 


SNAP? 


-1 0^ 


-1 Ql 


1 S 8QS 


A/If* trn in m 

lVlCldd.v^d.111 


Or, 


SNAP? 


-1 0^ 


-? Ql 


1 S QO^ 


A/Tf^crn in m 

lVlCldd.v^d.111 


Or, 




-0 7Q 


1 1 


1 ? 6QQ 


lVTpcrnf^n m 

lVlClddV^dlll 


Or, 


P177D 


-0 55 




n i6Q 


£i(Tn i n m 

IVlCldclv^cllll 


Or, 


SNAP? 


-1 OA 


-? Ql 


1 S 8QQ 
i j . oyy 


A/T tin in m 

lVlCldd.v^d.111 


Or, 


SNAP? 


-1 OA 


-? Ql 


1 S 868 
1J, ouo 


A/T f 2 * trn in m 

lVlCldd.V^d.111 


Zm 




v/. / O 


1 1 


1 ? 6QS 


A/T tin in m 

lVlCldd.v^d.111 


Zm 


pi 77n 

I III \_J 




J • ii 


n 140 


lVlCldd.v^d.111 


Zm 


SNAP? 


-1 0^ 

1 .WJ 


-? Ql 


1 S 874 


lVlCldd.v^d.111 


Zm 




w. / o 


1 1 


1 ? 676 

1Z.U / u 


A/1 poji i n m 

lVlCldd.v^d.111 


Zm 


P1 77n 

XI// LJ 


U.JZ 




1^ 1^0 


A/If* trn in m 

lVlCldd.v^d.111 


Zm 


SNAP? 


-1 fH 


-? Ql 


1 S 877 

LJ. OZ / 


A/T tin in m 

lVlCldd.v^d.111 


Zm 


p^of 


-0 7Q 


1 1 


1 1 67Q 


A /I £*rrn l 1 n m 
ivicyd.v^d.111 


Zm 


pi 77D 


S? 

-U.JZ 


^ zLzt 

- J .H-H- 


1 ^ 1 ^6 


A/1 poa in m 

lVlCldd.v^d.111 


Zm 


SNAP? 


1 .\J'-r 


-? Ql 


1 S Q00 




it 








1^ SOI 

1 J.JU 1 




8 


ru4iL 






1 ? 781 




F 








1 1 84Q 


OL'kJOZ.J 


I 








1 1 7^4 
ii./ jt 




z 








1 1 717 




11 








14 48S 




8 








1 ^ ^00 
1J. jUU 




V 


p^of 






1 7 844 


OJ.-'OOZ,. J 


I 








17 71 S 




z 








1 7 68 S 




11 


pi 77n 

X X / ID 






1 S OSO 

1 J . uJ W 




§ 


P1 77D 

XX// 






1 ^ 76^ 

1J. / UJ 




V 


P1 77D 

X X / ID 






1 ^ ^04 

1 J . JUT" 




I 


P177D 

I III D 






1^ 17S 

ij.i / j 


SDSS2.5 


z 


P177D 






13.138 


SDSS2.5 


u 


BD17 






10.492 


SDSS2.5 


8 


BD17 






9.651 


SDSS2.5 


r 


BD17 






9.357 


SDSS2.5 


i 


BD17 






9.263 


SDSS2.5 


z 


BD17 






9.251 



56 



M. Betoule et al.: Photometric Calibration of the SNLS and the SDSS. 



sures. The error resulting from aperture corrections is assumed 
to affect coherently the measurements of the solar analogs in a 
given band, and without correlation from one band to another. 
We account for the aperture correction residuals in the Landolt- 
based BD +17 4708 measurements independently as the ob- 
servation strategy was different and the derivations of the cor- 
rections were independent. The survey non-uniformity affects 
coherently the measurement of solar analogs that are related 
to the survey through the SNLS D3 field. The Landolt non- 
uniformity is included to account for a potential departure of 
the BD +17 4708 measurements from the average Landolt sys- 
tem. The errors on the color transformation of BD +17 4708 are 
assumed to be independent from one band to another. This may 
not be fully justified, but this error matters only in band z and 
potential correlations with other bands are negligible. 

In SDSS, the statistical uncertainty on the color transforma- 
tion between the PT and 2.5 m system affects coherently all the 
measurements in a given band. Star-to-star dispersion around the 
color transformation is assumed independent from one star to an- 
other. The non uniformity of the Smith catalog coherently affects 
the measurement of solar analogs as they lie at approximately the 
same hour angle. 

The resulting covariance matrix of measurements in 
Table E.9 (including measurement noise and systematics) is 
available from the SNLS/SDSS calibration webpage. 

As described in Sect. 9, when combining the available mea- 
surements to provide the final calibration, we account for other 
contributions to the covariance of the AB offsets: the HST-STIS 
spectrum measurement error, the STIS calibration error, and the 
uncertainty on instrument filters. A corresponding covariance 
matrix for each of these contributions is also available from the 
SNLS/SDSS calibration webpage. 

The estimation of the STIS measurement error is discussed 
in Sect. 9.2. The corresponding covariance matrix for a single 
epoch STIS measurement is displayed in Table E.10. 
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Table E.10. Covariance matrix of synthetic STIS magnitudes. 



cov(m,m') = 10 



25.8 17.2 9.8 -2.9 -2.1 -16.4 26.3 18.7 10.2 -1.4 -16.6 \ u M 

15.5 3.9 -4.4 -3.7 -15.6 16.8 16.3 4.2 -3.2 -15.9 g M 

26.0 21.0 20.8 26.0 9.9 3.8 26.6 20.8 26.9 r M 

35.5 33.5 58.6 -3.6 -5.0 21.1 32.9 60.2 i M 

31.8 54.7 -2.7 -4.4 20.9 31.2 56.2 U M 

119.1 -17.6 -17.2 26.1 52.6 122.5 z M 

27.5 18.3 10.3 -2.0 -17.9 u 

17.3 4.0 -3.8 -17.5 g 

27.3 20.9 26.9 r 

30.7 54.1 i 

125.9 ) z 



(E.l) 
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